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Further results of an extensive study of trajectories asymptotic to a known family of unstable 
periodic orbits in the earth’s dipolar magnetic field, carried out by means of Bush's differential 


analyzer, are presented in this paper. A detailed discussion is given of our methods of deter 


mining asymptotic trajectories by means of the differential analyzer and by numerical integra 


tion of a whole family at a time; comparison of the results obtained shows the absence of 


systematic errors of any consequence in the mechanical integrations and exhibits the precision 
attained with the differential analyzer. The families of asymptotic trajectories are then analyzed 
systematically in order to determine the main cones for latitudes up to 30°. This leads to the 
theory of the azimuthal effect and a study of the region in the vicinity of the zenith. 


N a preceding paper' to which reference should 

be made for a complete statement of the 
problem treated here and of our methods of 
attack, we obtained from the 
analysis of hundred asymptotic 
trajectories to a family of unstable 
periodic found by Bush's 
differential analyzer? and discussed fully the 
sections of the main allowed cones of cosmic 
radiation by the meridian plane, which in turn 
led us to the theory of the north-south asym- 
metry. Those results were presented at the time 
with reservations as far as a critical examination 
of their precision was concerned. We have now 
been able to complete the calculations announced 
in our preceding paper, to which we shall return 
below, with the result that suspected systematic 
errors are so small that they can well be neg- 
lected. Certain other points which were then 
summarily sketched will now be developed in 


gave results 


some three 
known 


orbits means of 


1G. Lemaitre and M. S. Vallarta, Phys. Rev. 49, 719 
(1936). 
*V. Bush, J. Frank. Inst. 212, 447 (1931). 


detail; in particular we shall present here a full 
discussion of our method of determination of 
asymptotic trajectories by means of the differ- 
ential analyzer and a fairly complete determina- 
tion of the main allowed cones for geomagnetic 
latitudes up to 30°. The last two sections will be 
devoted to the analysis of the azimuthal effect 
and to a study of the region in the vicinity of 
the zenith. 

DETERMINATION OF AsyMpToTic TRA- 
JECTORIES BY MEANS OF BusuH's 

DIFFERENTIAL ANALYZER 


1. THE 


The differential equations of motion to be 
integrated are (reference 1, Eqs. (4), (5)) 
d?x/do? = (1/(2y,)*)e?* —e~7+e-** cos? X, 
(1) 
d*\/do? =e-** sin \ cos A— (sin X/cos* A). 


Our problem is to find the solutions of these 
corresponding to tra- 
known family of 


differential equations 
jectories asymptotic to a 


unstable periodic orbits. 
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Before asymptotic trajectories can be studied, 
a full theory of the family of periodic orbits 
mentioned above must be available. The family 
in question was discovered by Stérmer*® who was 
also able to calculate two of its members by his 
method ef numerical integration. A complete 
treatment of the family has been given by one 
of us* so that any one of its members can be 
calculated without great trouble and with high 
precision. While the details of these calculations 
would be out of place here and consequently are 
left for a separate paper, the results are repro- 
duced in Fig. 1 where the small circles are the 
points on the periodic orbits of phase 0°, 15° 
30°, --- 90° which were computed. The equi- 
phase line 90°, i.e., the locus of the vertices of 
the periodic orbits which is important for our 
purpose as further discussed in a later section 
of this paper, has also been drawn in the figure. 

Solutions of the system (1) representing the 
desired asymptotic trajectories and valid in the 
vicinity of the periodic orbits have been given 
by Bouckaert;> other methods which we have 
recently developed and to which we return 
below are now available. In order to extend 
these solutions to the region of low energies and 
high latitudes, however, recourse must be had 
in general either to methods of numerical inte- 
gration such as the one outlined in a later section 
of this paper or to methods of mechanical 
integration such as that developed by Bush,’ 
because in general the system (1) cannot be 
integrated in terms of known functions. Since 
the knowledge of several hundred asymptotic 
trajectories is required for a satisfactory treat- 
ment of our problem the first line of attack 
involves a very large amount of tedious labor 
while the second can still be carried through in 
much shorter time. 

An outline of our method of determination of 
asymptotic trajectories by means of Bush’s 
differential analyzer has already been published.! 
Here we add only such details as are needed for 
a full understanding of our line of attack and of 
the results obtained. A first integration of the 
system (1) yields 


3 C, Stérmer, Zeits. f. Astrophys. 1, 237 (1930). 

4G. Lemaitre, Ann. de la Soc. Sci. de Bruxelles A54, 194 
(1935). 

5 L. Bouckaert, Ann. de la Soc. Sci. de Bruxelles A54, 174 
(1935). 
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Fic. 1. The family of periodic orbits. 


dx/da= S ((e**/16y1*) —e*+ cos? A)e**da, 
dd/do= S (e** — (1/cos* d)) sin A cos Ado. 


When written in the above form, suggested by 
S. H. Caldwell, the system (1) may be immedi- 
ately set in the differential analyzer now available 
at the Massachusetts Institute of Technology 
and solved mechanically. This differential ana- 
lyzer has ordinarily four input tables and six 
integrators, but an ingenious change in the 
machine connections devised for the present 
problem by S. H. Caldwell enables the use of 
the output table simultaneously as an input 
table so that in all five input tables and six 
integrators were available. Of these five input 
tables, two are needed to introduce in the 
machine the functions e* /16y;4—e? and e~**, and 
three more to introduce cos? \, 1 ‘cos* \ and sin \ 
cos \. Of the six integrators, two are required to 
integrate the product of the two functions under 
the integral sign in each of the two equations (2), 
and two more to integrate the derivatives dx/do 
and di/de so as to obtain x and X. The output 
table is then controlled so as to plot x as abscissa 
and \ as ordinate. 

The scale factors are determined from the 
range of x and \ which it is desired to explore. 
For the present investigation the scales were so 
chosen that x could vary from —1.0 to 0.7, \ 
from —0.7 to 1.0, dx/do and d\/do from —0.5 
to 0.5. The limits of the last two functions were 
found from a study of the energy integral 
(reference 1, Eq. (5)). 
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TABLE I. Data for some of the fundamental points. 

yi= 1 0.87 0.90 0.93 0.96 

x Mo o ” n v 
0.075 3.1587 1.166 1.022 0.850 . 
0.100 3.2127 - - - — 
0.125 3.2829 ‘ 
0.150 3.3670 1.211 1.062 0.890 
0.175 3.4647 
0.200 3.5790 0.721 
0.225 3.7114 1.281 1.133 0.959 ).738 
0.250 3.8608 
0.275 4.0352 
0.300 4.2342 1.371 1.266 1.051 0.821 
TABLE II. Values of the initial angle for y,;=0.91 for all 


equatorial runs. 


x v “ n 
0.075 —0.753 3.2265 0.968 
0.100 —0.531 3.2605 0.978 
0.125 —0.252 3.3056 0.992 
0.150 0.084 3.359 1.008 
0.175 0.380 3.431 1.029 
0.200 0.760 3.511 1.053 
0.225 1.276 3.596 1.079 
0.250 1.914 3.689 1.107 
0.275 2.745 3.788 1.136 
0.300 3.650 3.907 1.172 


The initial conditions are introduced into the 
machine knowing x, A and the initial slope tan 7. 
The knowledge of x and ) is sufficient to set the 
starting points on the input tables. To set the 
starting points on each one of the six integrators 
it is necessary to calculate the values of the 
functions e~**, 


e**/167:4—e7+ cos? A, e°*—1/cos* X, 


sinXcos\, dx/de and dd/de. 


The first four are computed from the known 
values of x and X, the last two from the initial 
inclination » by the relations 


dx/da=P'cosy; dd/do=P'sinn, (3) 
where P is given by Eq. (5).! 

For the efficient operation of the machine it 
is important that the calculations described 
above be carried out rapidly and accurately. 
For this purpose master sheets were prepared in 
advance giving the value of the required func- 
tions from x= —0.125 to x=0.350 by intervals 
of 0.025 along the equator where P reduces to 
e**/16y;4—(e-*—1)* and for every degree of 
latitude along the line 6=0 where P=e**/16y;'. 
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As already mentioned in our preceding paper,' 
for convenience trajectories were always started 
either from the equator or from the line @=0. 
Once the master sheets are ready it is easy to 
calculate the initial conditions while the machine 
is tracing a trajectory, i.e., in the interval of a 
few minutes. 

As already mentioned in our preceding paper 
one method for the determination of asymptotic 
trajectories consists in choosing a point on the 
equator within the periodic orbit and the initial 
slope arbitrarily. The initial conditions having 
been set on the machine, a trajectory is started 
in the direction towards the periodic orbit; the 
initial slope is then adjusted until the trajectory 
neither intersects nor falls short of the periodic 
orbit. More precisely let us suppose that the 
initial slope has been estimated not too far from 
the correct value, and assume tor the sake of 
exemplification that the trajectory cuts through 
the periodic orbit. The initial angle is then 
increased by equal amounts, say of 0.008 radian 
until the trajectory falls decidedly short of the 
periodic orbit. The angle is then decreased by 
steps of 0.004 radian until the trajectory again 
cuts through, then increased by steps of 0.002 
until the latter falls short, then again decreased 
by steps of 0.001 until it cuts through. In this 
way it is usually possible to determine the 
critical angle for an asymptotic orbit within 
around 0.001 radian by making from five to 
ten trials, each one taking from four to six 
minutes. That the critical angle is actually 
determined by this method with a precision of a 
few thousandths of a radian is confirmed by 
independent calculations of asymptotic tra- 
jectories to be more fully discussed in a later 
section of this paper. 

In order to be free from the need of finding 
the critical angle for every equatorial point in 
the manner outlined in the preceding paragraph 
a method of interpolation was devised which 
permits the calculation of critical angles at 
intermediate points once the values at certain 
fundamental points have been ascertained. Let 
e=1—y; and let 


Ho=lim (n/eé), (4) 


v1 


n being the initial critical angle for a given value 
of x. The quantity uo may be calculated from 





496 G. LEMAITRE AND 


TABLE III. Initial conditions for various values of y; and 
values of x <0.075 along the equator. 


y1 0.87 0.89 0.91 0.93 0.95 
=X n n n n n 
—0.100 . 0.780 
—0.075 0.890 0.747 
—0.050 0.981 0.866 0.721 
—0.025 0.965 0.852 0.707 
0 1.160 1.050 0.957 0.840 0.703 
0.025 1.155 1.050 0.956 0.840 0.705 
0.050 1.163 1.055 0.960 0.842 0.711 


Bouckaert’s formulas’ for values of x greater 
than about 0.1 and is given in Table I above. 
Let now uw be the corresponding quantity for 
7¥1~%1. We have 

p=n/e}, (5) 


where 7 is now determined by the trajectories 
traced by the differential analyzer as outlined 
above. Let us place 


M>=pMo— Ve. (6) 


Our process of interpolation consists in finding 
» for certain fundamental points and then 
calculating v from it by the use of the formulas 
written above ; v is then plotted to an appropriate 
scale as a function of y; for each chosen value of 
x through which we desire to have an asymptotic 
trajectory. For interpolated points vy is deter- 
mined from these plots knowing x and y, u is 
then calculated and finally ». Table I gives the 
data for some of the fundamental points. 

As an example of the use of the method of 
interpolation we give in Table II the values of 
the initial angle for y,=0.91 for all equatorial 
runs. 

Similar tables were prepared for all other 
values of y, which were studied (reference 1, 
p. 722). In all over one hundred asymptotic 
trajectories starting from equatorial points were 
thus determined. In order to verify interpolated 
values a few trajectories were continued as far 
as the periodic orbit with very satisfactory 
results. That the precision of interpolated tra- 
jectories is as high as that of trajectories through 
fundamental points will be shown in the next 


section. 

For points on the equator for which yo cannot 
be calculated from Bouckaert’s formulas, i.e., 
for very small and for negative values of x, the 
initial angle was determined from the condition 
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that an asymptotic trajectory must be tangent 
to the envelope of the asymptotic family which 
was known from those asymptotic trajectories 
already found. As a further verification the 
trajectories were continued in some cases as far 
as the periodic orbit, but because of the sharpness 
of the turns when such asymptotic orbits oscillate 
in the vicinity of the periodic orbit the ordinary 
condition for an asymptotic trajectory cannot 
be ascertained with precision by means of the 
differential analyzer. As an example we give in 
Table III the initial conditions for a few of the 
values of y; we have studied and values of 
x <0.075 along the equator. 

As regards asymptotic trajectories starting 
from the line 6=0 the procedure was similar to 
that described above: the criterion determining 
an asymptotic trajectory was that it must be 
tangent to the envelope of the asymptotic family. 
The interpolation was made directly on the 
initial angle, the procedure being to plot to an 
appropriate scale 7 as a function of \ for each 
value of y, from a few fundamental points and 
then to read off from the plot the value of » for 
intermediate values of \. Curves showing the 
relation between the critical angle and the 
latitude for all values of y; along the line @=0 
have already been published (Fig. 5, reference 1). 
As an example the values of the critical angle 
for fundamental points corresponding to y;=0.85 
are given in Table IV: 

Examples of families of asymptotic trajectories 
are given in Fig. 2. In this figure portions of 
trajectories intercepted by the earth are shown 
dotted, i.e., for these portions a parallel to the 
\-axis cuts a preceding part of the trajectory. 


TABLE IV. Values of the critical angle for fundamental points 


for y:=0.85. 
d (deg.) n m1 

— 23 2.170 
— 20 2.326 
0 1.253 1.885 

5 1.140 

10 1.035 
15 0.915 1.451 
20 0.794 1.235 
25 0.603 0.975 
26 0.562 0.915 
27 0.523 0.830 
28 0.515 0.750 
29 - 0.600 
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OF 


A few trajectories of the third kind and a self- 
reversing doubly asymptotic trajectory with 
(cf. sec. 4, 


self-reversal point 7> reference 1, 


Fig. 4) are clearly seen. 
TRAJECTORIES 


2. PRECISION OF ASYMPTOTIC 


outline of a 
method of calculation of asymptotic trajectories 


In this section we present an 
which we have developed recently in order to 
investigate systematic errors in the trajectories 
found by means of the differential analyzer. 
The details are reserved for a paper to appear 
shortly in the Annales de la Société Scientifique 
de Bruxelles. It 
errors in asymptotic trajectories determined by 


will be seen that systematic 
means of the differential analyzer are small 
enough that they can well be neglected, except 
for very high latitudes, and that the precision 
attained throughout is very satisfactory. 

The 


represented by a trigonometrical series in wr+ ¢, 


family of asymptotic trajectories is 


where 27/w is the period and ¢ the arbitrary 
phase of the periodic orbit. We have 


A 


x=D>[y:(e) sin k(wo+ ¢) 


mt 


+2,.(¢) cosk(wo+y)] (7) 
and a similar expression for \ with amplitudes 
u, v in which k& takes even integral values in the 
expansion for x and odd values in that for X. 
The asymptotic behavior of y(c) and 2(c¢) is 
exhibited by the expression for 2, : 


2. (0) =2.+2,/e% +2,/"e22o+..--, (8) 


2 - a o@ ; e v7 
where 2, 2%’, 2; 


are numerical coefficients of 
which the first corresponds to the periodic orbit. 
The variational equation then yields a set of 
linear and homogeneous equations the determi- 
nant of which must vanish. This condition, i.e., 
the determines the 
characteristic exponent 2. We have thus com- 


“secular equation,” real 
puted the first and second order terms in e®’. 
The second order term of 29 is rather large, so 
that a more convenient approximation had to 
be derived. This approximation is suggested by 
previous investigations of one of us and of 
Bouckaert® for the case where y; differs slightly 
from 1. Writing 


e zo(o) — e~70/ 1 +u), 





(9) 
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Fic. 2. A family of asymptotic trajectories. y; =0.93 
we have taken as approximation a solution of 


the equation 


(du/Qdo)?=u?(1+bu+cu?), (10) 


where the constants } and c are determined more 
or less empirically. By this means it was found 
possible to use the expansion up to a distance 
0.048 from the periodic orbit. Careful numerical 
the 
decimal. Difference tables of the functions y(c), 


checks proved it to be correct to sixth 
z(o) and the corresponding functions u(o), v(o) 
in the expansion for \ were then prepared for 
the six preceding values of ¢ with an interval 
corresponding to wio=15°, ie., Ao=0.293243. 
While the functions y, z, 
the corresponding trajectories perform many 


etc., vary very slowly 


very close oscillations in the vicinity of the 
periodic orbit so that their calculation by the 
usual methods of numerical integration would 
be extremely difficult, in any case requiring a 
much smaller interval Ac. For this reason instead 
of integrating numerically the differential equa- 
tions of motion in order to find the trajectories 
we have preferred to integrate the differential 
equations of the amplitudes y, 2, etc. 

The differential equations satisfied by these 
functions can be found as follows: replacing x 
and X by their series expressions and expanding 
in a trigonometrical series we may write 


40P/dx=>-[ Yi(c) sin k(wo+¢) 
k 


+Zi(c) cos k(wo+¢)]. (11) 


This defines Y and Z as a function of y, 2, » 
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Fic. 3. Calculated family of asymptotic trajectories. 
71 =0.929898. 





and v because of the equations of motion (1). 
For practical use of this equation we may take 
advantage of the fact that with the exception of 
Zo, wi and y; all the other terms are small and 
powers higher than the second may be neglected. 
In any case »; can be made to vanish by a 
suitable change of the phase angle. Once the 
phase has been shifted explicit expansion can 
be derived with coefficients depending on 2) and 
the new yw;’, and included in a convenient 
computation schedule. Finally the phase must 
be shifted back to its original value. 

The equations to be integrated are of the form 

(d?y,./do”) — 2kw(dz,/do) —kw*y,= Yx, 
(12) 

(d?z,. do*) +2kw(dy, ‘da) —kw*z,=Z,, 


and similar equations for pu, and ». Ordinary 
methods of integration such as Adams’ would 
have been impracticable as they would be 
tantamount to an attempt to calculate a sine 
from its second order differential equation using 
intervals of 15°, 30°, - 90° according to the 
value of k. Introducing the complex notation 
x=y+iz, X= ¥Y+iZ and taking the case k=1 
as typical, the above equations may be written 


(d?/da*)(xe'*) = Xe", (13) 
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which yields upon integration and elimination of 
integration constants by taking second differ- 


ences 

Aexei*t =A? ff Xe'*do’. (14) 
When w=0 we fall back on Stérmer’s method of 
numerical integration. The question now is to 
devise a suitable method of integration for the 
present equation. 

Writing wo=Ww(to+t), (15) 
where ¢ is the ordinal number of the integration 
points and w is the interval of integration, and 
supposing that the values of x,, X, are known 
up to n=—1, the question is to compute the 
value of xo as a function of x_,, x2, X_1, 6X_, 
= Y_,;—x_», &x_,= 6x_, — bx_», etc. The left-hand 
member of the equation 


Axe =u? ff Xe dt, (16) 


where a=iww is Xy»— 2e~*x_, +e7-**x_». X may be 
expressed as a function of Y_, and its successive 


differences by 


X (t) = D0 (1/k !)(t4+1)(¢-4+2)-- + (t+h)6*No. (17) 
k 


If D=0/da we have 
SS tre*'d? =D" ff e*'dt? = D"(e*'/ a?) (18) 
and the final relation is 
Xo = 2e~*x_1 — e774 x_2 + Dow?(1/k '!)(D+1)(D+42) 
k 
-++(D+k)(1—2e7-* +e? */ a?) 5k X_; 


== 2e~ "x. a sax ot > a,.5*X 1- (19) 
k 


Returning to real variables and writing a,=) 
+ic, we have finally the working formulas 


Vo= 2 COS ww: y_1—Y_2 COs 2ww 


+2s_1;sin ww—2_2 sin 2ww 


+w(>> b,6*V_1-—> c,6*Z_1) (20) 


k=0 k=0 


i) 


Z_1 COS wW—2Z_2 Cos 2ww 
—2y_; sin ww+y_.2 sin 2ww 


+w(>b.6Z 1+ c.6*+ V1). (21) 
k=0 k=0 





0. 
0. 
0.. 
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The coefficients };, ¢, have been computed up to 
k=4, i.e., up to and including the fourth differ- 
ence and for wi¢o=15°, 30°, --- 90° correspond- 
ing up to the sixth harmonic. 

The integration has been actually carried out 
for twenty steps until the region was finally 
reached where the change of the integrated 
functions became so large that the integration 
could no longer be continued without splitting 
the interval of integration. The family of twenty- 
four trajectories, symmetrical in pairs, has been 
computed and is shown in Fig. 3. For a table of 
the computed points reference must be made to 
the complete paper. The small circles shown in 
the figure are calculated points for which the 
phase is 0°, 15°, --- 90° and correspond with 
those shown on the periodic orbits (Fig. 1). 

The comparison with the asymptotic tra- 
the 
analyzer is shown in Fig. 4. It is seen that the 


jectories found by means of differential 
agreement is excellent particularly as regards 
the position of the cusp Co. One of the machine 
trajectories starting from x=0.275 and \=0° 
happens to have a value of g» equal to that of 
the calculated trajectory and both coincide, the 
others have different values of gp but that they 
all belong to the same family can be tested by 
analyzing their general shape and proved by 
the fact that they all are tangent to the same 
envelope. As a matter of fact the consistency 
between calculated and machine trajectories is 
the machine 
trajectories themselves. A more precise com- 


as good as that found among 
parison may be made by calculating the value 
of the inclination » of three fundamental tra- 
jectories at the equator. The result is shown in 


Table V. 

3. DETERMINATION OF THE ALLOWED MAIN CONE 
The theory of the main cone, or region of full 

light, has already been outlined in our previous 


papers. A sketch of our method for its determi- 
nation and its representation, once the families 


TaBLeE V. Calculated values of the inclination » of three 
fundamental trajectories at the equator. 

x nm (mach.) n (calc.) Diff. 
0.150 0.890 0.887 0.003 
0.225 0.959 0.953 0.006 
0.300 1.053 1.045 0.008 
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Fic. 4. Calculated asymptotic trajectories compared with 
trajectories obtained with differential analyzer. y; =0.93. 


of asymptotic trajectories have been found, has 
also been given in our preceding paper (reference 
1, p. 724). The coordinates of a representative 
point are sin@ and cos @sin 7; for positive 
particles 6 and 7 are counted positively eastwards 
and northwards, respectively. Let us assume an 
asymptotic trajectory passing through a point 
of coordinates x and \; 7 is then the angle with 
the zenith direction, i.e., with a parallel to the 
x axis and @ the angle between the trajectory 
and the meridian plane which can be immediately 
calculated from Stérmer’s formula (reference 1, 
Eq. (1)). The energy 7 in Stérmer’s corresponding 
to the given x is then given by Eq. (2).' We thus 
obtain a point of the representation of the cone 
for the given values of the energy and the 
latitude. 

For the systematic utilization of the asymp- 
totic the 
analyzer two lines of attack are open. Either 


trajectories traced by differential 
one may first calculate from Stérmer’s formula 
the lines 6=const. for all values of y; for which 
asymptotic trajectories were studied and then 
measure the values of the angle 7 along these 
lines noting at the same time the value of x. 
This is equivalent to determining sections of the 
cones of different energies and at different lati- 
tudes by planes parallel to the meridian plane. 
Suppose that the lines @=0°, 5°, 10°, --- have 
been computed for y;= 1.00, 0.99, 0.98, --- 0.78, 
and the angles 7 measured along each. Then 
plots may be readily constructed giving the 
angle » as a function \ for each value of 7 and, 
besides, the energy r as a function of x for each 
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value of y;. This is the method used for @=0° 
in our preceding paper! where diagrams giving 7 
as a function of \ for each value of y; are given 
(Fig. 5, reference 1). It has not been used for 
other values of @ because the computation of the 
lines @=const. is laborious. Once these have 
been found, however, the calculation of the cone 
is immediate. 

An alternative line of attack is to measure the 
inclination » along the lines \=const. for each 
value of y; which was studied, noting at the 
same time the value of x where the intersection 
between the asymptotic trajectory and the 
chosen line \=const. takes place. Graphs may 
then be drawn showing 7 as a function of r for 
each chosen value of A, for example \=0°, 5°, 
10° --+ and different values of y;.° An envelope 
is characterized on these nr-diagrams by a point 
of tangency of 7 as a function of r with a line 
parallel to the 7-axis, i.e., by a turning point for 
a maximum or minimum value of 7. A cusp is 
characterized by the coalescence of two such 
points; it is a point of inflection with a tangent 
parallel to the y-axis. These peculiarities of the 
nr-diagrams are reflected in the representation of 
the cone; their significance has already been 
taken up fully in our preceding paper.' If on the 
nr-diagrams a value r=const. is picked out it is 
readily seen that the boundary of the cone 
progresses eastwards (for positive particles) as 
yi increases. 

The main cones for \=0°, 20° and 30° are 
reproduced in Figs. 5, 6 and 7. The small circles 
are the points calculated by Bouckaert.® It is 
seen that except for high zenith angles and low 
energies the agreement is again quite satis- 
factory. The extreme northern part of the main 
cone where the shadow of the earth makes itself 
felt is given with reservations as a detailed 
analysis of the trajectories of the second and 
third kinds outside of the meridian plane still 


remains to be made. 


4. THe AZIMUTHAL EFFECT 


We begin by recalling that in addition to the 
main cones determined as described in the 





6 This laborious piece of work was carried out by Mr. L. 
de Borman and forms part of his master’s thesis at the 
University of Louvain, 1936, from which Figs. 5, 6 and 7 
are taken. 
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Fic. 5. The main cones. Energies in millist6érmers. Positive 
particles, northern hemisphere. 


previous paragraph there is also the region of 
penumbra. While the full analysis of the latter 
must be left to another occasion preliminary 
results obtained with the differential analyzer 
and earlier by one of us’ show that there are 
bands of light of complicated structure running 
alongside the boundary of the main cone alter- 
nating with bands of darkness and extending 
eastwards (for positive particles) clear over to 
Stérmer’s limiting cone’ of total darkness at the 
value of @ corresponding to the given latitude 
and energy and to y;=1. While the chief contri- 
bution to the azimuthal effect as described later 
in this paragraph undoubtedly comes from the 
main cones the theory here presented must be 
improved later, particularly as regards the east- 
west asymmetry, so as to take the penumbra 
into. account. Another correction must also be 
made to take care of the eccentricity of the 
earth’s magnetic center. This correction at a 
given point on the earth affects only the energy 
scale and can readily be computed by methods 
already given in a previous paper.’ 

A consideration of the main cones drawn in 
Figs. 5 to 7 shows immediately that the range of 


7G. Lemaitre, Ann. de la Soc. Sci. de Bruxelles A54, 162 
(1935). 

°C. Stérmer, University Observatory, Oslo, Publication 
No. 10 (1934). 

*M.S. Vallarta, Phys. Rev. 47, 647 (1935). 
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Fic. 6. The main cones. Energy in millistérmers. Positive 
particles, northern hemisphere. 


energy which can reach a point on the earth at 
a given zenith angle z is not the same at all 
azimuths. Consequently there must be a varia- 
tion of intensity at a fixed zenith angle and 
different azimuths (azimuthal effect). Special 
cases of the azimuthal effect are the north-south 
and east-west asymmetries. Suppose now that 
the circles z=const. are drawn on the diagrams 
representing the cones (Figs. 5 to 7). These 
intersect the boundaries of the cones of different 
energies at azimuths which can be read off 
directly on the figures. In this way the curves, 
Fig. 8, are obtained which give the minimum 
energy capable of reaching a given latitude at a 
zenith angle of 45° and azimuths shown by the 
curves. Similar curves for other zenith angles 
may be readily found. 

Unfortunately no systematic comparison with 
the experimental azimuthal effect is possible at 
the present time because no systematic measure- 
ments covering the complete turn of the horizon 
have been made. The most complete data 
available are those of Johnson'® who has made 
measurements at different azimuths and fixed 
zenith angles with a triple coincidence counter 
system covering more or less thoroughly the 
horizon from east to west through south. These 
measurements were made in Mexico (A= 29°), 


1° T. H. Johnson, Phys. Rev. 45, 569 (1934). 
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Panama (A=20°) and Peru (A=0°) at zenith 
angles 45° and 30° and heights of 6.8 m, 10 m 
and both 6.8 m and 6 m, respectively below the 
top of the equivalent water atmosphere. While 
no significant quantitative conclusions seem 
warranted at the present time there is good 
qualitative agreement between theory and ex- 
periment even without resorting to negative 
primaries. An important point brought out by 
the Peruvian experiment is that the addition of 
an atmospheric depth equivalent to 0.8 m of 
water does not change the relative order of the 
intensities with respect to one another. This 
must be expected if the azimuthal effect is 
essentially due to the action of the earth’s 
magnetic field on the primaries because, since 
the zenith angle is constant, the length of the 
atmospheric path traversed by the rays is the 
same at all azimuths. 

Besides the north-south asymmetry, already 
discussed extensively in our preceding paper,! 
an interesting case of the azimuthal effect is the 
east-west asymmetry first discussed by Rossi." 
Historically this was the first case of the azi- 
muthal effect to be discovered and has since 
been thoroughly studied experimentally ; perhaps 
we may be permitted here to voice the regret 





Fic. 7. The main cones. Energy in millistérmers. Positive 
particles, northern hemisphere. 


"' B. Rossi, Phys. Rev. 36, 606 (1930), see also Ricerca 
Scientifica 1, 561 (1934) for full references to his previous 
work. 
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Fic. 8. The azimuthal effect for z=45°. 


that it has been studied practically to the 
exclusion of all other azimuths. Fig. 9 shows a 
comparison of the energy ranges involved in the 
east-west and north-south asymmetries, where 
the former refers to positive particles and the 
latter to particles of either sign. It is immediately 
apparent that for purposes of gaining informa- 
tion as to the energy distribution of the primary 
radiation the study of the north-south asym- 
metry is superior to that of the east-west, first 
because the penumbra acts mainly on the east- 
west asymmetry but hardly on the north-south, 
second because the energy range is smaller and 
third because the north-south is independent of 
the particles’ sign. Except that the interpretation 
of the east-west asymmetry is made difficult 
because of the penumbra, a comparison between 
the east-west and north-south asymmetries fur- 
nishes a means, already pointed out in a previous 
paper, of determining the percentage of particles 
of either sign. A consideration of the curves in 
Fig. 9 shows that if the primary cosmic radiation 
is made up entirely of particles of one sign the 


+600 = -_—-- 






_MILLISTORMERS 


™ 


! 


_ ENERGY 


50 ; A P 90° 
ZENITH ANGLE. 


Fic. 9. The east-west and north-south asymmetries. Posi 
tive particles, northern hemisphere. 


north-south asymmetry must be considerably 
smaller than the east-west. A mixture of positive 
and negatives on the other hand, decreases the 
east-west asymmetry but leaves the north-south 
unaffected, provided the energy ranges covered 
by positives and negatives are not the same. 
Thus Johnson’s measurements in Mexico of the 
north-south and east-west asymmetries would 
seem to leave room for positive and negative 
primaries.” 

It should also be noted that once the cones 
have been calculated it is quite feasible to take 
the finite aperture of a multiple coincidence 
counter system into account. Failure to consider 
this point may easily lead into error when 
interpreting experimental results. 


5. THE MAIN CONE IN THE VICINITY OF 
THE ZENITH 
A particularly interesting region from the 
point of view of intensity measurements is the 
region around the zenith to which Zanstra'* has 
already called attention. We intend to devote 
this section to the study of this region. 


2 T. H. Johnson, Phys. Rev. 47, 91 (1935); 48, 290 (1935). 
13H. Zanstra, Naturwiss. 22, 171 (1934). 
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Fic. 10. Minimum energies in the region near the zenith. 


In the vicinity of the zenith the boundaries of 
the main cones run almost parallel to one another 
and intersect the meridian plane at a small 
angle y which increases with latitude. The value 
of the energy r for which the boundary of the 
main cone goes through the zenith gives the 
minimum energy that a particle must have to 
reach the earth in the vertical direction at a 
given latitude. This can be readily found by 
interpolation from the data contained in Figs. 5 
to 7 or from the equation of the line @=0 


(22) 


where the dependence of y; on longitude in the 
zenith direction is given by Fig. 5 of our pre- 
ceding paper.' This minimum energy is given in 
Fig. 10 together with other curves the signifi- 
cance of which will be taken up later. The angle 
¥ is plotted as a function of latitude in Fig. 11 
together with the gradient dr/dz at the zenith. 
It is seen that the gradient varies slowly with 
the latitude, nevertheless it is of importance 
when interpreting vertical intensity measure- 
ments made with triple coincidence counter 
systems of large aperture. The orientation of the 
counter system with respect to the direction y 
is also a significant datum of the experiment. 
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Fic. 11. The angle between the main cone and the meridian 
plane and the energy gradient at the zenith. 


The importance of the knowledge of the least 
vertical energy, emphasized by Zanstra, arises 
from the fact that, if the spectrum of the cosmic 
radiation is a continuous function of the energy, 
then intensity measurements in the vertical 
direction made beginning at the equator and 
extending sufficiently far north or south give a 
ready means of estimating with fair accuracy 
the energy lost by a cosmic-ray particle while 
traversing the atmosphere. It is sufficient to 
find the latitude beyond which the vertical in- 
tensity remains constant. Unfortunately neither 
the vertical intensity measurements of Auger 
Leprince-Ringuet'* nor those of Clay'® 
between A=10° and A=20° 
lead extend sufficiently far north or south to 
make this estimate possible. 

Another interesting question which can now 


and 


under 20 cm of 


be answered is as to the energy range which can 
reach a point of the earth of latitude \ coming 
from directions at a zenith angle z. For latitudes 
up to 30° the energy range involved can be read 
off directly from Figs. 5 to 7. The corresponding 
curves are plotted in Fig. 10. Their significance 
will be clear from the following example. Suppose 
we take \= 30° and ask what is the energy range 
which may arrive within the cone bounded by 
the directions at 45° with the zenith. The lower 
curve for z=45° then that 362 milli- 
stérmers is the minimum energy that a particle 


shows 


4 P. Auger and L. Leprince-Ringuet, Nature 133, 138 
(1934). 
1 J. Clay, Physica 2, 308 (1935). 
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must have to arrive at 45° from the zenith (and 
in a certain azimuth) while the upper curve 
gives that if it has a minimum energy of 540 
millist6érmers it may arrive from any direction 
within 45° of the zenith. Similarly 408 milli- 
stérmers is the least energy a particle must have 
to arrive at 30° from the zenith, while 510 
millistérmers is the least energy required to 
arrive from any direction within 30° of the 
zenith. In addition the limiting energies for 
which the main cone is completely open is 
given by the uppermost curve. This curve 
consists of two parts: first, the locus of the 
vertices of the periodic orbits (cf. Fig. 1) and 
then the limiting energy for the limit y, = 0.78856 
for which periodic orbits disappear. This curve 
takes into account only the trajectories of the 
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first kind.'® Finally the lowest curve gives the 
limiting energy for Stérmer’s limit y,=1 for 
which all directions are forbidden. 

It is a pleasure to renew the expression of our 
gratitude to the various persons, mentioned in 
our preceding paper, who have helped us to carry 
out the research of which this paper is another 
fruit. One of us is indebted to the Massachusetts 
Institute of Technology for continued support 
while this investigation was being completed at 
Louvain, and to the University of Louvain for 
providing numerous facilities during the period 
of his residence there. 

16 For trajectories of the second kind a somewhat higher 
energy limit can be readily found from the relation 
e*~=<2 cos \ above which the acceleration d?x/do? is neces 
sarily positive. For the same reason negative values of y; are 
of no interest. 
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Excitation-Curves for Fluorine and Lithium 


L. R. Harstap, N. P. HEYDENBURG AND M. A. Tuve, Department of Terrestrial Magnetism, Carnegie Institution o 


Washington 


(Received July 6, 1936) 


Observations on the alpha-particles emitted by lithium 
and the gamma-rays emitted by lithium and fluorine when 
bombarded with protons of energies up to 1000 kv have 
been improved by the use of a corona-free, 10,000-megohm 
voltmeter-resistor for voltage-measurements. By this means 
the accuracy of our measurements has been brought to 
two percent on an absolute scale and about one percent on a 
relative scale. Oscillograph-studies of voltage-fluctuations 
have shown that up to 1000 kv the voltage is constant to 
+1.4 percent ‘‘peak-ripple.” For the distribution-curve 
for voltage versus time the ‘‘half-maximum”’ width is about 
one percent at 1000 kv. Results thus far obtained for the 


INTRODUCTION 


UR work on proton-disintegrations last 

year' demonstrated the existence of sharp 
resonance-effects in nuclear disintegrations which 
called for voltage-control considerably more 
refined than that given by any apparatus then in 
existence. Such work depends essentially on the 
accurate reproducibility of specified voltages, 
and this requirement is not sufficiently well 


1 Hafstad and Tuve, Phys. Rev. 48, 306 (1935); also 
Tuve, Hafstad and Dahl, Phys. Rev. 48, 315 (1935). 


gamma-ray resonances are as follows: 


Voltage Element Half-width 
328 kv F < 4k 
440 kv Li 11 kv 
892 kv F <12 kv 
942 kv F <15 kv 


There is an indication of a weak multiplet structure in 
fluorine in the region between 500 and 700 kv with a broad 
but fairly prominent ‘‘resonance”’ at 650 to 700 kv. The 
existence of a resonance in lithium at 850 kv was not 
confirmed. 


satisfied even by rectifier and condenser installa- 
tions, unless means for voltage-measurements 
more accurate than spark-gap or particle-range 
determinations are provided. Unless very special 
precautions are taken, sphere-gap measurements 
are rarely reproducible to better than five per- 
cent. Even the official calibration-curves for 
sphere-gaps have been changed by more than 
ten percent during the past two years. The 
absence of reliable information on the range- 
energy relation for protons would prevent the 
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use of this method even though 
which would permit both rapid and accurate 


determinations of ranges were devised. 


apparatus 


In addition to the problem of specifying the 
mean voltage that of determining 
quantitatively the deviations from the mean 
during any given interval of time. This informa- 


used, is 


tion is essential in specifying such things as 
widths of resonance-levels, which have a theo- 
retical meaning and value only when it is 
demonstrated that the observed widths are not 
produced by inhomogeneities in the incident 
ion-beam. 


TECHNIQUE 


Tube 

Our apparatus is essentially as previously 
described. A new glass tube has been installed 
in which provision was made for permitting 
accurate axial alignment of the electrodes. In 
this tube the baffle-plates at each electrode have 
been eliminated so that the pumping speed 
through the tube is essentially that given by the 
12-inch diameter of the glass tube rather than 
that of the inside diameter of the electrodes. 
This arrangement will permit the use of much 
larger total ion-currents as soon as correspond- 
ingly fast pumps can be provided. At present 
the pumping speed at the base of the tube is 
approximately 30 liters per second, which is 
small compared to that used in most high- 
voltage installations. 

As was expected, the new tube, with its 
properly aligned electrodes, practically elimi- 
nated the annoying fluctuations in position of 
the spot which constituted the principal difficulty 
encountered in the work reported last year. 
The effectiveness of the focusing of such a tube 
is attested by the fact that with the same tube 
a concentrated ion-beam can be obtained at any 
voltage from 200 to 1200 kv. 


Corona-free, 10,000-megohm voltmeter-resistor 


For the highest possible accuracy in voltage- 
measurements, a wire-wound voltmeter-resist- 
ance would be desirable. Such units are available 
commercially for x-ray voltages, but the cost of 
such units arranged for 1000 kv would be 
prohibitive, and a corona-free design could hardly 


be attained by reason of the dimensions of the 
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voltmeter-resistor. 


corona-free 


Fic. 1. Construction of 


unit resistors, especially if a full-scale current of 
100 
per volt). However, a close approximation can 


microamperes were required (10,000 ohms 


be achieved by using modern metallized ‘‘grid- 
leak”’ resistors which give voltage-current curves 
that are nearly ohmic in character. Any deviation 
from linearity can be corrected for by a calibra- 
tion of such resistors against “purely ohmic” 
wire-wound units. 

In the assembly of such resistors in series to 
be used with very high voltages, the greatest 
technical problem is the elimination of corona. 
Of several possible solutions the simplest ap- 
peared to be that indicated by Fig. 1, which 
has now been proved to be entirely satisfactory. 
The individual resistors of 10 megohms each, 
rated at one watt or 1000 volts maximum, are 
connected in series in groups of 20, slipped into 
tube to 


a_ thick-walled rubber insure 


spacing, and then connected between brass disks 


proper 


12 inches in diameter, which are protected 
against edge-corona by j-inch copper tubing 
soldered to the circumference of each disk. The 
disks are separated by 2-inch Isolantite rods as 
spacers. Fifty such units make up the 10,000- 
megohm resistor. This construction limits the 
2-inch 


voltage-gradient to 20 kv section 


which is safely below corona-values, provided 


per 


no dust or lint is permitted to collect on the 
surfaces. To guard against the latter the entire 
assembly the ‘“Textolite’’ 
cylinder shown in Fig. 2. This cylinder was on 


was mounted in 


hand from previous work. As may be seen from 
Fig. 2, it is somewhat short for its purpose and 
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Fic. 2. 1200-kilovolt electrostatic generator showing tube 
and voltmeter-resistor. 


at present sparks to earth over the outside of this 
cylinder determine the maximum voltage, 1150 
kv, attained by this generator. Some difficulty 
was experienced when this resistance-unit was 
first mounted due to corona at the junction 
between the resistor and the sphere. This was 
largely eliminated when care was taken to design 
in such a way that the equipotential through the 
first corona-disk agreed roughly with the corre- 
sponding equipotential for the sphere itself. 

The resistance-unit shown in Fig. 2 measures 
only the voltage between the outside sphere and 
ground. With our somewhat inconvenient ar- 
rangement of concentric spheres the voltage 
between the one-meter and the two-meter 
spheres must be measured separately. A similar 
resistance-unit of six 20-kv sections was accord- 
ingly installed between the two spheres, the 
current being indicated by a meter installed at 
the surface of the outside sphere and read by a 
telescope. 

Voltmeter-calibration 


Through the courtesy of Dr. L. S. Taylor of 
the National Bureau of Standards we were able 


to calibrate our resistance against his wire-wound 
standards in 5-section, 100-kv units. Results on 
separate determinations agreed to within two 
percent. The voltage corresponding to any given 
current through the entire resistor is then 
obtained by adding the voltages of the separate, 
independently calibrated units, and the total 
voltage is obtained by adding the separately 
determined inside voltage (between spheres) 
and outside voltage (to ground). 

The possibility that part of the measured 
current at high voltages may be due to corona 
between sections of the voltmeter has been 
eliminated by three independent tests. First, 
during calibration, the units were carried to 40 
percent over-voltage before corona occurred, 
thus allowing a safe margin above usual opera- 
tion-conditions. Second, after installation the 
voltage-scale was tested for linearity by setting 
the mass-1 spot to a given deflection at the 
highest voltage and then reducing the voltage 
until the mass-2 spot was deflected to the same 
position. By this test the voltmeter read two 
percent low at the highest voltages attained. 
This error is therefore not only small but 
opposite in sign to that expected from corona. 
Finally, making use of an effect first reported 
to us by Herb but since frequently tested here, 
an effort was made to detect the change in 
current-reading which should be produced when 
any incipient corona was stopped by the intro- 
duction of CCl, vapor into the voltmeter- 
resistor casing. No effect was observed even at 
the highest voltages. 


Charging current 


In all work done so far our maximum target- 
currents have been limited by our pumping 
speed so that maximum charging currents have 
not been necessary and accordingly we have 
operated our belts at reduced speeds to reduce 
mechanical vibration. The added current-drain 
due to the voltmeter could be readily compen- 
sated by increased belt-speed, if desired. How- 
ever, during some tests on a disk-type electro- 
static generator? it was noted that charging 
currents could be greatly increased by intro- 
ducing a thin dielectric between the current- 


20. Dahl, Rev. Sci. Inst. 7, 254 (1936); E. H. Bramhall, 
Rev. Sci. Inst. 5, 18 (1934). 
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carrier and the earthed inductor, and making a 
simultaneous increase in spray-voltage. Adoption 
of this technique on our belt-machine doubled 
our previous charging currents so that we now 
operate at reduced speed and with only the 
“ingoing”’ belts thus the 
unsteadiness arising from incipient sparking at 


charged, reducing 
the belt-openings when operating at the higher 
voltages. 


Voltage-fluctuations 


As before, variation of voltage was accom- 
plished by moving a set of corona-points toward 
or away from the sphere. At all except the 
highest voltages, where sparking occurred, it 
was immediately obvious from the steadiness of 
the voltmeter-needle that no erratic voltage- 
fluctuations were occurring. It is important to 
note that this steadiness is a property which is 
apparently inherent in a corona-limited, con- 
stant-current machine, and is attained without 
recourse to smoothing condensers or compen- 
sating devices of any kind. 

Since the period of a portable microammeter 
is usually of the order of one second, there was 
of course a possibility of relatively fast fluctua- 
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Fic. 3. Voltage-fluctuations, 2-meter electrostatic generator. 
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tions which were smoothed out by the period of 
the meter. To eliminate this possibility observa- 
tions were made with a cathode-ray oscillograph 
connected across a one-megohm resistance in 
series with the voltmeter-multiplier resistance. 
When due for the 
impedance of the resistor-unit to high frequency 


allowance was made low 


pulses (capacitative impedance equals resistive 
impedance at 100 cycles), it was immediately 
from the mean of 


that no deviations 


greater than one or two percent were occurring. 


clear 


To obtain a quantitative measure of the fluctua- 
tions, however, a photographic was 
necessary. For this purpose the 10,000-megohm 


record 


resistor was connected to a single-stage, vacuum- 
tube amplifier operating a high-sensitivity Dudell 
oscillograph (700 cycles) which recorded on a 
moving film. The results for the voltage on the 
outside sphere are shown in Fig. 3. From these 
records the distribution-curves for 250 equally 
spaced ordinates measured from each of these 
curves were obtained. Similar records were taken 
for the voltage-fluctuations on the inside sphere. 
The voltage between the spheres showed a 
noticeable 60-cycle ripple which was absent on 
the outside sphere; from this and other evidence 
the voltage-fluctuations may be considered inde- 
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Fic. 4. Analysis of voltage-fluctuations on 2-meter electro- 
static generator. 
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pendent. In this case the widths of the two 
distribution-curves for the fluctuations should 
be added in the same way as the probable errors 
for independent error-curves, namely, W142 
=(W,?+ W,*)!. The data for the separate distri- 
bution-curves as well as for the total voltage are 
given in Fig. 4. This study shows that even at 
990 kv the maximum “‘spread”’ of the fluctua- 
tions is only 28 kv or +1.4 percent “‘peak-ripple.”’ 


YIELD-CURVES 


Accurate yield-curves for either particles or 
gamma rays can give much information of value 
in the development of a theory of the nucleus, as 
has been demonstrated especially by the work 
of Breit.* While yield-curves are desirable for 
all possible elements, we have thus far concen- 
trated on lithium, on which much work has 
already been done,‘ and on fluorine® in order to 
establish an absolute voltage-scale by means of 
its many resonances. Observations will be ex- 
tended to other elements as time permits. 


’ Breit and Yost, Phys. Rev. 48, 203 (1935). 

* Cockcroft and Walton, Proc. Roy. Soc. A137, 229 
(1932); Henderson, Phys. Rev. 43, 98 (1933); and Oliphant 
and Rutherford, Proc. Roy. Soc. A141, 259 (1933). 

> McMillan, Phys, Rev, 46, 868 (1934). 


ALPHA-PARTICLES FROM L1‘+H! 


This problem has long been studied under 
gradually improving experimental conditions. 
The recent precision work® on this problem and 
consequent calculations of theoretical expecta- 
tions’? emphasized the need for more observations 
at the higher voltages. 

Our most recent results for a thick target are 
given in Fig. 5, with results of other observers 
included for comparison. In order to obtain 
reproducible results it was necessary to adopt a 
target-technique in which a freshly prepared 
target of lithium metal, evaporated onto a 
copper plate in vacuum, could be moved across 
the ion-beam so as to present a fresh surface as 
often as desired. This technique has not yet 
been extended to thin targets because of the 
obvious difficulty of insuring uniformity. Un- 
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Fic. 6. Comparison theoretical and observed excitation- 


functions for lithium. 


6 Herb, Parkinson, and Kerst, Phys. Rev. 48, 118 (1935); 
Heydenburg, Zahn, and King, Phys. Rev. 49, 100 (1936). 
7 Ostrofsky, Breit, and Johnson, Phys. Rev. 49, 22 (1936). 














EXCITATION-CUR 


fortunately thick target data can be compared 
with theory only by assuming a law for the range- 
energy relation for the protons in lithium. This 
has been done following Herb*® and the results 
compared with theoretical values’? in Fig. 6. 
It may be seen that our recent results deduced 
from thick-target data agree well with prelimi- 
nary values obtained from thin-target observa- 
tions last year. Plotted in this way the data 
exhibit irregularities not present in the original 
data and lead one to suspect that the range- 
energy relation for protons in lithium may be 
something other the V! 
However, it is clear from these data that 


assumed. 
the 
cross section is not constant above 400 kv as 


than law 


reported by Henderson. 

Still more detailed lithium 
will be needed before the thin-target yield-curve 
“excitation-function”’ 


observations on 


or, more explicitly, the true 
for this process can be definitely established. 


GAMMA-RAyYs FROM L1’?7+H! 


Our work on the gamma-rays from lithium 
reported last year has been repeated (1) in order 
the 
voltage-scale and (2) in order to determine, if 


to establish resonances on an absolute 
possible, the widths of the resonance-levels. 
Observations were made, as before, using the 
Lauritsen type of electroscope. Representative 
results, the last of a series of curves, are given in 
Fig. 7. By comparison with our previous work 
it will be noted that the lower resonance occurs 
at practically the same position on this voltage- 
scale (440 kv) as on the rather arbitrary scale 
used last year (450 kv), again disagreeing 
markedly with the resonance-position of 650 kv 
reported by the Pasadena observers.* In Fig. 7 
it may also be noted that there is no confirming 
evidence for the resonance at 850 kv reported 
in our paper last year. It is probable that the 
gamma-rays observed were due to a contamina- 
tion by fluorine which has an extremely strong 
resonance-line at approximately this voltage. 
This possibility is favored by the fact that in 
order to obtain a sufficiently thin lithium target 
(less than 20 kv stopping power) the amount of 
Li present was reduced dangerously close to the 
contamination-levels of any materials giving 


® Crane, Delsasso, Fowler and Lauritsen, Phys. Rev. 48, 
125 (1935). 
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Fic. 8. Determination of half-width of 440-kv lithium 


gamma-ray resonance. 


large yields. In the present instance fluorine is 
especially to be suspected since CaF, had been 
used as a target in experiments just preceding 
the observations on lithium. 

The importance of obtaining estimates of 
absolute yields for curves such as the above is 
recognized, but the problem is too difficult to 
permit dependable values to be given at the 
present time. Comparisons with radium are 
misleading, for the gamma-rays involved in the 
present reaction are of extremely high energy as 
has been shown by the work of Lauritsen and 
his collaborators.* It that 
both energies and intensities of gamma-rays will 
have to be determined by the laborious cloud- 


appears at present 


chamber methods. 

In addition to the location of a level, the 
question of the width of the resonance-level is of 
theoretical importance. A special effort was made 
to set limits on the 440-kv resonance in lithium. 
In Fig. 8 are shown integrated curves for reso- 


‘“‘half-widths’’® with the ob- 


nances of various 


® It is customary to use the width at half-maximum or 
“half-width” as a measure of the sharpness of a resonance- 
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served points for the lithium-curve superimposed. 
From this comparison one can conclude that the 
“half-width” of the observed lithium-resonance 
curve is about 12 kv. However, part of this 
width must be ascribed to voltage-fluctuations, 
which produce a spread in the energies of the 
ions in the incident beam. Referring to Fig. 4 
one estimates this width to be about 4 kv, leaving 
a width of about 11 kv=(12?— 4°)! for the true 
resonance-width for this 440-kv process in 
lithium. This conclusion is open to the possible 
objection that the measurements of the voltage- 
fluctuations were not made at the same time as 
the gamma-ray measurements. However, the 
consistent behavior of the generator as indicated 
by the resistance-voltmeter and in our oscillo- 
graphic tests seems to preclude the possibility 
that the voltage-fluctuation was abnormal by a 
factor of three during the taking of the observa- 
tions on lithium. Further evidence that the above 
width is not due to the apparatus is given by 
the fact that a much smaller width was observed 
in fluorine as will be discussed below. 


GAMMA-RAyYs FROM F+H! 


Our results for this process are given in Fig. 9. 
The observations were primarily made to estab- 
lish the several fluorine resonances on an absolute 
voltage-scale and to test the “resolving power”’ 
of our apparatus in the higher voltage-regions. 
The points shown in Fig. 9 are numbered in the 
order in which they were taken, each point 
requiring from three to ten minutes for a com- 


curve; see, for example, Slater and Frank, Introduction to 
Theoretical Physics, p. 37. 


plete observation. The agreement between points 
1, 2, 3 and points 36, 37, 38 taken several hours 
later is particularly satisfying. The existence of 
the close doublet with components at 894 and 
942 kv has been checked by several inde- 
dently observed curves and appears to be real. 
The fact that these two lines are close probably 
accounts for the broad resonance-peak observed 
in this region in the observations made last year, 
when with less precise voltage-measurements our 
apparatus effectively had less ‘‘resolving power.” 

In the region from 400 to 700 kv the rapid 
rise in the curve suggests the possibility of 
several minor resonances. An effort was made to 
test this possibility, but no conclusive results 
were obtained due to difficulties in devising a 
suitable target. Crystals of CaF, were shattered 
by the localized heat of bombardment, and 
CaF, powder was blown about by electrostatic 
forces; hence, to obtain more precise data, 
experimentation with other compounds of fluo- 
rine will probably be required. 

Estimates of the “half-widths”’ of the 328-kv, 
892-kv, and 942-kv resonances made in the same 
way as for lithium give ‘‘observed’’ values of 
about 4 kv, 10 kv, and 15 kv, respectively. 
Since from the table in Fig. 4 these values 
appear to be practically equal to the “‘half- 
widths” of the voltage-distribution curves, the 
only permissible conclusion is that the “true” 
widths are probably smaller than these ‘‘ob- 
served”’ values. 

DIsCUSSION 
Process of gamma-ray emission” 


Any explanation of the observed resonance for 
the emission of gamma-rays when Li’ is bom- 
barded with protons has to make use of a 
temporary semistable state (virtual level) of 
Be’. The mean life of this level is determined by 
the half-value breadth of the level and is thus 
(1/4m)(510/11)(A/mc?) =3.7h/me=3X10- sec. 
The width of the level must be due to its damping 
by all possible modes of decay." 


10 We are indebted to Dr. Placzek and Professor Teller 
for valuable discussions of the problem of gamma-ray 
emission. We understand that these ideas will be published 
in a paper in the Physical Review which will appear shortly. 
The theoretical treatment given here in connection with the 
Li gamma-ray problem is due to Professor Breit. 

- ' for example, Breit and Wigner, Phys. Rev. 49, 519 
(1936). 
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EXCITATION 


The virtual level of Be® may either radiate 
directly into a lower level of Be® or it may decay 
by disintegration into two alpha-particles at 
least one of which will then have to be excited. 
It is hard to account for a number of gamma-ray 
wave-lengths using only quantum jumps between 
different levels of Be* because such levels are 
more likely to disintegrate into two normal 
alpha-particles than to radiate gamma-rays. For 
this reason Crane and Lauritsen have assumed 
that the Be® nucleus first of all disintegrates 
into a normal and an excited alpha-particle; 
the 
to radiate the complex gamma-ray spectrum. 


later excited alpha-particle is supposed 
According to a recent report® the spectrum 
consists of a single line having an energy of 
17 MEV and the main support for this picture 
of the process appears to be withdrawn. 

Some preliminary observations made in Wash- 
ington indicate that there may be also present a 
gamma-ray of 8-MEV energy. These observations 
were made on the magnetic curvature of Compton 
electrons expelled from a thin sheet of glass 
(about 500 kv thick) in a Wilson cloud-chamber 
by the Li+H! gamma-radiation. It was found 
MEV 
which definitely originated in the thin glass sheet 
22 fell in the region between 7.5 to 9.5 MEV. 
The only higher-energy Compton electrons were 
in the vicinity of 11 MEV, where three tracks 
were found. 

According to the above it is still somewhat 
the 
consists of a single line or not. If it does, the 
the 
virtual state radiates directly into the normal 
level of Be*®. Afterwards the Be*® nucleus may or 


that of 31 tracks with energies above 2 


uncertain whether gamma-ray spectrum 


simplest explanation is to suppose that 


may not dissociate depending on its still some- 
what uncertain mass. We failed to observe any 
irregularity in the excitation function for the 
ordinary 8-cm alpha-particles in the vicinity of 
440 kv either with thick or with thin targets. 
In the thick-target observations alpha-particles 
having a range greater than 7.6 cm would have 
been counted. The thin-target experiments were 
concerned with particles having greater pene- 
trating power than 1.2-cm alpha-particles. This 
indicates that the resonance to gamma-rays has 


2 Delsasso, Fowler, and Lauritsen, Seattle Meeting, 
Am. Phys. Soc. (1936). 
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to do with states of excitation of the Be*® nucleus 


which do not dissociate into normal alpha- 
particles. Quantitative measurements of gamma- 
ray intensities are at present still too uncertain 
to be sure of their number. Since, however, the 
probability of mechanical disruption is, under 
ordinary circumstances, many times greater than 
the probability of radiation it is probable that 
the virtual level of Be® is barred from disintegra- 
tion into alpha-particles through the operation 
of some rather rigid selection-rule. This conclu- 
sion appears to be very probable for any mechan- 
ism involving radiative transitions within the 
Be® nucleus. It does not follow on the hypothesis 
of disintegration into a normal and an excited 
alpha-particle because such a disintegration 
could conceivably be so much more probable 
than disintegration into normal alpha-particles 
that no change in the yield of the 8-cm alpha- 
particles would be expected. 

According to our present admittedly crude 
measurements of the yield of gamma-rays the 
collision cross section at resonance is roughly 
10-*? cm?®. Substituting this into formula (14) 
given by Breit and Wigner and using /I'=5 ky 
one finds that either AT,~1 volt and hT',~5 ky 
or else AT',~1 volt and Al',~5 kv. This means 
that 


accounting for 


there are essentially two possibilities for 
the order of magnitude of the 
(1) The damping of the virtual level is 


entirely to the possibility of 


yield : 
due practically 
break-up by dissociation into Li7+H!', the order 
of magnitude of hI, is then so small that it is 
most reasonable to explain the process as being 
due to the emission of gamma-rays to a lower 
level of Be®. (2) The damping due to disintegra- 
tion into Li’+H! could be supposed to be small 
(kAT,~1 volt) and all of the width of the level 
would be attributed so some other mode of 
disintegration. The high value of hI, would not 
fit in with the expected probability of gamma- 
radiation within the nucleus. This mode would 
then involve some such process as the disinte- 
gration into excited and normal alpha-particles. 
The 
somewhat improbable because it indicates that 


very small value of AT,~1 is however 


transitions from Li7+H! to the virtual level are 


much less likely than ordinary calculations 


indicate. [It 
interpret hI’, as due to protons and /T, as due 


would be possible similarly to 
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to alpha-particle disintegrations using a large I’, 
and a small IT’, arising from some selection rule. | 
A process with a small I, and a large I, cannot, 
however, be completely excluded, particularly 
since the gamma-ray spectrum is not yet certain. 
It would mean that the virtual level happens 
to be practically decoupled from Li?7+H! and 
at the same time it would have to dissociate 
more readily into a normal and an excited alpha- 
particle than into two normal alpha-particles. 
Thus even though a small I, and a large I’, are 
not logically excluded the number of conditions 
to be satisfied is greater and the radiative 
capture of Be* hypothesis fits in more naturally 
with present evidence. Experiments on the 
anomalous scattering of protons by lithium at 
the resonance voltage should differentiate be- 
tween the two alternatives. On the hypothesis 
of the radiative transitions within Be® a scatter- 
ing cross section of the order of 5X10-*4 cm? 
would be expected while alpha-particle disinte- 
gration would be associated with a much smaller 
scattering cross section. 

According to quantum mechanics energy-levels 
can be classified into even and odd and for 
closed systems there are no transitions between 
levels of these two classes; on the other hand 
electric dipole radiation jumps occur only be- 
tween levels of one class and those of another™ 
(Laporte’s rule). It has been mentioned that if 
radiation takes place between levels of Be® a 
selection rule is needed to prevent the disinte- 
gration of the resonance level into two normal 
alpha-particles. Since two normal alpha-particles 
obey Einstein-Bose statistics and have no spin, 
the resonance level could be assumed to be odd 
and disintegrations into normal alpha-particles 
would be then prohibited. At the same time it 
would be necessary to arrange for stability 
against decay into a normal and an excited 
alpha-particle. Such stability can be automati- 
cally provided by supposing that the first odd 
excited level of the alpha-particle lies more than 
17 MEV above the normal level. This is not in 
contradiction with Feenberg’s calculations on 
the position of the lowest 'P level in the alpha- 
particle and the °P level is probably above the 
'P level on account of the presence of exchange 
forces in the nucleus. On the above picture an 


18 E. Wigner, Gruppentheorie. 
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electric dipole radiation can take place from the 
virtual to the normal level of Be’. If the gamma- 
ray is actually quite monochromatic then one 
has to suppose that there is no other level 
between the resonance level and the normal state 
to which a transition with radiation can take 
place. It is very difficult, however, to be sure 
that even the 17-MEV component of the gamma 
ray is monochromatic and so far as we are able 
to tell at present there is just as much evidence 
for a gamma-ray multiplet having components 
between 13 MEV and 17 MEV as for an entirely 
monochromatic gamma-ray. Confining ourselves 
for the present to the possibility of the first 
radiative transition occuring in Be’ the following 
type of energy-level diagram of Fig. 10 can be 
considered as possible. The resonance level V is 
is responsible for the capture of the protons. 
It is supposed to be ‘‘odd”’ so as not to disinte- 
grate into normal alpha-particles. No irregularity 
in the excitation-curve for the 8-cm particles 
will be then expected because odd and even 
levels do not interact. The normal level of Be* 
is designated by N and is supposed to be even. 
The 17-MEV ray may arise through electric- 
dipole radiation in the jump VN. In addition 
there may be either odd or even levels N’, N” 
within say 3 MEV of N. If these levels are 
even, there may be electric dipole-transitions 
VN’, VN” accompanied by subsequent disinte- 
gration into short range alpha-particles. If the 
levels N’, N” are odd no disintegration into 
alpha-particles would be likely to occur from 
them and instead soft gamma-rays would be 
emitted, in such jumps as N’N. It is more 
probable that the levels N’, N”’ are even and a 
search for short range alpha-particles appears 
to be worthwhile. In addition to the group of 
levels N’, N” there appears to be no reason 
against admitting an odd level J approximately 
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between V and N. Magnetic dipole or quadripole 
radiation would then account for 8- MEV gamma- 
rays, an indication of which has been obtained. 
Since, it was supposed that the first odd level of 
the alpha-particle is more than 17 MEV above 
normal in order to account for the stability of V, 
the level J would also be stable towards me- 
chanical disruption and the absence of 4-MEV 
alpha-particles would follow. It is as yet not 
that such 
particles are not present in small amounts and 


completely certain experimentally 


the ‘‘odd”’ character of the level J is, therefore, 
not an absolute necessity. It should be noted 
that in thick targets with protons well above the 
resonance voltage the 4-MEV alpha-particles 
would be produced somewhat inside the target 
and would consequently be harder to detect. 

In the above discussion and diagram levels are 
designated as odd and even in accordance with 
what appears to be likely from present con- 
jectures about the normal state of Be®. Clearly 
the 
concerned with using selection rules. Thus the 


the only essential part of argument is 
level V could be supposed to be even with odd 
spin and would still be stable towards disinte- 
gration into normal alpha-particles. It appears 
to be premature to list all possible arrangements 
of levels which might be consistent with present 
experimental facts. A certain degree of proba- 
bility can be attached however to the particular 
assignment of levels given above making use of 
some recent and as yet unpublished calculations 
of Wigner and Feenberg. According to these one 
expects an even 'D level to lie about 3 MEV 
above N which should be a 'S level. The next 
even level ('G) would lie about 8 MEV above N 
and would be thus approximately in the position 
of J. In order to have an electric dipole transition 
to N one may suppose V to be an odd P level 
which would then be also likely to give radiation 
to N’='D. On this view there should be short- 
range (1.5-MEV) alpha-particles emitted from 
N’. Transitions to the 'G level would not occur 
except through the possible influence of higher 
multiple radiations. If transitions to the 'G level 
should occur some 4-MEV alpha-particles will 
be emitted because its mean life is short when 
estimated by means of formulas for barrier- 
penetration. 

If the virtual level V disintegrates by gamma- 
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FOR F AND Li 
ray emission all of its width must be attributed 
to the probability of dissociation into Li’+H!'. 
A theoretical estimate of the width can be made 
easily by regarding the proton as subjected to a 
central field due to the lithium nucleus. Using 
the formula" 


AT/T =! | G dp; G(ro) = G(ro) = Dp "0, 


where ro = nuclear radius, L = angular momentum 
of proton with respect to Li’, the graphs for 
©, ©,, given by Yost, Wheeler, Breit,'® and the 
WKB approximation for L=2, we obtain for 
the half-width 2A7T=400 kv for L=0, 20 kv for 
L=1 and } kv for L=2. In this calculation the 
integral in the denominator was approximated 
by G*p. The above estimates are likely to give 
the the 
into Li’+H! depends not 


too high values for width because 
disintegration of V 
only on the probability of the proton escaping 
but also on the chance of an internal nuclear 
rearrangement. There is thus no difficulty in 
accounting for the width using head-on collisions 
(L =0) and even those for L = 1, but it is probable 
that L=2 plays little part in the process. The 
normal configuration of Li’ is very probably odd. 
Since V must also be odd it is probable that the 
capture and dissociation occurs through the 
partial wave L=0. If, contrary to present views, 
the normal configuration of Li’? should be even 
then L=1 will have to be held responsible. 
This is less likely also since little leeway is left 
for the probability of internal rearrangement. 

Observations recently taken in this laboratory 
jointly with Dr. Rumbaugh of the Bartol 
Foundation indicate that there is stationary 
level in Li’ about 400 kv above the normal state. 
Presumably this level corresponds to an angular 
momentum of } and constitutes together with 
the normal state of Li’ a doublet p term. It is 
conceivable that the small cross section for the 
production of gamma-rays is partly due to the 
disintegration of the virtual level V into a slow 
proton with an excitation of Li’ to ‘‘p,”’ level. 
The probability of this happening at 440 kv 
incident proton-energy should be small, because 
the barrier-penetration factor for the 440 —400 

14 Breit and Yost, Phys. Rev. 48, 203 (1935). 

4 Yost, Wheeler and Breit, Phys. Rev. 49, 174 (1936) (the 
notation of this reference is used in the text above). 
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=40-kv protons decreases the probability of 
their interaction with the virtual level. It is 
possible, however, that at higher energies of 
incident protons evidence of inelastic collisions 
due to the excitation of the 400-kv level may 
be found. 

Considerations such as the above cannot at 
present be effectively applied to the fluorine 
resonances, which have been shown to be 
significantly sharp, since even the reaction or 
reactions involved are at present uncertain.®: '® 


VOLTAGE-SCALE FOR INTER-LABORATORY 
COMPARISON 


It has long been customary to measure 
voltages of ion-accelerators by giving the range 
in air of the ion-beam produced. However, the 
accurate specification of such ranges is by no 
means a simple matter, for such factors as 
straggling, ion-distribution, and current-density 
must be taken into account. The most readily 
determined specification of an ion-beam is 
probably a quantity which we might call the 
‘visual range,” this being merely the range of 
the particles in air as observed by the fluorescence 
produced when the beam is passed through a 
thin window into a darkened room. This quantity 
will be found to vary with current-density, room 
light, and different observers, and bears no 
definable relation to the ‘‘mean”’ range which is 
usually the quantity desired. In spite of these 
difficulties, however, it is repeatable to about 
+1 mm, and hence gives an accuracy of about 
+10 percent in voltage-measurements in the 
range in which this indicator can be conveniently 
used. For such purposes the relation R=2.5V}, 
where V is in megavolts, describes our results for 
the region between 500 and 1000 kv.'? The 
‘visual’ range will in general be longer than the 
“mean” range due to thin spots in windows, 
voltage-fluctuations (since the peaks will be 
measured), variations in current-density, and 
other such causes. 

For nonhomogeneous beams the visual range 
will of course give only the upper limit of 
energies. The effective voltage for disintegrations 


‘6 Feather, ‘Reports on progress in physics,” Proc. Phys. 
Soc., Vol. II. 

17 The voltage-scale used in this laboratory last year was 
based on a visual range of 2.85 cm at 1000 kv, 1.00 cm at 
500 kv. 


may then be any lower value, depending on the 
form of the ion-energy distribution-curve. Since 
in such cases sharp resonance-curves such as 
those described in this paper will no doubt be 
useful in calibrating ion-beams for both energy- 
and ion-distribution, we have endeavored to 
make our determinations as accurately as was 
immediately practicable. The original calibration 
against wire-wound standards insures that the 
measured voltages agree with the international 
scale to within two percent. This includes such 
errors as temperature-variation, which in the 
present design is relatively unimportant for the 
resistors are operated at only one-tenth of their 
rated capacity. While the location of a resonance 
is easily determined, the accurate specification 
of “half-widths,”” on which conclusions as to 
ion-distribution must be based, is difficult for 
two reasons. First, for sharp resonances, observa- 
tions must be made at intervals of about 5 kv, 
which at voltages over 500 kv requires that the 
mean energy be measured to less than one 
percent ; second, any asymmetry in the observed 
resonance-curves causes an ambiguity in the 
specification of ‘‘half-widths’” when plotted as 
in Fig. 8. This problem becomes even more 
difficult when the observed width approaches 
the “‘fluctuation-width” for in this case any 
slight asymmetry in the ion-energy distribution 
would also have to be taken into account. With 
these qualifications in mind the points at present 
available for the calibration of ion-beams of 
unknown energy or distribution are as follows: 


Voltage Element Half-width 
328 kv F < 4 kv 
440 kv Li 11 kv 
892 kv F <12 kv 
942 kv F <15 kv 


ACKNOWLEDGMENT 

The authors take pleasure in expressing their 
indebtedness to Professor G. Breit for his con- 
tinued interest and advice, to O. Dahl for 
invaluable aid in the design and construction of 
apparatus, and to Dr. R. G. Herb for his 
enthusiastic assistance in the construction and 
testing of the voltmeter-resistor while a guest at 
our Laboratory. Our thanks are especially due 
to our Director, Dr. J. A. Fleming, for his 
constant encouragement and unfailing support. 








te 
m 


in 


an 


chi 
me 
tic 








SEPTEMBER 15, 1936 PHYSICAL 


The Precise Measurement of Three Radium B Beta-Particle Energies 


F. T. RoGers, Jr., 


(Received July 1, 1936) 


The Hp of the three most intense Ka B B-particle lines 
have been measured by a precise magnetic spectrograph 
method and found to be 1406.0, 1671.1, and 1931.5 gauss 
cm. An explicit analysis has been made of the Cotton 
balance equations for J/ in this experiment and on the 
basis of this a complete analysis was made of the maximum 
relative error to be expected in the accepted values of Ip. 
It is shown that the maximum relative error to be expected 
in the least Hp is 1 part in 3000 and that within this 
possible error the greatest H/p agrees with the value obtained 


INTRODUCTION 


N 1924 Ellis and Skinner! measured I/p (// is 
the magnetic field necessary to produce a 
radius p of curvature in the path of a 8-particle) 
for many of the 8-particle energies of Ra B; for 
the three most intense lines of Ra B they got, 
correctly to within what they thought was 1 part 

in 500, 

IIp = 1410, 1677, 1938 gauss cm. 


In 1934 Ellis? repeated this work with more pre- 
cision and got 


Tp = 1400.4, 1665.9, 1925.5 gauss cm 


(in which the last figures are not entirely signifi- 
cant) for these same lines. In 1934 also, Scott*® 
performed a precise experiment on the most in- 
tense of these lines using the Cotton balance 
method of measuring the magnetic field rather 
than an inductive method; he got for this most 
intense line 
ITp=1931.8 gauss cm 


and hoped it to be correct to 1 part in 10,000. 

The present experiment was undertaken to 
check Scott's work and method and to extend the 
measurements to other lines in the Ra B §-par- 
ticle spectrum. The values of 


7p = 1406.0, 1671.1, 1931.5 gauss cm 


were obtained, the last of which is in good agree- 
ment with Scott's value for the most intense line. 
The Ho's for the other two lines as found by this 


1 Ellis and Skinner, Proc. Rov. Soc. Al05, 165 (1924). 
2 Ellis, Proc. Roy. Soc. Al43, 352 (1934). 
Scott,’ Phys. Rev. 46, 633 (1934). 
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by Scott. It is suggested that the accuracy of | part in 
10,000 hoped for by Scott is rather too optimistic. The 


energies of these three 6-particles are 1.512, 2.044, and 
2.610 in units of 10° electron volts. It is shown that the 
least of these cannot be expected to be any more accurate 
than 1 part in 1000. The conclusion is that except for 


systematic errors which might have been present and 


undetected throughout this experiment and that of Scott, 
these experiments have established the values of these 
Hp's to within 1 part in 3000. 


experiment yield slightly different relative values 
of Hp from those previously obtained. 

A more thorough investigation of the theory 
of the Cotton balance method of measuring // 
was carried out for this experiment. On the basis 
of it an analysis was made of the errors which 
might be present in the finally accepted values 
of /Ip. The conclusion was reached that in this 
experiment an accuracy of 1 part in 3000 
attainable if the value of // used in the finally 
accepted //p's is determined from the mean of at 


1S 


least forty separate, equally weighted determi- 
nations of the mass used with the Cotton balance 
and if the values of p used finally are the means of 
seven separate, equally weighted determinations. 

The good agreement (within 1 part in 6400 in 
Hp) of the present experiment with that of 
Scott would seem to indicate that this type of 
measurement has at last established the absolute 
values of the J/p’s for the three most intense lines 
in the Ra B 8-particle spectrum to within 1 part 
in 3000 except for possible systematic errors 
which may have been present and undetected in 
both experiments. 

The energies corresponding to the three J/p's 
we have determined are found to be 
V=1.512X10°, 2.044 10°, 2.610 10° electron 


volts 


and can be considered as accurate to not better 
than 1 part in 1000. 


APPARATUS 


The apparatus used in this experiment was, 
except for a few small changes, that used by 
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Scott in his work and described by him in his 
paper.’ It consists of a large permanent magnet 
with cobalt steel poles and with a 1.52 cm air 
gap faced on each side by soft iron disks 15 cm 
in diameter and 7.5 cm thick; the faces are quite 
nearly parallel and vertical. The flux density in 
the air gap can easily be adjusted to the desired 
value to within 1.5 gauss by means of a con- 
trolled current in two solenoids around the pole 
pieces. An investigation of the direction of the 
field in the gap showed that it deviates from the 
normal to the pole faces by less than 0° 05’. 
Computation by the series for cos 0° 05’ shows 
that the horizontal component of // differs from 
IT by less than 1 part in 500,000. Thus the force 
which the Cotton balance measures, i.e., the 
horizontal component of //, is sensibly //. 

The field in the air gap could be explored with 
search coils and galvanometers calibrated as 
fluxmeters. J/ as a function of the distance from 
the center of the field was determined with a 
coil and galvanometer designed to give an ac- 
curacy of one percent. Small variations of J/ 
near the center of the field were obtained with a 
1000 turn coil of diameter 0.9 cm and an 18 ohm, 
high sensitivity galvanometer, the combination 
having a sensitivity of 0.0241 gauss per cm. 

Supported above the air gap was a new Becker 
precision balance of which one pan had been 
removed. In place of the pan was hung a Pyrex 
glass plate about 30 cm long and 3 cm wide with 
its longer edges optically worked to be parallel 
to within 1/500 mm for 15 cm from one end. 
This plate hung from the lower ends of two 
parallel aluminum strips which hung in turn 
from a support-assembly attached to a knife edge 
of the balance. The plate-assembly was free to 
swing in two dimensions and was protected from 
external air currents by a surrounding brass case. 
Two straight silver strips about 1/10 mm thick 
were cemented with very dilute white shellac 
around the lower end and the two long edges of 
the plate; they were joined by a flat fold at the 
middle of the lower edge and formed the con- 
ductor of currént in the magnetic field. This 
lower edge was placed along the x axis of Fig. 1 
while measurements of /7 were in progress; the 
center line of the plate was along the y axis. 

Measurements of plate and strips could be 
made with a Brown and Sharpe gauge which was 





Fic. 1. PP’ are photographic plates; O is the source of 
8-particles. 


graduated in 10~*’s of an inch and which could be 
reliably estimated to a 10° of an inch, with a 
Gaertner comparator and invar steel standard 
scale 10 cm long, and with a new Brown and 
Sharpe 31 mm standard disk. The invar scale 
was stated by the Bureau of Standards to be 
accurate to 1 part in 10°. All measurements of 
the plate, strips, and photographic plates are ex- 
pressed in terms of this standard scale when used 
as absolute lengths. 

When the current in the strips was reversed, 
the change in force on them was measured by 
balancing it with a mass m on the balance. 
m was nearly 400 mg and was obtained in the 
main by combining the 200, 100 and 100 mg 
weights of a new and previously unused set of 
platinum, assayers’ weights by Becker; each of 
these weights is quoted by the manufacturer as 
being accurate to within 0.005 mg. To get m 
really accurately, small balance deflections (oc- 
curring when the current in the strips was re- 
versed) could be measured by the motion of a 
spot of light reflected from a small mirror on the 
balance beam to a scale about two meters from 
the balance. The location of the spot of light 
could be read to within 1/5 mm with a magnify- 
ing glass. When fully loaded in the experiment, 
the balance had a sensitivity of 0.409 mg/cm, 
which was entirely sufficient to get m accurately. 
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Current was supplied to the strips around the 
plate by two 160 ampere-hour lead storage bat- 
teries connected in parallel. The current was sent 
through two standard resistors of 4.0001 inter- 
national ohms at 25°C and having negligible 
temperature 25°C. These re- 
sistors were connected in parallel by heavy cop- 
per bars and the current through them adjusted 
until a galvanometer in series with a standard 


corrections near 


cell (the two being connected across the resistors) 
showed no deflection. A measurement of // using 
one resistor balanced by two standard cells 
vielded the same value of // as this arrangement, 
so the resistance between the resistors and the 
copper bars was entirely negligible. The resistors 
were calibrated to 0.005'percent by the Bureau of 
Standards in February, 1934, and, except for 
Scott's work then, have not been used until the 
present experiment. The standard cell along with 
two others was calibrated to 0.01 percent by the 
Bureau of Standards at the same time and was 
found to have an e.m.f. of 1.01884 international 
25°C; 


volts at 
these cells had not changed their relative e.m.f.’s 
by more than }X10~° volts, so it was assumed 
that their absolute had altered 
appreciably. The galvanometer used to indicate 
the the current was of 
sensitivity 1.8410-° amp. cm, readable to 1/5 
mm easily. The current I could be adjusted by 
means of a variable resistance network to within 
}Xmm scale deflection, or to within 10~° amp. 
The Ra B £6-particles (from the active deposit 
from the radon of 23 mg of radium) were emitted 


at the time of this experiment 


values not 


existence of correct 


from the surface of a platinum wire placed at 0 
(Fig. 1); they traveled in opposite directions as 
allowed by slits in the camera and fell on photo- 
graphic plates PP’. The distances between the 
outer edges of the three sets of clearly defined 
“images’’ on the plates could be measured with 
the Gaertner comparator and standard scale. 
The camera and its enclosing box are of brass 
tested previously for magnetic effects and found 
to have none; Scott fully describes them in his 
paper.’ The vacuum in the box was got by a 
Cenco Hyvac pump running continuously during 
the three and one-half hour exposures. 

All observations and measurements were made 


in two rooms which were maintained at tem- 


peratures between 23° and 25°C. 
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THEORY 


Preliminary measurements have shown that 
when the glass plate is in position for field 
measurements, the center of the silver strip 
around the glass plate lies, in the large, along a 
curve C of the form (exaggerated in the diagram) 
of abcdefgh in Fig. 1; this curve is sufficiently 
nearly symmetric to be treated for our purposes 
as symmetric. In Fig. 1 the magnetic field J/ is 
perpendicular to the plane of the paper and 
directed into the paper ; x and y axes perpendicu- 
lar to // are chosen with their origin in line with 
the centers of the pole faces and half way between 
them, the x axis horizontally ; ds with x compo- 
nent dx is an element of length along C. The 
width / of the current loop is measured at some 
convenient location (as bg) along the glass 
plate. 

To get the field //) at the origin by the Cotton 
balance method, note that the force acting on 
the current J (e.m.u.) in the strip is, besides 
gravity, 1 f, HIdx; when J is reversed, the change 
in force on the strip is measured by the weight 
of the mass m which will just counterbalance it; 


then 


mg=2I | IIdx. 


Since JJ and C are necessarily continuous, we 
can write 


f-fefefe pes 


and since subsequent experiment shows that 


ad . 


fife fe fic, 


a c e J d 


approximate evaluation of these integrals is 


satisfactory. Let 


IT=H1)+h(x) along de, 
Il =H F(y) along ac and fh, 
IT=ITy along cd and ef; 
then 
mg = 21} Ho | Fi ydx+Ho | dx 


wh 


+ | [1y+h(x) ldx+I1y dx+Ho| Fiy)dx} 


“d e I 
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Also let C be represented from a to c by x=f(y), 
and let 


ad 
dx=}Al; 


¢ 


then in view of the symmetry of C with respect 


to the y axis, 


mg = 2 211 F(y)f' (y)dy+HoAl 


+Hdl+ | hiss} 


9 


By solving for the coefficient of J, 


lf mg v! 
lem | ~ 211, F(y) f’ty)dy 
2] ae 


— I]],Al— | had) 


eo ) 


here m, g, 2 and J are measurable with high ac- 
curacy ; the last three terms in the right member 


1 
l 


which is the expression adopted for J/) to be 
used in J/p. 

To get p, the distance d between the outer, 
clearly defined edges of corresponding images on 
the photographic plates is measured; from it is 
subtracted the diameter D of the platinum wire 
supporting the source at O. Then 


p=(d—D)/4. (2) 


I1¢ for a B-particle energy is the product of the 
appropriate value of expressions (1) and (2). 

The relative error to be expected in a single 
determination of //) follows from two principles 
concerning errors: the relative error in a 
product of quantities having known relative 
errors is the sum of the relative errors of the 
multiplicands, and the absolute error of a sum 
of quantities having known absolute errors is the 


* Hartree, Proc. Camb. Phil. Soc. 21, 746 (1923). 
5Gibbs, Adjustment of Errors in Practical Science 
(Oxford Press, 1929), pp. 76-78. 


are of the nature of corrections to the first and 
(as previously stated) are fairly small, so in them 
only the approximate value of /7) need be used. 
To this Cotton balance equation for /7) must 

be added for this experiment the Hartree cor- 
rection for non-uniformity of /7 along the 8-par- 
ticle paths. In Fig. 1 two semi-circular paths for 
the same §-particle energy are Okj and Omn; 
measuring @ along them as shown in the figure, 
let 

H = Hy+h,(@) along Okj, 

I= Hy+h2(@) along Omn. 


Then according to Hartree,* the first order cor- 
rection to be applied to //) to give the constant 
value of H which, if it were acting all along the 
B-particle path Okj, would produce the ob- 
served deviation Oj is 


rf h,(@) sin 6d0. 
0 


Applying to the above equation for //) the mean 
of the two Hartree corrections for the two paths 
of the same energy of 8-particle, we get 


meg w of /2 1 * ad 
ITy= —2Ho| F(y)f’ (y)dy — Hy Al — h(x + ( f h,(@) sin 6d0+ he(@) sin wid), (1) 
21 ve “—i/2 4 “e vo 


sum of the absolute errors of the summands. 
Let e;, 7=1 through 10, be the absolute errors 
which might possibly be present in single de- 
terminations of d, D, 1, m, g, I, Ilo JS Ff'dy, IvAl, 
Sh(x)dx, and Sh(@) sin 0@d@ respectively ; these 
must be estimated from the precisions of the 
observations on the several quantities entering 
into expressions (1) and (2). Then since (by the 
definition of relative error) there corresponds to 
the relative error in a quantity an absolute error 
whose magnitude is equal to the product of the 
magnitude of the quantity into the relative error, 


(mg 21) (€4 m+ és ete IT) +2eztestes 


is the absolute error to be expected in the brack- 
eted part of expression (1); to it corresponds 
quite nearly the relative error 


e,/m+e5/g+eg/I+(2er+estes)/(mg,2/). 


Hence the absolute error in the first term of the 
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right member of (1) is approximately 


€3 €4 €5 €—§ 2Z€7+egstes 
11o( + + + + sewrrees. ) ; 
Lm g I mg/2I 


added to the absolute error 
3( Sé10t+  €10) = €10 a 


in the second term and divided by //o, we have as 
the relative error to be expected in a single 
determination of //9: 

€ €6 Zeztestes €10 


3 €4 €5 
— pp ' . 
; mg f mg/2I 2H 


e.= 

The relative error to be expected in a single 
determination of ¢ is obviously 

€2=(€:+€2)/(d—D). (4) 

If €1, 


tively as the maximum absolute errors to be ex- 


+= 1 through 10, are estimated conserva- 


pected in their parent magnitudes, the corre- 
sponding maximum relative error to be expected 
in an accepted //p is 

e3=Ci ter. (5) 

To get the energy £E of a 8-particle of known 
IIlp, combine with the relativistic expression 

E=m,c(1/(1—v?/c?)!—1] 
for E the relativistic force equation 
Hev = mv" /p = myv?/p(1 —v?/c?)! 


so as to eliminate v. Here v is the velocity of the 
8-particles, mo their rest.mass, e their charge, and 
c the velocity of light. The resulting expression 
E=myoc[(A+1)!—1], A =(1pe/myo)*(1/c)?, 

is the one commonly used. If e/mp is in e.m.u. 
and if e is in e.s.u., then m=(mo/e)(e/c); also if 
E is expressed (as V’) in electron volts, its value 
in ergs must be multiplied by 10~* c/e; thus if 
IIp is in gauss cm, the energy in electron volts is 

v=(mo/e)?[(A+1)'—1]X10-. (6) 
This can be put more conveniently for computa- 
tion if we define ¢ by 

@=tan-! A!=tan (Hpe/my)(1/c), 

then 


v= (Mo e)c?(sec 6—1) K10 ". (7) 
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If €11 and €i2 are the absolute errors to be ex- 
pected in e/mo and c respectively, the relative 
error to be expected in (A +1) (see Eq. (6)) is 


A 2€11 2€12 
(20.4 + ). 
A+1 e/My Cc 


Since the extraction of the square root of a quan- 
tity with a known relative error yields roots with 
relative errors only half that of their square, the 
absolute error in (A +1)! is 


A €11 €i2 
(a+ + ). 
(A+1)3 e/mo C 


which is also to be expected in [(A+1)!—1]. 
Hence by addition of the relative errors in 
mo/e, c, and [(A+1)!—1], the relative error to 


be expected in V is 


€11 2€12 


+- 


€; Mo c 


A €11 €12 
+ (« + ; + ). (8) 
(A+1)—(A+1):} e/my ¢ 


The values of €:; and e€;2 are preferably the 


eg 


maximum errors judged to be possible in the ac- 
cepted values of their parent magnitudes; if 
they are not, but are the probable errors quoted 
with the accepted values, e, is of the nature of 
an hybrid which is only a lower bound of the 
maximum relative error to be expected in V. 

It should be remembered that expressions (3), 
(4), (5) and (8) do not take into account any 
systematic errors which might have been pre- 
sent and undetected in the determinations of the 
quantities appearing in expressions (1), (2) and 
(/). 


MEASUREMENTS 


The experiment consisted in measuring the 
quantities appearing in expressions (1) and (2) 
so as to determine, respectively, 7 and the p's 


‘For if a is a quantity with a small known absolute 
error Aa, then 
(a+ Aa)!=a!+ Aa/2a!+---=a!+ Aa/2a}. 
Thus the relative error in a} is 
(1/a4) x (Aa/2a*) =1/2(Aa/a), 


or is one-half the relative error in a. 
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for the three lines. J/) was determined from forty 
determinations of m and each p was determined 
seven times in order to reduce the maximum 
possible relative error to be expected in //p as far 
as is reasonably possible. The accepted //p for an 
energy is the product of the //) value (from the 
forty determinations of m) as corrected for the 
proper path by the Hartree correction, into the 
appropriate radius of curvature of the path. 
The correction quantities were in general de- 
termined first. The error quantities ¢; were 
estimated as each quantity was measured. All 
measurements were repeated at least once, 
usually more than once. Data are usually recorded 
with more figures than are strictly significant, 
in order to prevent the entrance of errors in 
computation. 

The functions h(x) and h(@) were determined 
from changes in J/ revealed by the 1000 turn 
coil and the sensitive galvanometer. Table I 
is the data obtained for /h(x)dx at H7)>=1250 
gauss. Integration was effected numerically by 
assuming that A(x) is linear between adjacent 
experimentally determined points; the result is 
that 

ali2 
h(x)dx = —0.319 gauss cm. 


be 


i/é 


It is estimated that «=0.01. Table II is the 
data used in calculating /h(@) sin 6d6 for 
I7y>=%1250 gauss and p=1.54 cm; A(@) are given 
in gauss. Integration is exactly as before and the 
results for this radius and for the two other radii 
used are tabulated in Table III; the Hartree 
correction as tabulated is the mean of the two 
corrections for the two paths from 0. It is esti- 
mated from experience with these that €;9= 0.01. 

The value of / was measured at the location 
bg (Fig. 1), which is between y=8.0 and y=9.0. 
It was measured first with the comparator and 
standard scale and found to be 3.07918 cm. 
After the experiment was terminated, / was 
measured again, but with the 31 mm disk and 
the Brown and Sharpe gauge. The value of 
3.07891 cm was obtained in this way. The finally 


accepted value of / is the mean of these two: 
/=3.07905 cm. These values of ] are got by 
subtracting from the overall width of the glass 
plate and silver strips at bg, half the thicknesses 


TABLE I. Data for fh(x)dx. 


x h(x) 
(cm) gauss) 
—1.5 —0.018 
—1.0 —0.106 
—0.5 — 0.076 

0.0 —(0).000 

0.5 —0.113 

1.0 —0.187 

ta — 0.263 


TABLE II. Data for Sh 0) sin 6d0, p=1.54 cm. 
p=1.54 cm 

0 hy(@) he( 0) 
0 0.000 0.000 
30 —().027 —0.263 
60° —0.115 — 0.566 
90 —0.110 — 0.667 
120 —0.234 —().501 
150 —0.130 —0.314 
180 — 0.189 —0.150 


TABLE III. Hartree corrections. 


p s Sh; sin 6d0 | } fh sin 0d0 Correction 
cm) gauss) gauss) (gauss) 
1.12 — 0.098 —0.352 —().225 
1.34 —0.128 —00.418 —0.273 
1.54 —0.122 —().471 —().296 


of the silver strips at that location. To accom- 
pany / we assign ¢,=0.00015. 

To evaluate f’dy and A/, the total width of the 
glass plate and of the silver strips in place was 
measured at many locations along the plate with 
the 10-* inch gauge; as above, from these widths 
were subtracted half the thicknesses at the same 
locations of the silver strips. In Table IV are 
representative data from the large amount 
actually obtained by these measurements. As 
before it is assumed that f(y) is linear between 
points determined experimentally. Since / was 
chosen as the width of the current loop between 
y=8.0 and y=9.0, Al is obtained by subtracting 
the width of the current loop at y=8.5 from 
the width at y=0.0; from Table IV we have 
Al = 3.08251 — 3.07905 = 0.00346 cm. 

The function F(y) was obtained with the less 
sensitive fluxmeter arrangement; its values are 
given in Table V. For purposes of integration, 
F(y) was also assumed to be linear between ob- 
served points. It is found that 2/Ff’dy= 
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TABLE IV. Width of current loop 


Width of strips | } sum of strip | 2f(y) =width of 


y and plate thicknesses current loop 
cm) (cm) cm) (cm) 
0.0 3.09195 0.00944 3.08251 
0.5 3.08914 0.00945 3.07969 
1.0 3.08901 0.00945 3.07956 
2.0 3.08853 0.00947 3.07906 
23 3.08835 0.00948 3.07887 
2.4 3.08876 0.00949 3.07927 
2.8 3.08905 0.00950 3.07955 
2.9 3.08862 0.00951 3.07911 
3.0 3.08849 0.00952 3.07897 
4.3 3.08966 0.00952 3.08014 
4.5 3.08953 0.00953 3.08000 
5.5 3.08878 0.00952 3.07926 
6.0 3.08877 0.00951 3.07926 
6.7 3.08903 0.00949 3.07954 
7.0 3.08911 0.00950 3.07961 
8.5 3.08853 0.00948 3.07905 
9.5 3.08837 0.00946 3.07891 
11.5 3.08861 0.00950 3.07911 
12.0 3.08873 0.00949 3.07924 
12.4 3.08854 0.00948 3.07906 
12.6 3.08876 0.00946 3.07932 
13.2 3.08905 0.00944 3.07963 
13.5 3.08861 0.00944 3.07917 
15.0 3.08861 0.00945 3.07916 
17.0 3.08869 0.00945 3.07924 
17.5 3.08903 0.00943 3.07960 
18.0 3.08854 0.00943 3.07911 
0.0 3.08847 0.00942 3.07905 
25.0 3.08872 0.00944 3.07928 


—().00354 using the entire set of data of which 
Table IV is a representative portion; the Table 
IV data give a slightly lower value than this. 

The actual values of the corrections //)Al and 
211, Sf Ff'dy are obtained directly as 


I, Al = 1250 X 0.00346 = 4.33 gauss cm, 
2H S Ff'dy = 1250 X ( — 0.00354) 
= —4.43 gauss cm. 


It is estimated that ¢ =«3=0.05. 

the in getting J are 
calibrated to only 0.005 percent and since the 
standard cells, to only 0.01 percent, the maximum 
error to be expected in J is at least 0.015 percent. 
The error inherent in adjusting for J with the 
resistance network is at most 1077 e.m.u. By 
Ohm's law J is 0.509403 international ampere ; 
the conversion’ to e.m.u. gives J=0.0509433 
e.m.u. Thus ¢; = 0.0000076+0.0000025 = 0.00001 ; 
the second term is added because of the un- 
certainty in the factor of conversion from inter- 


Since resistors used 


national to electromagnetic units. 


? The Smithsonian Physical Tables, eighth edition (1933), 
p. 81, gives one international ampere equal to (0.99995 
+0.00005) absolute ampere. 
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Because three weights are required to produce 
the 400 mg of m, e, must contain 0.015 mg. 
Because of errors possible in locating the ends of 
the balance deflections when J is reversed, which 
come to 2X0.02 0.41 =0.0164 mg, e, is larger 
than 0.015. Since m was measured nine times 
in each of the forty determinations of it, we 
take® e,= (0.015+0.0164 /9!) /40!=0.0032 mg. 

The value of g in the locality of the experiment 
is 979.28 cm/sec.” and ¢,=0.005 cm/sec.’. 

By observing the balance deflections when / 
was reversed and when the 400 mg were placed 
on the balance, the correction to be applied to 
the 400 mg to give m accurately could be found. 
These deflections were observed nine times for 
each determination of m. The forty mean ob- 
servations of the balance deflection are recorded 
with their mean in Table VI. The evaluation of 
II) by taking for m the mean of a large number 


TABLE V. Function F(y). 


¥ 
(cm) F(y) 
0.0 1.00 
5.6 1.00 
7.5 0.91 
8.5 0.50 
9.5 0.26 
10.5 0.17 
11.5 0.13 
12.5 0.11 
15.0 0.07 
17.5 0.04 
22.5 0.02 


TABLE VI. Observations of balance deflections for m. 


DEFI DEF! DeFl Der! 
TRIAl cm) TRIAL (cm) TRIAI cm TRIAL m 
1 1.53 11 1.52 21 1.52 31 1.52 
2 1.52 12 1.54 22 1.50 32 1.52 
3 1.52 13 1.53 23 1.52 33 1.52 
4 1.53 14 1.54 24 1.51 34 1.52 
5 1.52 15 1.54 25 1.51 35 1.51 
6 1.53 16 1.50 26 1.51 36 1.51 
7 1.51 17 1.51 27 1.52 37 1.54 
8 1.52 18 1.51 28 1.53 38 1.51 
9 1.53 19 1.52 29 1.52 39 1.52 
10 1.54 20 1.52 30 1.52 40 1.50 


Mean = 1.520 


* Analogously to the customary procedure used in the 
operation of classical error computations, we assume that 
the relative error of the mean of m determinations of a 
quantity varies inversely as nm. See, e.g., Bartlett, The 
Method of Least Squares, third edition (1915), p. 44, 
Eq. (46). ° 
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of observations which extended over several 
days is valid because Scott*® showed in his experi- 
ment that over a period of weeks the value of [/» 
remained sensibly constant. The value of m to 





1 0.400622 X 979.28 
Hy= ( 


3.07905 2 0.0509433 

Values of d were read from the photographic 
plates with the Gaertner comparator corrected 
by the standard scale; they could be repeated to 
within 0.0001 cm usually. Since each was done 
seven separate times, ¢;=0.0001/7!=0.000038. 
The value of D was found with a direct-reading 
10-4-inch Brown and Sharpe micrometer caliper 
and is 0.02405 cm to within e.=0.0001 cm. In 
Table VII are presented the seven values of each 
p got by Eq. (2) and their three means. Since // 
is constant, these should be constant within the 
proper groups in the table. 

Knowledge of the 7) from Table VI, of the 
mean p’s of Table VII, of the Hartree corrections 
of Table III, of the various correction terms, 
and of the estimates of the possible errors is 
sufficient to provide the desired values of //p and 
of the possible errors present. 


RESULTS 


Tables VIII, LX, X, and XI contain the results 
of the experiment. The accepted /Hp’s of Table 


TABLE VII. 


p Mean p 
Trial (cm) (cm) 


1.12442 
1.12446 
1.12442 
1.12440 1.12440 
1.12442 
1.12440 
1.12426 


SOUS wr— 


1.33642 
1.33642 
1.33650 
1.33645 1.33644 
1.33657 
1.33635 
1.33642 
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1.54471 

1.54476 

1.54468 

1.54462 1.544 
1 

1 

1 
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54478 
.54470 
.54471 
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be used is thus 0.4000000+1.52 0.000409 
= (0.400622 gm. By use of Eq. (1) and the pre- 
viously evaluated correction quantities it is seen 
that, except for the Hartree corrections, 


+4.43— 33+0.32) = 1250.70 gauss. 


VIII are from the value of //) from Table VI 
and corrected by Table III and from the p’s of 
Table VII. In view of the good agreement (1 
part in 6400) of the present value of //p for the 
most intense 8-particles and that of Scott, it 
seems safe to consider the //p’s listed in Table 
VIII as sufficiently adequately determined in 
absolute value (within the possible limits of 
relative error given below) to serve as rather 
reliable standards of 8-particle J/p’s. 

In Table LX are the possible relative errors 
necessary to the computations based on Eqs. 
(3), (4) and (5) for the least //p; they are col- 
lected from the section on Measurements. Thus 
the values of J/p in Table VIII cannot be con- 
sidered as any more accurate than 1 part in 
3000, even though each measurement used in 
getting //p is more accurate than 1 part in 3000. 
Evidently some 60 percent of the uncertainty is 
accounted for by the uncertainty in the calibra- 
tions of the standard cells and standard resistors 
and in the uncertainty in the conversion factor 
for international and electromagnetic units of 
current. Within this limit of 1 part in 3000 the 
present experiment yields the same value of //p 
for the most intense Ra B £-particles as Scott's 
experiment, though it is seen that Scott's hopes 
for an accuracy of 1 part in 10,000 are rather too 
optimistic. 

The energies corresponding to the //p’s we have 
established are obtained from Eq. (7) using® 


TABLE VIII. Hop’s of three intense Ra B B-ray lines. 


Hp 
Line (gauss cm) 
Least intense 1406.0 
Of medium intensity 1671.1 
Most intense 1931.5 


®The value of e/mo is that given by Birge, Phys. 
Rev. 49, 204 (1936). c is taken from The Smithsonian 
Physical Tables, eighth edition (1933), p. 74. 
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TABLE IX. Maximum relative errors to be expected in least Hp. 


Magnitude 
parts per 


Error Magnitude 10,000) 
€;/l 0.00015 /3.079 0.49 
e,/m 0.00002 /.4007 0.08 
5/2 0.005 /979.3 0.05 
e6/I .00001 /.05094 1.97 
21 (2ert+ estes) (mg 0.16/3850 0.41 
€10/2 0.01 /2500 0.04 
e,/(d—D) .000038 /4.496 0.08 
eo/(d—D) .0001 /4.496 0.22 
Sum =e; 3.34 


€/ My= 1.7562 X10? e.m.u./g, 
c= 2.99796 X 10" cm/sec. 


with €,,;=0.00027 X10? and €;.=0.00004 X10". 
They are in Table X. 

In Table XI are the numbers necessary for the 
calculation of the lower bound e; of the relative 
error to be expected in the least V according to 
Eq. (8). They follow from the preceding two 
paragraphs. Eq. (8) becomes by substitution of 
these numbers 


1.54 0.24 
&= + 
10,000 10,000 
0.678 / 3.34 1.54 0.24 1.1 
+—_( +) a 
0.383 (10,000 10,000 10,000 1000 


Thus the energies V cannot be considered as more 
accurate than 1 part in 1000. In fact, since e€;, 
and €,2 are only classical probable errors, e, may 
be even larger than 1 part in 1000, the increase 
depending only upon how much larger than the 
the absolute 


errors possible in e/my and ¢ are. 


probable — errors corresponding 
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TABLE X. Energies of three intense Ra B B-rays 
Hp V 
gauss cm) electron volts) 
1406.0 1.512 «10 
1671.1 2.044 « 105 
1931.5 2.610 x 108 


TABLE XI. Quantities for calculation of es for least V 


Magnitude 
parts per 
10,000) 


Quantity Magnitude 


€1:/(e/my) 0.00027 / 1.7562 1.54 
€12/¢ 0.00004 (2.998 0.12 
e 3.34/ 10,000 3.34 
1.406 < 1.756\? . 
= vf 
4 ( 2.998 ) 20.019 


It is also noteworthy that even though IIp is 
accurate to within 1 part in 3000, the major part 
of the 1 part in 1000 accuracy allowed for V' is 
accounted by possible errors in J/p. 

Finally, it should be noted that the improve- 
ments of this experiment over that of Scott lie 
in the facts that: a new balance of higher sensi- 
tivity, a new set of weights for m, straight silver 
strips, the standard scale, and a new gauge 
graduated in 10~*’s of an inch were available ; 
twice as large a current, determined by a far 
simpler potentiometer, was used in the Cotton 
balance ; it includes measurements for the three 
most intense Ra B lines rather than for only one; 
and it provides a conservative, explicit estimate 
of the precision of which it is capable, 1 part 
in 3000 in //p, 1 in 1000 in V. 

The writer also wishes in conclusion to record 
his appreciation both for the suggesting of this 
experiment and for the valuable aid given him 
during the experiment, by Dr. H. A. Wilson. 
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Tests of the Validity of X-Ray Crystal Methods of Determining e 
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(Received June 15, 1936) 


The procedure in determining the electronic charge, e, 
from measurements of x-ray diffraction angles and the 
density of crystals has been criticized on the ground that 
while the density measurements are averages of a large 
volume, say 10 cc of crystal, the x-rays are diffracted from 
a surface layer as determined by ‘‘extinction’’ only about 
4 10-4 cm thick, a very small fraction of the total volume 
in an entirely specialized location. Partial answers to this 
criticism by Y. C. Tu and by Allison and Armstrong have 
rendered its validity much less probable. The present work 
with extremely fine powdered calcite samples whose 
density was determined with a pycnometer closes the ques- 
tion completely, we believe, leaving no ground whatever 
for the criticism and placing the x-ray determinations of 


I. CRITICISMS OF THE X-RAY METHOD 
OF DETERMINING e 


T now seems highly probable that the excellent 
and concordant ruled grating measurements! 
of the wave-lengths of a number of cha.’ >teristic 
x-ray lines establish in absolute centime anits 
a reliable scale of wave-length values which 
combined with the very precise relative scale of 
x-ray wave-lengths in terms of Siegbahn’s con- 
ventional unit fixes the absolute values of all of 
the carefully measured x-ray lines to a precision 
of say one part in two thousand. The usual pro- 
cedure in determining the electronic charge e then 
involves a measurement of the diffraction angle @ 
in Bragg reflection from a crystal surface and a 
computation of the interplanar spacing d of the 
crystal lattice in absolute units by means of the 
familiar formula (with refractive index correction 
of Siegbahn) 


n\ = 2d sin 6(1—6/sin? @). (1) 


By such methods the absolute dimensions of the 
crystal lattice may be determined. If, as in the 
case of calcite, the lattice is not cubic the absolute 
volume of the unit cell may be determined either 
(1) by combining the value of d with x-ray or 


1 The following is far from an exhaustive list of such 
work: Compton and Doan, Proc. Nat. Acad. Sci. 11, 598 
(1926); J. Thibaud, J. de phys. et rad. 8, 13 (1927); E. 
Backlin, Inaugural dissertation (Uppsala, 1928), Zeits. f. 
Physik 93, 450 (1933); J. M. Cork, Phys. Rev. 35, 1456 
(1930); J. H. Bearden, Phys. Rev. 37, 1210 (1931). 


eon a very solid experimental foundation. Our work yields 
a value for e of 4.799+0.007 abs. e.s.u. based on the scale of 
ruled grating wave-lengths in which Cu Ka has a wave- 
length (Bearden’s value) of 1.540610~-§ cm. Efforts to 
remove the above objection by studies on macroscopic 
calcite samples which failed to withstand our own criti 
cisms are described since they emphasize the pitfalls that 
a conclusive test must avoid. Several new unexplained 
crystal diffraction effects discovered in this search for a 
conclusive test are also described. Agreement to one part in 
six thousand between Bragg and Laue diffraction angles 
was established for two different sets of Miller indices in 


macroscopic calcite. 


optical goniometric measurements on a macro- 
scopic crystal to establish the characteristic 
angle 8 of the rhombohedron or (2) x-ray de- 
terminations of the interplanar spacing may be 
made for several Miller indices, two being theo- 
retically sufficient to establish the necessary 
angular relationships. The first mentioned pro- 
cedure is by far the more common. The absolute 
volume of the unit cell for calcite generally ex- 
pressed as 8d*#(8) is then combined with a 
measurement of the macroscopic density p of the 
crystal employed and by means of the formula 


d=(Me)'{2Qp%(8) ]}-! or e=2d°Qpb(8)M- (2) 
J ¢ 


the value of e is computed. M is the molecular 
weight of CaCO ;; and Q the Faraday constant. 


Three criticisms enumerated 

(1) Criticisms have been raised from time to 
time that the macroscopic density which is an 
average value for a crystal volume of, say, 10 cc 
may not apply to the extremely thin superficial 
layer (5X10-* cm thick in the first order) in- 
volved in Bragg reflection. Allison and Arm- 
strong”? have endeavored to answer this criticism 
by measuring the reflection angles in different 
orders up to the fifth where the penetration, they 


2 Allison and Armstrong, Phys. Rev. 26, 701 (1925). 
See also Allison’s excellent discussion in Compton and 
Allison, X-Rays in Theory and Experiment, page 698. For 
experimental support of the shallow penetration depths 
given by extinction formulae see recent work by W. Kossel 
and H. Voges, Ann. d. Physik 23, 677 (1935). 
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compute, should be of the order of 1.2 10-* or 
40,000 planes. Y. C. Tu® finds that calcite, rock- 
salt, KCl and diamond agree very well in giving 
the same value of e. Criticisms of the x-ray 
method of determining e nevertheless have con- 
tinued on the above grounds and it must be 
admitted that the superficial layers of a crystal 
constitute such a specialized part of its volume 
that every effort should be made to extend the 
X-ray measurements so as to make the latter as 
representative as possible of the entire volume 
whose density has been measured. Unless this is 
done there is certainly a logical hiatus in the 
above reasoning by which we arrive at the value 
of e and it is the closing of this gap by a solid 
experimental procedure which we here chiefly 
wish to report. 

(2) It is not out of place to mention that the 
methods we have employed also conclusively 
answer the criticism sometimes raised that the 
crystal may contain minute cracks, gaps or 
crevices in its structure. 

(3) Another type of criticism which we have 
never heard raised but which we believe might 
have occurred to the truly critical minded is a 
question as to the precise validity of the Eq. (1) 
relating 6, \ and d, when the variables are re- 
garded as having an absolute significance. A 
really critical defense of this equation requires an 
appeal to the dynamical theory of x-ray diffrac- 
tion. Darwin, Ewald, Prins, and Laue have 
all contributed brilliantly to this subject. Never- 
theless the dynamical theory for the space grat- 
ing, taking account of mutual to and fro inter- 
action between all of its parts, becomes so com- 
plicated that none of the treatments can quite be 
regarded as straightforward rigorous solutions of 
the wave equation for the case of a medium with 
a threefold periodic distribution of refractive 
index with the same beautiful generality and 
finality of the corresponding treatment by 
Eckart* on diffraction from ruled plane gratings. 
The effect of extinction about which there is yet 
very little good experimental information prob- 
ably operates to produce very rapidly changing 
intensities and amplitudes of the electromagnetic 
radiation from point to point within the crystal. 
Such an assumption as the immediate reorganiza- 


‘ 


Y. C. Tu, Phys. Rev. 40, 662 (1932). 
C. Eckart, Phys. Rev. 44, 12 (1933). 
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tion of secondary dipole wavelets to form plane 
waves inside the crystal, so essential to Ewald’s 
theory, or the assumption of constancy of phase 
displacement for all atoms might not be quite 
correct for instance. The precision measurements 
in the third order here described using calcite 
powdered to the order of 10-* cm and smaller 
seems to us to answer completely any objections 
founded on doubts as to the rigor of the dy- 
namical theory since the simple elementary 
theory of Laue or Bragg must certainly apply to 
crystals so small that the radiation is negligibly 
extinguished in passing through them. 


II. SEARCH FOR THE Most CONCLUSIVE METHOD 


The idea of the powdered calcite method of 
answering the above objections only came to us 
as an evolution from two other methods which 
after much research we felt did not meet critical 
objections which we ourselves could raise. It is 
probably valuable therefore to describe as briefly 
as possible the pitfalls discovered in this way the 
better te emphasize what must be avoided if a 
truly r able conclusion is to be drawn. Further- 
mo. ‘ran across certain interesting new and as 
yet unexplained phenomena, well worth report- 
ing for themselves, in this search for a conclusive 
method. 


Two crystal spectrometer studies 

Our first attempt consisted in a comparison, by 
means of the two crystal spectrometer, of the 
calcite reflection angles for the two cases which 
we called briefly Bragg reflection and Laue re- 
Hection. We reasoned that in the Laue case with 
the incident beam entering one face and the re- 
flected beam emerging from an opposite face 
after reflections by internal planes parallel to 
neither of these faces the whole interior structure 
of the crystal must have served for reflection of 
the radiation and the reflection angle should 
therefore be representative of the entire crystal 
just as was the mean density measurement. This 
method of course fails to meet criticisms (2) and 
(3) above however. 

On account of the construction of the only two 
crystal spectrometer then available,’ not the 
entire Bragg or Laue angle @ was measured but 


>]. W. M. DuMond and A. Hoyt, Phys. Rev. 36, 1702 
1930). 
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Fic. 1. Line structures obtained with etched calcite. 


only the angular interval between two widely 
separated lines of the molybdenum K spectrum. 
Crystal No. 1 of this spectrometer, a large and 
very perfect specimen of split calcite obtained 
from Steeg and Reuter in Hamburg was regarded 
as the collimator. Crystal No. 2 was very thin 
and could be oriented so as to reflect the col- 
limated beam from No. 1 either in Bragg reflec- 
tion from its large flat surface or in Laue reflec- 
tion with the beam traversing its thin dimension 
reflected from internal planes parallel to a cleav- 
age surface, in this case the thin edge of the flake. 
All appropriate alignments of the crystals by 
optical reflections from the faces and from the 
thin edge were carefully attended to. Gliding of 
the reflecting region across the crystal faces in the 
exploration of a spectrum was prevented by the 
mechanism of the spectrometer and x-ray table 
support. The effort was to see whether the angular 
interval between two given lines would be iden- 


tically the same for the Bragg case as for the 
Laue. One of the first calcite samples used for test 


was rendered thin by etching down a rectangular 
portion in the center of a thicker (5 mm) slab 
with dilute HCI squirted through a fine nozzle. 
The etching did not proceed very uniformly and 
had to be controlled by stopping off various 
regions from time to time with wax. The thin 
window so formed varied in thickness from 0.5 to 
0.2 mm. The unetched cleavage face on the 
reverse side could be used for Bragg reflection. 

The etched crystal yielded the surprising line 
structures of Fig. Prominent facts and infer 
ences about these curves are: (1) The “‘lines’”’ are 
about 3 times as broad as those usually observed 
with good calcite though the specimen was 
initially of the highest quality. (2) The curious 
rugged profile of each line, suggesting the pres- 
ence of several simple component lines partially 
resolved, is reproduced for quite different wave- 
lengths (Mo Ka; Mo K£,, 3). (3) The invari- 
ability of pattern for different wave-lengths indi- 
cates that these component peaks are caused by a 
variety of orientations in the crystal rather than 
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by a variety of grating constants. (4) A composite 
structure appears even in Bragg reflection for the 


unetched and apparently quite undamaged cleav- 
age face on the opposite side to the etched surface. 
This may be caused either by released strains and 
disorientations of parts of the crystal where 
etching occurred extending themselves elastically 
through the thin dimension of the crystal to the 
unetched face or by the penetration of radiation 
from the unetched face through to the disturbed 
etched region sufficiently to permit considerable 
reflection there. (5) The rocking curves in the 
(1-1) order both in Bragg and Laue reflection 
also exhibit a composite broad structure but the 
structures in Laue and in Bragg reflection for 
(1,1) and (1-1) 
a different characteristic profile. 


orders, four cases in all, have each 


These results which we cannot explain are the 
more surprising since it has recently been shown 
that the rocking and reflection curves can actually 
be markedly narrowed for calcite by appropriate 
light etching. This etched sample was obviously 
useless for our primary purpose. These results 
made us perhaps overcautious in regard to the 
possible existence of internal stresses in calcite 
which may have been released and thus rendered 
observable by the nonuniform etching process. 

The work with the two crystal spectrometer 
was continued with thin cleaved plates of calcite 
for the crystal No. 2. Several good thin samples 
a millimeter or less in thickness were obtained. A 
great deal of experimentation was required to 
work out the best mounting for such thin and 
easily deformed flakes and much trouble was en- 
countered with slight angular shiftings of the 
flake in its mount. In order to get reliable results 
which would average out any slight effect of 
curvature of the flake it was eventually found 
essential to take measurements of the angular 
interval between lines from both sides of the 
crystals both for the Bragg and Laue cases. 
Slight differences of temperature in different 
parts of the spectrometer, caused by a lamp used 
to read the angle scale, were found to cause some 
of the anomalies. Some mysterious erratic shift- 
ings of the line positions with time were never 
completely explained but fortunately disappeared 
in our later and more careful work. Our first 
observation before all of these difficulties had 
been found and corrected happened to show by 
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pure accident a difference between the angular 
intervals for Laue as compared to Bragg reflec- 
tion indicating a difference in grating constant 
for these two cases of 0.2 percent and of the right 
direction to explain the, at that time, current 
magnitude of the discrepancy between x-ray and 
oil drop values of e. The possible importance of 
such a result made it imperative to check the two 
crystal spectrometer work by an independent 
method which would permit of measuring with 
precision the entire diffraction angle in Laue 
reflection. We therefore constructed an especial 
precision photographic spectrometer, to be de- 
scribed below, which yielded Laue diffraction 
angles with identically the same thin calcite 
crystal that had been used on the two crystal 
spectrometer in complete agreement with the 
most reliable determinations of Bragg diffraction 
angles to well within the estimated error of the 
method (one part in six thousand). 

To eliminate all uncertainties we continued the 
two crystal spectrometer comparisons and found 
that with a sufficient number of precautions and 
especially that above mentioned of taking an 
average of the angular interval measurements 
from both faces of the crystal to eliminate 
minute curvature effects we could obtain com- 
plete and satisfactory agreement between the 
Bragg and Laue angular intervals. The slight 
shifting effect of differential absorption across the 
width of a line in Laue reflection was established 
as negligible. Lines a2 and 82 of molybdenum 
were used. An aluminum filter was introduced in 
the case of the a» line as a convenience to make 
the intensity comparable to that of the 82 line 





Fic. 2. Two crystal spectrometer arrangement for compar- 
ing Laue and Bragg diffraction angles. 




















528 |. W. M. DuMOND AND V. I BOLLMAN 
M / 
1 108 Mo a, 
Unfiltered Filtered thru Al 
reg Chmm , Lave 125mm 
AN, 
(] \ 
ie re |) 
3 pA + 
Tt bX > a 4025246" [Lave revettian — a \ 7 
J NR 0 5296 5" Beam enters face 2} - E- = ty 
+ Sh et* FTX 
a ae oe Y—7-—_ _. ae \ oe 
acaground a nO Wo cs a \ am 
/ : LOANGIOMG 
laine elias . -4§—4- — siete 
¥T * [ Am 
~ i+ 40752°45.5' Lae (O0CTIO7 pf. \ 
i Bean ervers facel / * 
z K, : oe \ 
1A ye ™" , / ‘ xy 
ott + * as _ rr a 4 Ab \ a. 
an eee, al P \ —— f 7 5 _s 
Tv + ee j 
— a a i = + +/ | \ \ Ya 
Back * ry 2 = ZAAG Wd 
—————— *- 1 > = \ 
sf = 40=52'47/° Bragg reriettion it me Ke 
P ran "a Bean rettected trom racel w } \ \. 
oY \ 4. / 
/ 
we fit™ x [|e 
a / \ + d t ° DY 
— } d x ot ai ieee $+ ——e 4 ——___ \ yr v 
Kground —/ + vu cane 4 2 = 
—— — f : ; | \ Basgrouwd 
/ = \. 40 = 52° 45.1" Bragg reriection a = 7 
a yo retlected from fae 2 / \ 
4 ; \. 
YY \ : J j . 
——— = + 4 ; ‘ he, 
, gl ya 
FacAground a ential fat 
eHlection angles -/0 Sttonds— Fagrovr) 
| | l | | | l | l l — 





Fic. 3. Angular intervals between Mo Ka, and Mo Kf; lines 


obtained on two crystal spectrometer comparing Laue and 


Bragg cases. 


without changing the electrometer sensitivity. A 
thinner filter was appropriate for the Laue case 
than for the Bragg since the crystal itself is a 
strong filter. Fig. 2 shows the two crystals in 
place in the two crystal spectrometer while Fig. 3 
shows the final results of the angular interval 
comparisons with this method. 

In the Laue reflection work with the two crys- 
tal spectrometer it was found that the symmetry 
and width of the lines depended greatly on the 
region of the crystal used for reflection in the case 
of all samples. Regions too near the point of 
support were always ‘‘bad” but other regions 
with no such easily explained cause also showed 
distortions of line shape. Such distortions often 
seemed present at, as nearly as we could estimate, 
the same general regions in Bragg reflection. A 


suspicion could have appeared at this point that 
our Laue measurements were not as representa- 


tive of the interior portions of the crystal as we at 
first thought. At this point however we merely 
felt the strong need for some statistical method 
which would free the x-ray measurements of the 
individual peculiarities of different crystals and 
of different regions in the same crystal. 


Study of Laue reflection with precision photo- 
graphic transmission spectrometer 


Fig. 4 is a photograph of the precision photo- 
graphic transmission spectrometer. Its chief 
novelty is the crystal and slit mounting. The 
crystal and slit are held on small independently 
movable carriages sliding in ways on the top of a 
brass disk which turns on a well made conical 
pivot whose axis is accurately normal to the plane 
bed plate surface of the instrument. The plate 
holder is carefully made so as to maintain the 
photographic emulsion also accurately normal to 
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the bed plate surface. By means of appropriate 
slides and leveling screws the slit was placed 
accurately parallel to and in coincidence with 
the axis of rotation of the pivot by reversals 
under a microscope with cross hairs and the thin 
crystal flake by means of other slides and leveling 
screws had its vertical thin edge and its flat face 
both made parallel to the pivot axis by optical 
reflections. One of the slides permitted moving 
the crystal flake so that one face stood very close 
indeed to the jaws of the slit. By further optical 
adjustments the pivoted disk was so oriented that 
a plane through the slit parallel to the reflecting 
planes to be used in the crystal flake would strike 
the photographic plate accurately normally and a 
small lead shield with a fine slit was placed at this 
point on the plateholder opening, just where the 
foot of this perpendicular fell on the plate, to 
serve as a fiducial mark. The pivot was firmly 
clamped during x-ray exposure. Initially, spectral 
line photographs were taken with the crystal on 
the x-ray tube side of the slit. By swinging the 
spectrograph with respect to the x-ray tube, Ka 
lines of the Mo spectrum could be formed on the 
plate by Laue reflection both from the right and 
from the left side of the lattice planes inside the 
crystal flake, the transmitted beam being shielded 
in each case from fogging the plate by a movable 
lead slide. The fiducial center mark was also 
recorded on the plate by a short auxiliary ex- 
posure through the slit and crystal with direct 
radiation. In all cases the foot of the perpendicu- 
lar from slit to plate proved to be so accurately 
centered between the right-hand and left-hand 
spectra that no correction for geometrical asym- 
metry would have been significant. The dis- 
tance from slit to plate was carefully measured 
with a cathetometer and found to be 606.5+0.1 
mm. The cathetometer scale was checked on the 
same comparator used for measuring line separa- 
tions on the plate. Fig. 5 shows the geometry of 
the set-up. 

Very beautiful and finely resolved lines of the a 
doublet were obtained on several plates. Comp- 
ton’s precise measurements’ of the Bragg reflec- 
tion angle in calcite with a precision Société 
Génevoise two crystal spectrometer yielded an 
angle of 6° 42’ 36” for the a; line of molybdenum. 
If this same angle holds for Laue reflection the 


*A. H. Compton, Rev. Sci. Inst. 2, 365 (1931). 
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separation 6 of the right- and left-hand lines on 
our plate should be 142.710 mm. On two photo- 
graphic exposures we obtained 
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142.719 mm. 
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These Laue reflection data are in better agree- 
ment with Compton’s Bragg data then the pre- 
cision of our cathetometer measurement war- 
ranted. These cathetometer measurements, good 
to 1 part in 6000, constitute the determining 
element in the accuracy of this method and the 
agreement between Bragg and Laue reflection 
angles was thus established to that relative 
precision. 

This was identically the same crystal flake used 
on the two crystal spectrometer and it was the 
above result that completely proved the unreality 
of the preliminary apparent discrepancy between 
Laue and Bragg angles obtained at first with the 
two crystal method. 

As a further check we also made a measure- 
ment with the photographic transmission spec- 
trometer using (110) planes (these indices are 
referred to cleavage axes). These planes are 
shown in Fig. 6 which is a projection of a unit 
cleavage cell looking along one of its edges. It is 
evident that 

doo /dyj9 = 2 Cos $a. 


Bearden gives a= 105° 3’ 29”. The wider separa- 
tion of the right- and left-hand spectra in this 
case required moving the plateholder to a new 
position closer to the slit. The new slit to plate 
distance (482.55 mm-+0.03) was determined 





Fic. 4. Photograph of precision photographic transmis- 
sion spectrometer. Slit and crystal are here in position of 
Fig. 7 or Fig. 9a. 
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Fic. 5. Geometry of precision photographic transmission 
spectrometer as normally used for determining diffraction 
angles with slit on plate side of crystal. The radiation left 
the nearly horizontal target surface in the x-ray tube at a 
small glancing angle of about 5°. The diameter of the round 
focal spot furnished a ‘sufficiently extended source just to 
permit the simultaneous reflection of the a and 8 doublets 
with careful adjustment. 


as before and the computed separation of the 
lines reflected from (110) based on Compton's 
Bragg (100) reflection data was 6= 138.613 mm. 
On two plates we obtained 


6,= 138.621 
62= 138.614 


bm = 138.617 mm. 


This would furnish a complete check for this 
piece of calcite not only of the lattice parameter 
but of the angles characteristic of the calcite 
rhomb for reflection from internal planes were it 
not for the results of the next experiment. 


Laue reflection probably not a valid test of the 
interior of a crystal 


As a last precaution we thought it would be 
safer to make quite sure that, in the above work 
with the two crystal spectrometer and the 
photographic spectrometer, the reflection was 
indeed coming from the interior of the crystal. 
To our astonishment and chagrin the contrary 
proved true. 

We turned the crystal and slit mounting of 
the precision photographic transmission spec- 
trometer round through 180° so that now the 
crystal was between the slit and the photo- 
graphic plate. The geometry of the new arrange- 
ment appears in Fig. 7. Now instead of a fine line 
on the plate for each wave-length we expected to 
get a broad band of width d; and d_ for the right- 
and left-hand spectra, respectively. The trian- 
gular region OAB in the crystal we thought 
would be a field in which multiple to and fro 
selective reflection would take place so as to fill 
out the spaces between lines A,B, and AB, with 
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more or less uniform distributions of widths d; 
and do. Instead of this lines only appeared at Aj, 
B,, Az and By, with negligible intensity between 
them. Fig. 8 shows an enlargment of a very short 
section of these doubled lines. For simplicity of 
discussion we will here describe the phenomenon 
entirely in terms of the appearance of a single 


monochromatic line, say a.’ 


’ Entry and exit surface Laue reflections not quite 


parallel 

From the thickness of the crystal and the 
angle y Fig. 7 we computed the separations d, 
and dz to be expected as 0.0232 cm and 0.0217 
cm, respectively. The observed separations of the 
line pairs A,B, and A2Bz turned out to be not 
quite equal to these expectations at all distances 
but, to approach these values the closer the plate 
came to the crystal. The beams A ,B,A2B2 which 
one would have expected to be propagated in 
parallel pairs turned out to have the members of 
each pair apparently very slightly convergent. 
This might have been an effect on the plate of the 
slightly increased overlapping of the members of 
a pair at greater distances as the breadths of the 
slightly with increasing 


members increased 
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Fic. 6. Projection of unit cleavage cell of calcite looking 
along a direction parallel to one edge to illustrate the (110) 
planes in relation to the cleavage planes. 





7 The effect of enhanced Laue reflection at a ground or 
otherwise disturbed crystal face is of course well known 
and understood. The crystals here used however had excel- 
lent cleavage surfaces with narrow rocking curves and we 
were quite unprepared for the almost complete absence of 
reflected intensity between A and B. 
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distance. The tendency for doublet separations 
to appear than they should as the 
members overlap more and more is a well-known 


smaller 


illusion. We tend to doubt this explanation and 
to think the convergence real. 


Fluctuations of intensity along the lengths of the 
spectral lines 


An even more surprising peculiarity of the 
spectral lines taken with the crystal in this posi- 
tion and one which we believe has never before 
been reported, is a more or less periodic fluctu- 
ation of the intensity vertically along the length 
of the line. Both these phenomena were observed 
in varying degrees in all the thin flake samples of 
calcite we examined, some of which had very 
marked fluctuations of intensity amounting to 
complete extinction of the line at roughly uni- 
formly spaced intervals. A careful study of these 
fluctuations at different slit-plate distances with 
the x-ray tube at different vertical elevations 
with respect to the slit and plate convinced us 
that each of these fluctuations is definitely asso- 
ciated with some local condition permanently 
assignable to different localities in the crystal 
and is not as we thought at first some sort of 
Fresnel interference effect. It seems strange 
however that the extinguished points on the 
two lines forming a pair such as A,B, seemed to 
occur at about the same vertical heights with 
more than accidental agreement. If, as appears 
probable, the two members of each pair are 
merely strongly enhanced reflections from the 
disturbed front and back surfaces of the crystal 
it is hard to see why these disturbances should 
be localized at corresponding points on the two 
sides of the thin flake. No such fluctuation of 
the length of the 


intensity along lines was 


observed with the slit between crystal and 


photographic plate. 


Faint radiation between the exit and entry Laue 
reflections 


We had to settle the question as to whether 
the radiation was absent between the members 
of a pair such as A,B, because of some sort of 
destructive Fresnel interference or whether its 
absence there was simply because the interior of 
the crystal failed to reflect. This we did by intro- 
ducing auxiliary slit jaws J;, Jz with micrometer 
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Fic. 7. Geometry of precision photographic transmission 
spectrometer as used to study entrance and exit surface 
Laue reflections. The slit is now on the opposite side of the 
crystal from the plate. 


screw movements close behind the crystal as 
shown in Figs. 9a and 9b, so as to cut off at will 
either line A; or B,; or both. It was found that 
with either member of the pair cut off the re- 


With both 


members cut off so as to have a slit that would 


maining member was unchanged. 


only permit radiation C,C;’ or C2C»’ reflected 
from the interior of the crystal to reach the 
plate, some fifteen times as much exposure was 
necessary to obtain lines, and these lines were 
hazy, ill defined and unsuitable for precise 
measurement. 

More work should certainly be done in the 
study of these strange effects for themselves but 
for the point at issue the methods so far described 
of studying Laue reflection in macroscopic calcite 
may be criticized (1) as really failing to measure 
the interior grating constant of the crystal which 
is involved in the density determinations, (II) 
as being too highly dependent on peculiarities of 
individual crystals or individual regions in one 
crystal. 


Ill. THE PowbER METHOD A CONCLUSIVE 
ANSWER TO ALL THESE CRITICISMS 


Since one is apparently restricted to investiga- 
tions of a very superficial layer when using 
macroscopic calcite either in Laue or Bragg 
reflection a conclusive answer to all these criticisms 
should come from a study of finely powdered calcite 
of grain size smaller than the extinction layer. If one 
investigates the diffraction angles and density 
of crystals of size 10-* cm on a side, so small that 
extinction is negligible in them, it is obvious that 
a statistical result truly representative of the 
entire sample is obtained both for density and 
diffraction. If it 
balanced forces near the surface of a macroscopic 


be a valid criticism that un- 


crystal render the density there somewhat dif- 
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Fic. 8. Photographic reproduction of Laue reflections from 
exit and entry surfaces of a thin calcite crystal. 


ferent from the mean then one might expect the 
density of the finely powdered calcite to be dif- 
ferent from the mean density of macroscopic 
calcite. Jt is important therefore to measure the 
density of the identical calcite powder used in the 
diffraction experiment. This was done with a 
pycnometer with included thermometer. 


Pycnometric density determination on powdered 
calcite 

The density of the pycnometer water was 
checked by weighing a known volume which had 
been measured out with a standardized ‘‘aus- 
guss”’ pipette. The laboratory distilled water was 
0.08 percent lower density than that given by the 
handbooks and it was necessary to redistill it in 
an all glass still and boil it to remove occluded 
gas. With these precautions the average devia- 
tion from the standard table values of six deter- 
minations of the density of this water was only 
three parts in one hundred thousand. The density 
of both macroscopic and microscopic calcite 
crystals was determined in the same pycnometer 
bottle as a check, the macroscopic samples being 
slivers as large as would go in the neck of the 
bottle. The macroscopic samples weighed 12 
grams while the powdered samples ranged from 
21 grams to 10 grams in amount. The powder 
was made by grinding in an agate mortar and 
rejecting the coarser grains by precipitation in 
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water. The impalpable powder thus obtained had 
some visible grains under an immersion micro- 
scope of a size which we estimated as 2X 10 cm. 
The method of precipitation of course would give 
grain sizes ranging down below this upper limit 
into the submicroscopic. In measuring the dif- 
fraction rings from the powder the few isolated 
specks caused by reflection from the fortuitously 
oriented larger grains were ignored and the 
measurements were taken to the centers of the 
smoothly filled out parts of the rings only. It 
seems likely therefore that the grains actually 
involved in the measurements were very con- 
siderably smaller than these few largest ones seen 
in the microscope and that they constituted the 
bulk of the samples. 

The results of our density determinations on 
macroscopic and powdered calcite are compared 
in Table I with the excellent data obtained by 
Bearden® on macroscopic calcite by weighing in 
fluids. The probable errors of our data are com- 
puted from internal consistency of the observa- 
tions. The reproducibility of the measurements 
of the fine powder is indeed remarkable. In filling 
the pycnometer bottle containing the calcite 
powder with water it was found that pre- 
evacuation of air by means of a vacuum pump 
decreased the tendency for small air bubbles to 
be trapped by the powder. Stirring by rotation 
of the bottle during filling also helped. 

It is easy to show that the solubility of cal- 
cite in water introduces an error of only one part 
in forty thousand in the density determination. 


The powder diffraction measurements 


To obtain the maximum possible precision in 
the absolute determination of the size and shape 
of the lattice a Seeman-Bohlin spectrograph was 
used and reflections were chosen for such primary 
wave-length and Miller indices that the reflected 
beam would return almost directly backward to 
a region FF’ (see Fig. 10) not far from the pinhole 

TABLE I. Density of calcite. 


No. ot 
Orn- DEN- PROBABLI OBSER 
STATI METHOD SERVER SITY ERROR VATIONS 
Powdered Pycnometer | Bollman | 2.71022 +0.00035 6 
Macroscopic | Pycnometer | Bollman | 2.71085 +0.00019 3 
Macroscopic | Immersion Bearden 2.71030 +0.00003 


8 Bearden, Phys. Rev. 38, 1389 (1931). 
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at S through which the primary rays had entered. 
This means that @ in the equation m\=2d sin @ 
is nearly 90°. Considering \ as a known constant 
from grating wave-length determinations if we 
differentiate this equation and replace @ by 
a=42—6 (since the measurements of the diam- 
eters of diffraction rings centered around the 
entrance pinhole are nearly proportional to a) 
we obtain 

a tan a(Aa@ a)=Ad/d (3 


or for small a@ Ad /d = a?(Aa/a). (4) 


Thus if in radian measure a is of order 0.1 the 
relative error in determining the lattice constant 
d is only one one-hundredth as great as the 
relative error committed in measuring on the 
film. 

Nickel K radiation was chosen for the work. 
A rather long wave-length is necessary in order 
to have reasonably intense reflections in calcite 
at angles @ near 90°. On the other hand too long 
a wave-length limits the depth of penetration by 
extinction whereas the whole point of the experi- 
ment is to have the penetration into the small 
particles as great as possible. In fact the two 
requirements, low extinction and strong reflec- 
tion, are more or less mutually contradictory. 
The obvious solution is to tolerate the weakest 
possible diffraction ring intensities that will yield 
acceptable precision in measuring their locations 
in order to get as little extinction as possible in the 
passage of radiation through the particles. The 
reason for choosing nickel for the target was 
because we computed that the Ni Ka lines would 
be just reflected in the third order from the (110) 
planes (cleavage system) of Fig. 6 at an angle @ 
slightly less than 90°. Measurements eventually 
were discarded on this reflection as it was found 
to be too faint and diffuse on account of the 
enormous dispersion at this large angle but 
three planes with other Miller indices were found 
to yield all the information needed from the 
nickel radiation. The target of the tube was made 
by nickel plating an old water-cooled copper 
target. A cylindrical Nonex glass envelope about 
1} inches in diameter provided with a thin 
concave glass window in the end of a short side 
tube was blown around the cathode and target 
and a thoroughly outgassed charcoal trap in a 
side tube was kept in liquid air during operation. 
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The spectrograph had a lead alloy lining to 
reduce the general scattering. Diffraction ex- 
posures taken without the calcite powder but 
with the x-rays simply incident on the lead 
showed strong diffraction rings inconveniently 
located close to the calculated positions of dif- 
fraction rings from calcite. These were com- 
pletely eliminated by covering the lead in the 
radiated area with several sheets of gold foil 
which in a test exposure showed no rings what- 
ever in the region where our calcite rings 
appeared. The spectrograph had a tight cover 
provided with a gasket and could be evacuated 
with a fore pump. This precaution we found 
necessary to reduce the fogging and general 
background of the film. Even so the calcite 
diffraction rings were very faint and are unfor- 
tunately unsuitable for reproduction by half-tone 
process. 

Fig. 11 shows the x-ray tube and the spectro- 
graph with cover removed. The calcite powder 
which was made to adhere to the cylindrical 
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Fic. 9a. Auxiliary slit jaws used in study of entrance and 
exit surface Laue reflections. 9b. Photograph of transmis- 
sion spectrometer slit and crystal mounting with auxiliary 
jaws. 
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Fic. 10. Geometry of Seeman-Bohlin spectrograph. 


surface with a tiny quantity of collodion can be 
seen as a white patch on the inner wall of the 
spectrograph. Part of the film holder and the 
slit both on the opposite side of the spectrograph 
can also be seen. The calcite powder sample had 
a thickness of 0.25 mm to correspond with the 
thickness of the double-coated x-ray film. The 
diameter of the spectrograph D=6.905 inches 
was carefully calipered inside at a number of 
points. The greatest deviation from uniformity of 
D was one part in 7000. 


Identification of powder diffraction rings with 
Miller indices 


It is of course absolutely essential to identify 
beyond all doubt the Miller indices corresponding 
to all powder rings measured. If a reciprocal 
lattice for calcite is constructed the coordinates 
of each point of which correspond to the Miller 
indices of a set of planes and a sphere is struck out 
centered at the origin with radius R= (2 sin 8)/d 


a(1+2 cos? y—3 cos* y)! 


dnihehs = 


*For the structure of calcite and its representation in different systems of reference see Ewald-Hermann, 
Structurbericht, p. 292. 
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[(hy+hee+hs?) sin? y+2(Aihe+hihs+hehs)(cos* y—cos +) } 
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Fic. 11. Nickel target x-ray tube, Seeman-Bohlin spectro- 
graph and pycnometer bottle. 


Ewald has shown that any points of the reciprocal 
lattice lying on this sphere will correspond to the 
Miller indices of planes which can reflect radi- 
ation \ under angle @. Two such reciprocal lat- 
tices were graphically constructed. The points 
lie in planes equally spaced in a vertical direction. 
In each plane the points lie at the corners of 
lozenge shaped cells having angles of 60° and 
120°. Fig. 12 shows a projection of the reciprocal 
lattices looking along the vertical axis. A sphere 
of radius R= (2 sin @)/X was struck out in one of 
these for @=75° 24’ and in the other for 
6=83° 39’. This sphere intersects the equally 
spaced horizontal planes in circles as shown in 
the figure. Thus if a circle passes through a 
lattice point in the projection that point may 
correspond to a possible reflection. The method 
of constructing the reciprocal lattice by starting 
for convenience with the calcite lattice in the 
hexagonal system gives all possible calcite Miller 
indices and some which are not a part of the 
calcite system. These can be easily discarded 
since they do not yield integers for the Miller 
indices in the rhombohedral system. This method 
is thus a completely exhaustive one for enumerating 
the Miller indices of all calcite planes which can 
be responsible for reflections of \ at the angle @. 
The following formula’ gives the grating space 
dy,n,r, for any set of Miller indices expressed in 
the rhombohedral crystallographic (not cleavage) 


system 


’ 
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TABLE I]. Designation of all possible planes that have values 
of d appropriate to the regions in which powder rings were 


observed. 
DeprTH oF 
Srrvc- PENETRA- 
Hex Roms. | Gratinc | Weicrt TURE TION TO } 
Inpices, | INDICES, SPACE Factor, | Factor, INTEN- INTENSITY 
HiheHe¢ hihohs) dayhohs P F SITY cm 
476 0.8532 12 3.66 161 
374 8528 12 4.38 | 230 
163 8557 12 8.04 774 0.585 X 10 
233 8458 6 11.3 762 0.416 X10 
141 8586 6 4.38 115 
932 8619 é 1.83 20 
677 8349 6 11.3 762 0.416X 10 
357 8389 12 0 0 
733 8326 6 0 0 
146 8338 12 3.66 161 
| 005 &366 6 0 0 
| 330 8304 6 7.58 343 | 0.62 1073 





in which to an accuracy sufficient for the identi- 
fication of diffraction rings the constants can be 
assumed to have the value a=6.361A, y =46° 7’. 
A knowledge of the relative intensities to be 
expected from the various planes served as a 
final identification test. These were computed 

from the formula 
I=KPF°, (6) 


in which K is simply a proportionality constant, 
P is the permutation factor entering because 
some planes admit of reflection with the crystal- 
line particle in more orientations than do other 
planes, and F is the crystal structure factor. In 
computing the crystal structure factor, values of 
the atomic structure factors for Ca, C, and O 
were obtained from the work of Pauling and 
Sherman.’ Table II gives all of the possible 
planes expressed both in hexagonal and rhom- 
bohedral crystallographic indices that have 
values of d appropriate to the regions in which 
powder rings were observed. Referring to the 
rhombohedral indices the computed intensities 
immediately identify (163) and (233) as the 
planes giving the two observed reflections in the 
region 75° 24’ while in a similar way (677) was 
identified with the reflection in the 83° 39’ 
region. As before stated the reflection (330) 
which we have referred to in Fig. 6 did not 
present sufficient intensity and contrast to be 
used with certainty. 


Depths of penetration into the crystal particles 


The depths of penetration to } intensity in 
Table II were computed from the formula for 





‘© Pauling and Sherman, Zeits. f. Krist. 81, 1 (1932). 
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the extinction coefficient yu, 
be=4rF5/Z X sin Oo, (7) 


in which Z is the number of electrons in the 
calcite unit cell, 6=10~° is the difference between 
the refractive index and unity for the wave- 
length A=1.658A, 6=90° approximately. Note 
that the largest crystals observed with the 
microscope in our powder sample have dimen- 
sions only half the shallowest half penetration 
depth computed above. Thus in the worst cases 
the beam only lost about a quarter of its original 
intensity in passing through a single crystalline 
particle. 
Measurements of diffraction angles 

Five different exposures were taken but since 
not all of these showed the same degree of con- 
trast some of the lines could not be measured 
with certainty on certain of the exposures. 
Table III giving the observed distances 2s 
(Fig. 10) measured on the film indicates the 
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Fic. 12. Projection of calcite reciprocal lattice and 
circular intersections of propagation sphere with the lattice 
planes for two radii. 
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TABLE III. Observed values of the distance 2s measured on 
the film. 


Fitm No 

MILLER SPECTRAL 1 2 3 4 
INDEX LINE DISTANCES 2s (INCHES MEAN 
hihohs a) 7.05 7.07 7.06 7.05 7.05 7.056 
(163 ae 6.79 6.75 6.75 6.80 6.75 6.768 
3 3 3) a 5.70 5.67 5.685 
Xs as 5.39 5.30 5.345 

79 a 3.78 3.79 3.80 3.75 3.76 3.776 
tors ae 3.30 3.30 3.26 3.29 3.30 3.290 


inner consistency here. The resulting probable 
relative error in d is of course far smaller. 
Now by means of the formula 


d=X/(2 cos 2s/4D), (8) 


which can be easily shown to follow from the 
geometry of the spectrograph and assuming the 
wave-length values for Ni Ka, and Kaz as 
measured by Erikson" based on Siegbahn’s con- 
ventional value of the calcite lattice constant 


dy = 3.02904A at 18°C for calcite, 


ha, = 1.65450, Aa, = 1.65835, 


we can compute the lattice constants of our 
planes in terms of the Siegbahn conventional 
scale. These are given in Table IV. 
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TABLE IV. Lattice constants of planes in terms of Siegbahn 
conventional scale. 


Indices 
Cleavage 


Indices 
Crystallograph 


Angstroms 


Rhomb. d d Cell 
hihehs From @, From a» daverage  CilaC3 
163 0.85500 0.85472 0.85486 6,14,2 
233 .84510 .84494 .84502 12,88 
677 .83504 .83510 .83507 48,8 


Cohen” has analyzed the sources of error in 
precision work with the Seeman-Bohlin spec- 
trograph. Using his formula and also taking into 
account the inner consistency of the measure- 
ments on 2s we estimate very conservatively that 
the total error from both systematic and random 
causes cannot exceed 0.05 of one percent in the 
above results of Table IV. The agreement 
between the results on a; and ag is in accord with 
this estimate. 


Determinations of the angle 8 of the calcite 
cleavage cell 

The above data are all that we need to deter- 
mine the angle 8 characteristic of the calcite 
cleavage cell. This is the angle which enters into 
the function (8) in the formula (2) by which e 
is computed. If cyce and c; are the Miller indices 
in terms of the cleavage cell we have 


c(1+2 cos* B—3 cos? B)! 


dcie2e3 = 


(9) 


[ (c1?7+¢2? +c?) sin? B+ 2(€:Co+6263+0€1C3)(cos* B—cos B ) } 


The ratio of the d’s for any two different sets of planes is then sufficient to calculate 8. If the 
subscripts 7 and & refer for brevity to two such different sets of planes we have 


cos B= 


[ne — Ndi} = | Ene —n bu? P+4[a bi? —ay ]ime—mibi2 }}! 


ae (10) 


2[a,—a di? | 


where );,=d;/dy, m=c?>+c?+<¢;?, 
The three sets of planes observed in our powder 
camera thus permitted three ratios };, to be 
formed from which three independent deter- 
minations of 8 were computed as follows 


Pairs of 
Cleavage cell indices Values of 8 Mean value of 8 
(4,8,8)/(12,8,8) 101° 53’ 40” 
(4,8,8)/(6,14,2) 101° 52’ 51” 101° 54’ 18.448” 


(Probable error from 
inner consistency) 


ut Eriksson, Zeits. f. Physik 48, 360 (1928). 


(12,8,8)/(6,14,2) 101° 56’ 24” 


N=2(C\Co+CiC3+C2C3), Ap=N~—My, 


nN; —M;=a;j. 


This should be compared with Bearden’s 
value of 8 for macroscopic calcite” 
B=101° 54’ 4"+10”. 
Determination of d for the cleavage planes in 
microscopic calcite 
We are now in a position to compute the value 
of d for the cleavage planes in our microscopic 


2 M. U. Cohen, Rev. Sci. Inst. 6, 68 (1935). 
18 Bearden, Phys. Rev. 38, 2089 (1931). 
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calcite for comparison with macroscopic values. We proceed by first computing the lattice para- 


meter C of the cleavage cell by means of the formula 


C;=d; 


Three different values distinguished by subscript 
i were obtained from the data for the three dif- 
ferent sets of planes observed with the result 


C.=12.8238A, C3=12.8194A, C,=12.8159A, 


all of course on the Siegbahn scale. The mean 
value is 


C=12.8197+0.0016A (Siegbahn scale) 
B=101° 54’ 18’°+.48”. 


Remembering that the cleavage dq is really d soo 
in the 32 molecule cell we have 


12.8197(1+2 cos* B—3 cos? B)! 
dioo= 
(16 sin? g)? 
12) 


= 3.02823A on the Siegbahn scale. 


The corresponding value for macroscopic calcite 
on the same scale is dso=3.02904A. 


Value of electronic charge from grating wave- 
length data 


Transformation from the Siegbahn scale to the 
absolute scale of lengths as determined by ruled 
gratings was effected by a comparison of the 
crystal results of Wennerléf™ on Cu Ka, radi- 
ation Ae; = 1.537395 (Siegbahn scale) with Bear- 
den’s measurement of the same wave-length with 
a ruled grating \, = 1.5406 X 10-* cm. Using then 
the value of d on the absolute scale together with 
M =100.078+0.005, Q=2.89270+0.00021 x 10" 
abs. e.s.u., p=2.71022+0.00035 g/cc, (8) 
= 1.09602 + 0.00048 one obtains by formula 2 


e=4.799+0.007. 


IV. CONCLUSIONS 


The results of this investigation we believe 
conclusively and permanently answer all three of 
the objections raised to the x-ray crystal method 


“ Wennerlof, Ark. Met. Astr. O. Fys. (A) 22, No. 8 
(1930). 


(m; sin? B+n;(cos? B—cos 8))? ¢=2, 3,4 


(1+2 cos* 8B—3 cos? B)! B=101° 54’ 18” 


of determining e. The angles of diffraction are 
representative of the entire powder sample since 
the half depth of penetration computed from the 
extinction coefficient is double the dimensions 
of the largest crystals used in the sample. The 
density determination therefore applies to iden- 
tically the same calcite as has been used to 
diffract the x-rays. 

The agreement between the densities of micro- 
scopic calcite and macroscopic calcite as well as 
the agreement of absolute size and shape of the 
lattice cell for these two is surprising evidence of 
the permanence of the lattice structure against 
mechanical maltreatment by powdering. 

The elementary diffraction theory of Laue or 
Bragg must apply to the case of diffraction by 
crystals so small that the radiation is negligibly 
extinguished in passing through them and hence 
any critical objection to the crystal determina- 
tion of e based on doubts as to the absolute 
rigor of the dynamical theory of x-ray reflection 
are invalidated by our results. 
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Deep Terms in the Spectra of Pb V and Bi VI 


G. K. SCHOEPFLE, Department of Physics, A. and M. College of Texas 


(Received March 7, 1936) 


In earlier investigations the Pt I-like isoelectronic sequence was extended to Pb V and Bi VI, 
establishing terms due to the 5d%6s, 5d%6p, and, later, 5d°7s and 5d’ configurations. This paper 
extends further the previous results by Mack and by the author to include terms assigned as aris- 
ing from 5d%8s, 5499s, 5d°7p and 5d°8p configurations. A total of 69 term values involving 205 
lines are listed for Pb V, and 40 terms involving 94 lines for Bi VI. The wave-length measure- 
ments for the lines involved are those made by Arvidsson and range from 1400A to 200A, and 
in more than half the cases the discrepancies are not more than +1 cm“. The absolute values 
of the 5d°; )26p3/2 J=3 terms are computed by a Ritz formula. 


NDEPENDENT extensions of the Pt I-like measurements used being those made on a 

isoelectronic sequence to PbV and BiVI_ grazing incidence spectrograph at Upsala by 
were made and announced almost simultaneously Arvidsson.° 
by Mack! and by the author.’ As the results In Table I the 69 term values assigned to 
were not in complete agreement, a more detailed Pb V are given, which involve a total of 205 lines 
study has been made by Mack’ and also by the for over half of which the discrepancy is not 
author,‘ including a study of the spectra of more than +1 cm. The wave-lengths of newly 
Pb IV and Bi V5 The results are now almost in @Ssigned lines are given in Table II. Similarly the 
agreement, differences being discussed in the 40 terms assigned to Bi VI, which involve a total 


comments to Tables I and III. The wave-length of 94 lines for over three-quarters of which the 


£ ; ee, ‘ P ss Pe 
region concerned is from 1400A to 200A, the discre Pan OSs mere than +1 cm“, are given 
in Table III. The wave-lengths of newly assigned 


lines are given in Table IV. The terms assigned 
as arising from 5d", 5d%6s, 5d°7s, and 5d%6p con- 
figurations are as assigned by Mack. By a Ritz 


1A. T. Goble and J. E. Mack, Phys. Rev. 42, 909 
(1932). 

?G. K. Schoepfle, Phys. Rev. 43, 374 (1933). 

3]. E. Mack and M. Fromer, Phys. Rev. 48, 357 
(1935). 

*G. K. Schoepfle, Phys. Rev. 45, 747 (1934). 

5G. K. Schoepfle, Phys. Rev. 47, 232 (1935). 


formula, the absolute value of the 5d°;,26p3/2J = 3 


® Gustaf Arvidsson, Ann. d. Physik 12, 1 (1932). 
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z indicates the classification according to the state of ionization by Arvidsson 
‘able I. b = broad, d =diffuse, a = found in arc. 
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284.329 
58.636 
387.870 
389.745 
390.882 
395.280 
395.770 
403.217 
404.644 
409.728 
410.705 
416.950 
418.264 
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TABLE II. Newly classified lines of Pb V. 
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strengthened by the expected transition 5d°;,.6p 


SPECTRA 
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AND 
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TABLE III. Energy levels and transitions in Bi VI. 
M indicates lines listed by Mack; a number locates the line in Table IV. Mack's choice of 298862, J 


is 


»— 5d°3,2.7s designated as 26, apparently overlooked by him 
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TABLE IV. Newly classified lines of Bi V1. 
z indicates the classification according to the state of ionization by Arvidsson.* The number in the column 
labeled “Designation” refers to Table III. b=broad, d=diffuse, a =found in arc 
INTEN DESIGNA- DISCREF INTEN DESIGNA DIscREP 
2 SITY A(va TION ANCY SITY A(va v(va TION ANCY 
2 189.274) 528334 46 0 VI 4 710.32 140782 12 +2 
VI 2 214.318 | 466596 45 +4 VI 5 713.08 | 140237 6 0 
1 230.710 | 433445 44 —2 Vi 2 723.56 | 138206 38 —1 
VII 2 276.397 | 361798 34 0 \ 10 738.17 135470 21 0 
VI 1 277.591 | 360242 40 0 , » . e 24 = 
VI 278.095 | 359589 39 0 VI | 10b | 743.13 | 134566 MI +3 
VI l 292.080 | 342372 30 +] , ace , 20 0 
VI 3 295.530 | 338375 41 0 VI ‘ 1OB.44 | 150134 M2 0 
VII 0 300.225 | 333084 31 +3 1 773.76 | 129239 23 +1 
VI 4 304.278 | 328647 35 +6 V 3 775.43 | 128961 13 —3 
VI 6 306.718 | 326032 39 0 V 1 777.05 | 128692 9 —5 
VI 5 309.976 | 322606 22 0 IV 2 790.425 | 126514 14 —2 
VI 2 314.982 | 317478 25 +1 | 809.82 123484 15 +] 
VI 3 326.695 | 306096 32 —2 V 3 850.66 | 117556 18 0 
VI 3 337.590 | 296217 33 —1 IV 4 930.65 | 107452 8 —1 
1 345.479 | 289453 27 +5 1 993.10 | 100695 2 —1 
IV 8 376.660 | 265491 28 0 VI 7b 1004.75 99527.2 16 —6.2 
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\ + 619.13 | 161517 4 0 IV 2 1131.39 88386.9 1 +9 
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\ 4 693.29 | 144240 7 —2 V 0 1394.46 71712.3 11 +.3 
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term is computed as 562,600 cm~ for Pb V and 
as 761,800 cm~ for Bi VI. The ionization poten- 
tials are 69.40 volts and 93.97 volts respectively, 
which values are 1.3 percent and 6.6 percent 
higher than those obtained by the extrapolation 
method by Mack. In Fig. 1 a Moseley diagram 


is given. 
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The Pressure Effect on the C Band of Strontium Hydride 


R. F. HUMPHREYS AND W. R. FREDRICKSON, Syracuse University, Syracuse, N. Y. 


(Received June 6, 1936) 


At low hydrogen pressures the emission spectrum of the 3808A band of SrH shows a sharp 
break-off at R(18), while at high pressures the lines extend to R(49). The presence of lines with 
high K values allows an accurate determination of the constants of the upper state. Analysis 
shows the presence of an uncoupling term —0.1415 (K—1) in the upper state. Perturbations 
due to resonance and predissociation are present. The electronic configurations and probable 


dissociation energies are considered. 


F the known five band systems which the 

spectrum of SrH exhibits the C system at 
3808A is particularly interesting, because its 
main band, the (0,0), has a line structure that is 
intense and well defined to a certain point and 
then abruptly breaks off, leaving no further trace 
of lines to be found. An analysis of this band has 
shown it to be a C°S—>N* transition.' The only 
other band observed in the system is a weak 
(0,1) band, which permitted a sufficiently com- 
plete analysis to show the same kind of abrupt 
breaking-off in rotational structure. This phe- 
nomenon has been attributed! to predissociation 
of the upper (“‘C’’) state due to the influence of a 
neighboring electronic state (probably the D 
state). 

An analogous situation had long been known 
to exist in the C system of CaH.? Later Grund- 
strém and Hulthén’ obtained its spectrum with 
the pressure of the hydrogen surrounding the 
calcium greatly increased. The bands were no 
longer predissociated, but were more or less 
normal in appearance. Bearing in mind the close 
analogy between the spectrum of CaH and SrH 
it was thought worth while to investigate the 
effect of pressure on the SrH C band. 


1W. R. Fredrickson, M. E. Hogan, Jr., and W. W. 
Watson, Phys. Rev. 48, 602 (1935). 

? Mulliken, Phys. Rev. 25, 509 (1925). 

3’Grundstrém and Hulthén, Nature 125, 634 (1930). 





EXPERIMENTAL PROCEDURE 

The spectrum was obtained from an arc be- 
tween a copper cathode and a copper cup con- 
taining metallic strontium as an anode. The arc 
burned in an atmosphere of hydrogen at pres- 
sures varying from a few mm of Hg to 700 mm of 
Hg. The cathode was water cooled, and a con- 
tinuous stream of water was directed on the glass 
walls of the tube containing the arc. With this 
arrangement it was possible to operate the arc 
continuously for two or three hours. 

A 120-cm concave grating with a dispersion of 
5A/mm in the second order was used in a Row- 
land mounting. Exposure times were of the order 
of four hours for second-order plates. 


ANALYSIS 


When the C band is photographed with a 
hydrogen pressure of a few mm of Hg its rota- 
tional structure breaks off at K‘=19. When the 
pressure is increased to 400 mm or more the band 
no longer exhibits this predissociation but dis- 
plays a fairly normal structure. In addition to the 
(0,0) very weak traces of the (1,1) and (2,2) 
bands appear. The assignment of the K values is 
given in Table I. Combination differences permit 
these assignments up to K = 28, beyond which the 
P branch becomes too weak to observe. For 
K> 28 the R lines were fitted to a Fortrat. para- 
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IN BAND SPECTRA 


TABLE I. Assignment of frequencies in the C*=—>N°Z (0,0) band of SrH (cm units). Two lines blocked together indicate the 
boundary of a wide line. 


Intensity v (obs.) P(K) R(K) 
2 26,245.01 9, 13 
9 Bad 8, 9, 10, 11, 12, 
254.45. 13, 14, 15, 16 
9 bad 7,9 
259.62 11, 17 
8 262.45 6, 18 
3 268.18 5, 19 
7 271.98 4, 20 
3 278.13 3 
3 279.85 21 
3 284.62 2, 22 
2 291.26 1, 23 
2 298.24 24 
1 306.85 25 0 
2 314.79 26 1 
4 323.58 27 2 
2 333.36 28 3 
3 343.08 4 
3 353.21 5 
1 357.66* 
1 362.24 6 
0 368.16 7 
1 374.48 7 
378.73 . 
4 fod ’ 
4 391.82 8 
4 403.95 9 
2 407.65 9 
0 413.23* 
5 417.77 10 
0 425.17 11 
6 432.01 11 
9 445.93 12 
0 454.19" 
5 461.26 13 
3 474.27 14 


* Unassigned lines. 


bola and thus given their K values. The resolving 
power of the grating was not sufficient to dis- 
tinguish the spin doubling nor to resolve the ab- 
normally broad lines in the head of the P branch 
and in several other places throughout the R 
branch. 

A study of the photographic plate and of the 
assignments given in Table I reveals the presence 
of several perturbations in the band, evidenced 
by anomalous displacements, irregular intensity 
distribution and abnormal doubling of certain 
lines. In order to obtain information about these 
perturbations an attempt was made to write a 
theoretical equation for the lines of the R branch. 
A deviation of the observed position of a line 
from its theoretical would thus indicate a per- 
turbation. To write such an equation necessitates 
the determination of the coefficients of the energy 


level equation 


T(K)=To+BoK(K+1)+DoKX(K+1)? (1) 


Intensity v (obs P K R K) 
494.01 ‘ 

‘ fer] 15 
4 510.73 16 
531.13 ‘ 

Berd 17 
4 547.38 18 
2 569.16 19 
2 586.09 20 
1 605.97 21 
1 626.09 22 
1 647.17 23 
1 668.15 24 
1 689.44 25 
0 711.54 26 
1 734.40 27 
2 757.22 28 
1 781.66 29 
0 800.19 30 
0 816.72 31 
0 846.97 32 
0 868.61 33 
0 888.85 34 
0 909.68 35 
0 934.56 36 
arc line 963.32 37 
arc line 990.18 38 
0 27,017.39 39 
0 030.17 40 
0 049.48 41 
0 076.44 42 
0 105.00 43 
0 140.19 44 
0 165.16 45 
0 201.09 46 
0 231.69 47 
0 262.44 48 
0 293.47 49 


for the two states in question, and thus determin- 
ing the coefficients of the line equation. For the 
lower state By” and D,’’ have been calculated‘ 
(Table II). The levels of the upper state were 
found with the aid of combination differences and 
were observed to be somewhat irregular. In view 
of this, the procedure used by Grundstrém® for 
CaH bands was followed. An additional function 


f(K) was assumed to exist, so that 


T’'(K)=T)'+ Bo K(K+1)+D)'K*(K+1)*+f(K). 
TABLE II. Constants for the N and C states of SrH. 


Normal state C state 


By’ =3.9416 
D,' = — 1.3986 X10~* 
we = 1206 (experimental) w, = 1349 (experimetal) 
Xe = 17 (experimental) X <2, = 24 (experimental) 
ro =2.067A 
a(uncoupling) = —0.1415 


By’ =3.6344 
Do'' = —1.287 X10-* 


* Watson, Fredrickson and Hogan, Phys. Rev. 49, 150 
(1936). 
5 Grundstrém, Zeits. f. Physik 95, 574 (1935). 
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Fic. 1. A, 7’(K)/K vs. K*®. The circles are obtained from 
upper state levels; the crosses represent additional levels 
giving rise to the anomalous doubling. The y intercept of 
the straight line equals 2By’ and the slope equals 4D)’. 


From this we may obtain 


A,T’(K)/K=(T'(K)—-T’(K—-1)]/K=2B,’ 
+4Do'K?+(f(K) —f(K —1)]/K. 


If Ai\7’(K)/K is plotted against K? 
expect to get a straight line only when [/f(K) 
—f(K —1) ]/K is zero or negligible. Such a plot is 
The straight line occurs from 


we should 


given in Fig. 
K*=400 on. Therefore, by substituting the term 
values for levels from K=20 to K =30 into Eq. 
(1), the values of By’ and Dy’ were obtained 
(Table II). For points at small K? values in Fig. 1 
there seems to be a general sinking below the 
straight line (except, of course, for points marked 
anomalous doublets 
This sinking of 


with x’s, which 


and thus cannot be considered). 


represent 


the points suggests the possibility of an addi- 
tional term, a/K, where a@ is a constant. Thus 
[f(K)—-f(K-1)]/K=a/K or f(K)=a(K+n) 
+const. where » is any integer. The assumptions 
are made that n=1, the sign is negative and the 
const. =0. These assumptions are quite arbitrary 
and are chosen to produce the best agreement 
between the theoretical 
numbers of the levels. The upper state, therefore, 


observed and wave 


may be represented by 
7'(K) =7T)'+ By’ K(K +1) 
+ Dy'K*?(K+1)*+a(K - 
If we combine the values of B and D for the two 
states we get the following equation: 
R(K) = 26,298.53 +-7.5760(K +1) 
3072( 


+0 K+1)?—2.6856 x 10-4(K +1)? 
—11.1610-*(K+1)4+eaK. 


By substitution in several presumably 
turbed lines the average a is found (Table II). It 
of “uncoupling” of 


unper- 


is a measure the orbital 


AND 
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TABLE III. Probable configurations and dissociation products. 


Probable Dissociation Dissociation 
State Configuration Products Energy 
N-*> 5so0*4doa Sr, 'S+H, 28 12,000 cm 
A—?Il 5so*5 pr Sr, *P+H, 2S 13,300 
B-—*> 5sa*5 po Sr, °P+H, 2S 12,300 
D—-*z Ssa5 po" Sr, *P+H, 2S 10,000 (?) 
F-—"Il  SsoSpoSpr —- Sr, 3P+H,, 2S 7,800 
C—*> 5sc°6se Sr, 3S +H, 2S 15,000 
angular momentum in the C state. The wave 


numbers of the line of the R branch were calcu- 


lated from the above equation and then sub- 


tracted from those of the observed lines. Per- 
turbations were thus found to occur in four 
regions: R(7) to R(11), R(17) to R(18), R(30) to 


R(40) and R(40) to R(49). For high K values the 
differences between calculated and observed posi- 
tions become large due to a lack of powers above 
the fourth in the theoretical equation, but these 
differences are readily distinguishable from per- 
turbations by their uniformity in trend. 

In the region K=7 K=11, the perturba- 
tions are of the resonance type, but those in the 
regions following K =17, 30 and 40 must be due 
to the overlapping of levels which are above the 
dissociation limit of their own state and conse- 
quently would cause predissociation of the C 
state if the hydrogen atmosphere were not under 
pressure. Unfortunately, the D state has not 
been sufficiently analyzed to permit a quantita- 
tive investigation of these perturbations. 


DIssOcIATION PRODUCTS 

In a recent paper® Grundstrém has called at- 

tention to the analogous spectra of the hydrides of 

the alkaline-earths. If one makes the assumption 
logical by comparison with MgH and CaH 

that the C state of SrH predissociates into the 


’P Sr atom (the dissociation product of the dis- 


turbing D state) and a normal 7S H atom, then 
the relative positions of the levels of the SrH 
molecule and the Sr atom can be obtained. 


Again, by analogy, one would assume the product 
of dissociation of the 8S Sr atom, 
which would give a value of the energy of dissoci- 
ation of approximately 15,000 cm~'. To check 
this the ordinary extrapolation methods of Birge 
To obtain w, and x,.w, 


C state to be a 


and Rydberg’ were used. 


6’ Grundstrém, Zeits. f. Physik 99, 595 (1936). 
‘ Jevons, Report on Band Spectra of Diatomic 
p. 193. 
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for the normal state the C (0,0), (0,1) and 
D (v’,2), (v’,3), (v4) bands* were used. With 
these and the C (1,1), (2,2) bands, with heads at 
26,380 cm and 26,495 cm™, respectively, the 
vibration constants of the C state were calculated 
(Table II). These gave for the dissociation value: 
by Rydberg’s niethod, D=10,000 cm; linear 
extrapolation, D= 18,200 cm™. Since neither of 
these agrees with the position of an atomic state 
of Sr, one concludes that the *S, occurring at 
about the average, is more likely than any other 


§ See reference 4. The (v’,4) of the D system, located at 
16,272 cm™, has not yet been reported. 
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to be the state of excitation of the atom after 
dissociating from the C state. This would indicate 
an electron configuration® of the molecule in the 
C state of 5sc*6sc. The dissociation values for all 
the known states of SrH, obtained from products 
of dissociation, are listed in Table IIT.'° 

We wish to thank Dr. A. J. King for generous 
supplies of pure metallic strontium. 


’ This configuration differs from that assigned by Grund 
strém (reference 3), but the difference might possibly be 
due to a typographical error. 

10 These values differ somewhat from those listed by 
Grundstrém (reference 3) though the methods used were 
identical. 
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Mass-Spectrograph Analysis of Beryllium 


WALKER BLEAKNEY, JOHN P. BLEWett, RuBBY SHERR AND ROMAN SMOLUCHOWSKI, Palmer Physical Labor 


atory, Princeton 


University, Princeton, N. J. 


(Received July 10, 1936) 


A new source of ions suitable for positive ray analysis was developed in the form of a molecu- 


lar beam from a hot oven. Ions were produced in the beam by electron impact and focused on 


the entrance slit of the mass spectrograph by suitable electric fields. The method was applied 
to the study of beryllium but only Be® was found. The experiment placed an upper limit on the 


existence of Be* of one part in ten thousand. 


LTHOUGH most of the lighter elements 

have been satisfactorily investigated, anal- 
yses of the isotopic constitution of beryllium have 
not yielded consistent information. Watson and 
Parker! have reported the existence of Be® in an 
intensity ratio Be® : Be*= 2000 : 1 from a study 
of the band spectra of BeH. On the other hand, 
Olsson? failed to observe any Be*H lines on 
heavily exposed photographs of the bands. 
Aston’s* analysis by means of accelerated anode 
rays, using BeF, has given no conclusive results 
because of the difficulty of obtaining sufficient 
intensity. Furthermore, while Be* appears pre- 
sumably as a disintegraion product of at least 
five nuclear reactions,‘ two, or possibly three, of 
these would be satisfied by the existence of a 
stable Be® nucleus, although an alternative ex- 
planation has recently been suggested. From 
these considerations it would seem of some im- 


1 Watson and Parker, Phys. Rev. 37, 167 (1931). 
? Olsson, Zeits. f. Physik 73, 732 (1932). 

* Aston, Mass Spectra and Isotopes. 

* Wells and Hill, Phys. Rev. 49, 858 (1936). 


portance to be able to set an upper limit for the 
relative abundance of Be’ in natural sources of 
beryllium. Our investigation was undertaken 
with this end in view. 

The apparatus that 
Blewett.® The glass tube carrying the ionization 
chamber and electron filament was modified in 
order to permit the introduction of a small 
furnace in a manner shown in Fig. 1A. The 
furnace O can be inserted through either of the 


used is described by 


ground glass joints as indicated in the figure. It 
consists of a cylindrical spiral 7, Fig. 1B, made 
from a strip of 1- or 2-mil tantalum, 2 cm long, 
8 mm wide spot welded to two nickel rods which 
are supported by tungsten leads. The lip extend- 
ing over the side serves to partially collimate the 
vapor stream. It has been found desirable to 
use a radiation shield N, usually of sheet ni- 
chrome, to protect the glass walls and to reduce 
sputtering as much as possible. 


A 1-kw 120/20-volt transformer supplies the 


. Blewett, Phys. Rev. 49, 900 (1936). 
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Fic. 1. A. Electrode arrangement of the ion source. B. Enlarged cross section of the furnace. 


current (between 15 and 25 amperes) which 
passes directly through the tantalum strip. The 
metal under investigation is placed in contact 
with the furnace and vaporizes when it reaches 
the proper temperature. Ions are formed in the 
molecular beam entering the ionization chamber 
C, Fig. 1A, by electrons accelerated from the 
tungsten filament F. The ions are then acceler- 
ated into the analyzing chamber through a slit 
system S in the usual fashion. 

To test the practicability of this arrangement, 
the furnace was packed with copper filings and 
heated to about 1000°C. Peaks of the order of 
several thousand divisions on the detecting in- 
strument were easily obtained. It may be noted 
in passing that the presence of argon in the tube 
augmented the size of the Cu* peaks to an 
astounding degree, an effect probably caused by 
collisions of the second kind.* This supposition 
was substantiated by the disappearance of the 
effect below the ionization potential of argon. 

Beryllium proved slightly more troublesome 


6 Duffendack and Black, Phys. Rev. 34, 35 (1929). 





because of the degassing of the glass walls. A 
small electric fan, however, served to keep the 
walls as well as the waxed joints quite cool. Peaks 
as high as 13,500 divisions on our scale (2.7 107"! 
amp.) were obtained and a search made for Be’*. 
The ionization efficiency was a maximum be- 
tween 40 and 45 volts, so that in this range the 
presence of O,,**, whose appearance potertial is 
about 50 volts, was not possible. The background 
was negligible and the peaks sufficiently sepa- 
rated so that the tail of the mass 9 peak disap- 
peared completely before the mass 8 position 
was reached. Repeated runs failed to reveal the 
presence of any ions at the mass 8 position. We 
believe this experiment places an upper limit on 
the existence of Be® of 1 part in 10,000 of Be’. 
A small peak was observed at the mass 10 posi- 
tion, but since its size bore no definite relation 
to that of the mass 9, it was attributed to the 
presence of a small amount of BeH*. 

The success of the present method leads us to 
hope that it will now be possible to investigate 
heavier metals which have heretofore been difh- 
cult to handle by other means. 
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The Optical Properties of Nonpolar Liquids 


Hans MvuELLER, The George Eastman Laboratory, Massachusetts Institute of Technology, Cambridge, Mass. 


(Received April 27, 1936) 


A theory of the refraction, Kerr effect and light scattering 
for liquids with nonpolar axially symmetric molecules is 
developed. The molecules are assumed to possess a limited 
freedom of rotation. This assumption alone does not alter 
the classical equations. It is necessary to consider the 
quasi-crystalline grouping of neighboring molecules which 
produces an anisotropic polarization field. The anisotropy 
of the Lorentz forces is calculated by combining the clas- 
sical cavity method with Ewald’s lattice theory. It is 
shown that the procedure of Raman and Krishnan is not 
justified. With a low potential barrier the quasi-crystalline 
grouping gives an increase of the refraction, Kerr effect, 
depolarization and intensity of the Rayleigh scattered 
light. This is true also for an arbitrarily high potential 
barrier, provided the anisotropy of the Lorentz force is 
sufficiently large. These results agree with the observed 
increase of the refraction in compressed gases and in 
some solids, and with the observed increase of the de 
polarization in compressed gases and in liquids near the 
critical temperature. The observations in liquids are ex- 


T is a well-established fact that the classical 
theories of refraction, light scattering and the 
electro-optical Kerr effect fail when applied to 
liquids. In liquids the molecules appear to have a 
slightly decreased polarizability and a much 
smaller optical anisotropy than in the gas. Since 
for most substances the absorption-spectrum is 
essentially the same in the gaseous, liquid and 
solid phase, these changes cannot be due to a 
change of the molecular structure. They must be 
connected with the fact that in liquids the mole- 
cules cannot rotate and translate as freely as in a 
gas. In a qualitative way the apparent decrease 
of the optical anisotropy can be understood by 
considering the anisotropic Lorentz force. Raman 
and Krishnan' have given a theory of this effect, 
but Stuart and Volkmann? have criticized their 
assumptions and we shall point out the inade- 
quacies of their arguments. 

Recent work has shown that in a liquid the 
positions and orientations of the molecules are 
not random. The structure of a liquid can be de- 
scribed by considering an arbitrarily chosen 


1C. V. Raman and K. S. Krishnan, Proc. Roy. Soc. A117, 
589 (1928); Phil. Mag. 5, 498 (1928). 
(19333 A. Stuart and H. Volkmann, Zeits. f. Physik 83, 461 


wn 


plained by considering a moderate anisotropy of the 
Lorentz force and a potential barrier of finite height. In 
this case the molar refraction is slightly smaller than that of 
the ideal gas. No large deviations from the Lorentz- 
Lorenz equation can be expected even if the anisotropy is 
considerable. The refraction changes with pressure and 
temperature and its dispersion differs from that of the gas. 
The depolarization and the intensity of the scattered light 
are considerably smaller than, and the Kerr constant is 
only a fraction of, the values derived from the classical 
equation. The classical relations between depolarization 
and intensity and between Kerr constant and depolariza- 
tion are only approximate. The theory is applied to benzene. 
It is assumed that the molecules have the same polariza 
bilities in the liquid and gaseous state. All observations 
can be explained by assuming a potential barrier of about 
8.8 kT and an anisotropy of the Lorentz force which is in 
good accordance with the results of the x-ray analysis of the 


liquid structure. 


molecule. Due to the dense packing and the 
intramolecular forces the nearest neighboring 
molecules tend to assume definite preferred posi- 
tions with respect to the molecule considered. 
They form a quasi-crystalline or cybotactic 
group. For anisotropic molecules these groups 
are anisotropic. Hence the neighboring molecules 
exert forces on the central molecule which pre- i 
vent its free rotation. Quasi-crystalline grouping 
and hindered rotation are closely interrelated, 
but the calculation is*simplified by considering 
their effects separately. 


THE INFLUENCE OF HINDERED ROTATION 

We consider a molecule with axial symmetry. 
b,, bg =}; are its polarizabilities. If the molecule is 
polar its permanent moment uz is parallel to the 
axis of symmetry. Following the method of 
Debye’ we can assume that the axis of each mole- 
cule has at every instant a definite preferred 
direction. The energy V(Q) is required to deflect 
the axis by an angle 9 from its preferred position. 
For nonpolar symmetric molecules, e.g., COs, 
C.He, the potential hill is symmetric with respect 
to two minima at 0=0 and O@=r7. For more 


3 P. Debye, Physik. Zeits. 36, 100, 193 (1935). 


+ 
~ 





548 HANS MUELLER 














— — 








2 4A lo js 


i 





Fic. 1. The functions go(Vo/RkT) and gs(Vo/kT) for a sym- 
metric potential barrier V = Vo(1—cos? 6). 


complex molecules and probably for most polar 
molecules © =0 gives the only minimum. Using 
classical statistics one finds the average polariza- 
bilities 8 of the molecule. If the field is acting 
parallel to the preferred direction 


B, = bo+ (b; — be) go+ (g2— 217) u?/RT (1) 
and for a field normal to the preferred direction 


B2= Bs = 3[(b1 +2) — (b1 — 2) g2 
+(1—ge)y2/kT). (1’) 


Here we have introduced a set of functions 


| exp (— V/kT) cos" @ sin 0d0 


g,(V/kT)= . (2) 


Sn 
ot 


exp (— V/kT) sin OdO 


ei) 


If ’=1\(1—cos @), g, can be expressed in terms 
of Langevin’s function L(x)=cotgh x—1/x, 
2i=L(Vo/kRT), go=1—2L(Vo/RT)RT/ Vo. For a 
symmetric potential hill go,,:=0. The simplest 
example of a symmetric potential barrier is 
V = V.(1—cos? ©) which gives 


1 


22= (3s hexpx { exp rat) —$x', 
“0 


+ 


x=Vo/kT, and gy=g2—(3g2—1)/2x. 


The latter two functions are plotted in Fig. 1. go 
has the limiting values } for free rotation and 1 
for no rotation. g,; varies between 0.2 and 1. 

In the classical theory it is assumed that the 
field acting on any molecule is parallel to the ap- 





plied field E and has the value E+4rJ/3, where I 
is the polarization. Since all angles between the 
field and the preferred directions occur with equal 
probability the average molar polarization P can 
be calculated in the usual way, giving 


P=(e—1)M, (e+2)p=(42 3)A4(81+2£62) 
=(4r 9)A[(b,+2b2)+ (1 —g;*)? kT |. 


Hence with the classical assumption the molar 
refraction R=(47r/9)A(b,+2b.) will not be 
changed whatever the form of the potential 
barrier may be. Furthermore, if the potential hill 
is symmetric, the dipole-contribution also will not 
be altered. Hindered rotation reduces the dipole 
term only if (©) is asymmetric. We shall show 
below that analogous conclusions hold also for 
the Kerr constant and the scattering of light. 

Experiments, however, show that neither polar 
nor nonpolar liquids satisfy the Lorentz-Lorenz 
equation. According to Stuart and Volkmann® the 
differences between the theoretical and the 
measured values of the Kerr constant and the 
depolarization are equally large for either type of 
molecule. Hence Debye’s theory of hindered rota- 
tion, while it may explain the dielectric constant 
and the Kerr effect of some special polar liquids, 
is not sufficient to account for the optical prop- 
erties of most liquids. 

Our further considerations are confined to non- 
polar liquids. For these the classical assumption 
F=E+4rlI/3 leads always to the Lorentz- 
Lorenz equation. Hence it is necessary to abandon 
this assumption and replace the isotropic by an 
anisotropic Lorentz force. 


THE ANISOTROPIC LORENTZ FORCE IN CRYSTALS 


Inasmuch as the structure of a liquid has a 
closer resemblance to that of a crystal than to 
that of a gas, it is natural to approach the prob- 
lem of the inner field by considering first an ideal 
crystal. The theory of the Lorentz force in a 
lattice was developed by Ewald‘ and generalized 
by Born.’ The writer® has given a simplified 
derivation. 

We consider an infinite lattice whose elemen- 
tary cell contains only one atom situated at the 


*P. P. Ewald, Ann. d. Physik 49, 1, 117 (1916). 

5M. Born and M. Goeppert-Mayer, Handbuch der 
Phystk, Vol. 24, p. 770. 

®H. Mueller, Phys. Rev. 47, 948 (1935). 
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OPTICAL 


origin of the cell. If x, y, s are the axes of the 
optical index ellipsoid, the Lorentz force acting 
on any atom can be written 


F,=E,+4L,l,. (3) 


The ‘Lorentz factors” L,, L,, L. can be derived 
from Ewald’s lattice potential’ 


y=y-—1/r, L.=A/41(0p/ dx?) 9. 
Poisson's equation requires 
L,+L,+L.,=1. (4) 


The Lorentz factors depend on the shape, but not 
on the size of the elementary cell. Previously the 
author’ has given the values of L,=L, and L, for 
tetragonal lattices. Fig. 2 gives these factors for 
simple hexagonal lattices. They depend on the 
ratio c/a, where c is the lattice spacing in the 
direction of the hexagonal axis, and a is the 
distance to the nearest neighbors normal to the 
axis. The results are very similar to those for 
tetragonal lattices. 

We observe again that if c/a is moderately 
large L,=kce/a and L,=1—2kc a. with 
¢ a>1 we have a layer lattice, this result has the 


Since 


following simple interpretation: In a layer lattice 
the Lorentz force can be calculated by consider- 
ing the atoms continuously distributed in all 
layers except the layer in which the atom is 
situated. For a field E, the continuous distribu- 
tion gives the field 47/7. and hence contributes 1 
to the factor L,. The dipoles in the same layer 
give a negative field which can readily be calcu- 
lated by direct summation. Thus one finds 
k=0.359 for tetragonal and k=0.370 for hex- 
agonal lattices. These factors agree with the 
slope® of the L.-curves if c a>1.5. 

Similarly for small values of c/a one has 
L.=k'(a/c)*?, which means that in a chain lattice 
it is necessary to consider the discontinuous struc- 
ture along only a single chain of atoms. k’ is 
therefore given by a simple infinite sum and has 
the values 0.3826 and 0.3313 for tetragonal and 

? The lattice potential y is the electrostatic potential of 
a charge distribution consisting of discrete unit positive 
charges at the lattice points and a uniform constant nega- 
tive charge density —1/A, where A is the volume of the 
elementary cell. If the charge at the origin is missing the 
potential is y=y—1/r and the derivative must be taken 
for the values x = y =z =0. 

* In the paper mentioned (reference 6) we estimated from 


the graph for ‘tetragonal lattices 2k =0.72. The direct sum- 
mation gives 0.718. 
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Fic. 2. The Lorentz factors for simple hexagonal lattices 
and the approximations for small and large values of c/a. 
Dotted curve: Raman and Krishnan's approximation. 


hexagonal lattices respectively. This is a good 
approximation if c a<0.6. If ca differs little 
from unity one can use a similar approximation 
by calculating the field of the dipoles within a 
rectangular box. The validity of the latter 
method was demonstrated in a previous paper.* 

In evaluating the field of the dipoles in a single 
layer or chain one again can consider the atoms at 
a large distance as continuously distributed. 
Hence the Lorentz factors are determined if we 
know the distribution of the few nearest neigh- 
bors only. Even in a crystal the atoms more than 
a few atomic distances away can be looked upon 
as forming a continuum. 

This result justifies the well-known cavity 
method for calculating the inner field® and fur- 
nishes the basic principle for evaluating the 
Lorentz force in liquids. 


illustrates another important point. 


Fig. 2 
Following a suggestion of Havelock,”” Raman and 
Krishnan! assume that for an anisotropic struc- 
ture an ellipsoidal cavity can be found which has 
the property that the atoms within give no con- 
tribution to the Lorentz force. This assumption 
is incorrect because L can take negative values, 
whereas the cavity method leads always to posi- 
tive values. Without justifying themselves Ra- 
man and Krishnan postulate furthermore that 
the hypothetical ellipsoid has an axial ratio c/a. 
Hence they assume a Lorentz factor 

’ This result shows also that in general it is advisable to 


use an anisotropic cavity for an anisotropic crystal. 
‘0 T, H. Havelock, Proc. Roy. Soc. A80, 31 (1908). 
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where ¢ is the eccentricity @=1—(a/c)*. If a>c 
one takes & = (a/c)*—1 and obtains 


1+é 
$,°-—— (1-a e) tan '). (5’) 
é 


s, and s, are determined by the relations s,=s, 
and s,+s,+s,=1. In Fig. 2 the dotted curve 
gives the values of s.. Raman and Krishnan’s 
assumption s,=L, is obviously untenable. 

The use of the s factors is however justified if 
the cavity method is correctly applied. The 
cavity can have any size as long as it includes at 
least the nearest 20 or 30 neighboring atoms. We 
find then F,= E,+4rs,I,+ F,’. Comparison with 
(3) gives 


F,' =4r(L,—s,)I,=4ry-lz, (6) 
where y.+7,+7:=0, hence y=y:=—2y7:= 
—2y,=L:-S:. 


F’ is the field of the finite number of dipoles 
within an ellipsoidal cavity, the axes of which are 
parallel to the optical axes of the crystal. The 
field F’ does not depend on the size of the ellip- 
soid, but only on its eccentricity. Consequently 
the field of the dipoles within a uniformly 
polarized shell, bordered by two concentric 
ellipsoids of equal eccentricities is zero at the 
center. 

The above results hold also for molecular 
lattices. In a forthcoming paper it will be shown 
that they are also valid for “‘diffuse”’ lattices in 
which the atoms or molecules oscillate about their 
equilibrium positions. 


THE LORENTz FORCE IN LIQUIDS 


From the work of Stewart,' Warren,” Bernal 
and Fowler," Debye,’ Frenkel" and others, the 
following assumptions concerning the structure of 
liquids seem reasonable. For rod- or disk-shaped 
molecules the preferred positions of the nearest 
neighboring molecules are arranged in a lattice 
with axial symmetry. In many liquids the x-ray 





1G, W. Stewart, Trans. Faraday Soc. 29, 982 (1933). 

12 B. E. Warren, Phys. Rev. 44, 969 (1933). 

18 J. D. Bernal and R. H. Fowler, J. Chem. Phys. 1, 515 
(1933). 

44 J. Frenkel, Acta phys. chem. 3, 663 (1936). 
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analysis” points to a simple hexagonal lattice." 
The preferred orientation of the molecule con- 
sidered and of its nearest neighbors will have the 
direction of the symmetry axis of the lattice. The 
position and orientation of any distant molecule 
is arbitrary. 

The environment of a molecule can be divided 


into three regions: 


1. The cybotactic region, with a pronounced 
lattice structure. Within it all molecules have the 
same preferred orientation. The dimensions of 
this zone must be much smaller than the wave- 
length of light. We assume that it includes at 
least 20 to 30 molecules and that it has the shape 
of an ellipsoid of rotation whose eccentricity € is 
smaller than the eccentricity of the ellipsoid 
representing the shape of the molecule. 


2. The continuous region. This includes all dis- 
tant molecules whose positions and orientations 
are independent of those of the molecule upon 
which we fix our attention. As far as this molecule 
is concerned the distant molecules act as if 
forming a continuous distribution of matter. This 
region, extending to infinity, starts at a certain 
boundary, which is assumed to be an ellipsoid, 
larger but of the same eccentricity as the one 
limiting the cybotactic zone. 


3. The intermediate region between 1 and 2. Here 
the molecules are only slightly influericed by the 
intramolecular forces exerted by the central 
molecule. This region may be divided into a series 
of concentric ellipsoidal shells of constant ec- 
centricity «. Every molecule within one shell will 
be influenced in the same manner by the central 
molecule. 


The molecule upon which we fix our attention 
may be any arbitrary molecule. The symmetry 
axis of its cybotactic lattice can take on any 
direction with equal probability. We assume that 
no outer influence, not even that of a strong 
electric field, can change the distribution of the 
preferred directions to any marked degree.'® The 
preferred orientation of a molecule changes 
slowly with time, but this change is assumed to 


5 This is the simplest possible assumption, but it is not 
essential for our calculation. 

6 Investigations of the light scattering and the x-ray 
diffraction of liquids in a strong electric field have furnished 
no evidence to the contrary. 
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be negligible during the period of a light wave.'? 

Under the influence of an electric field E the 
liquid acquires an average polarization J in the 
direction of E. Within a cybotactic group, how- 
ever, the polarization is not in the direction of the 
field. Hence we introduce the concept of the 
“local polarization J.’’ J varies continuously 
throughout the liquid. It depends on the local 
preferred direction. 

The results of the preceding section give for the 
Lorentz force in a liquid 


~ 
— 


F,=EFE,+49s,I1,+4ry-,J-. ( 


4rs,I, is the field of the distant molecules which 
act as if they were a continuum. 47y,J, is accord- 
ing to Eq. (6) the field of the dipoles in the cybo- 
tactic group. The dipoles in the intermediate 
range give no field, because, as we have seen, the 
field of a uniformly polarized ellipsoidal shell is 
zero at the center. The three equations (7) give 
the components of F parallel and normal to the 
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preferred axis. Hence 


J,= NB:2F;; Jy=NB2F,; J.=NB6,F,, (8) 


where 8, and #2 are the average polarizabilities of 
the molecule defined by Eq. (1). N=pA/M is the 
number of molecules per cc, p, the density, A, 
Avogadro's number and M the molecular weight. 


Introducing Eq. (8) in Eq. (7) gives 
F,=1(E+4/3rl),, (9) 
where, since J =(n?—1)E/4r, 


1,=1,=[1—3/2(s—4)(m?—1)/(n?+2)] 
(1+2ryNB8z2), 
1,=[1+3(s—4)(m?—1)/(n?+2) ] 


(1—4ryNB;). (10) 


By combining Eqs. (9) and (7) one gets three 
equations for J,, J, and J,. Finally, since J is ob- 
viously the average value of the local polarization 
J, I=J, one can eliminate J and calculate the 
molar refraction R, of the liquid. 


} 1T By (1 —4ryB8,N)+2p2 (1+2ry62N) | 


R,=(n*—1)AM,(n*+2)p=-27A 
3 


This expression can be simplified. The denom- 
inator corresponds to the correction given by 
Raman and Krishnan.'* But since the boundary 
of a cybotactic group has a smaller eccentricity 
than the molecule, our value of (s—34) is much 
smaller than theirs. In fact this correction is so 
small that it can be neglected in most cases. 
Therefore we take s=}, which means that the 
cybotactic groups are approximately spherical. 
Even if they have elliposidal shape, the assump- 
tion s= 4 does not alter the results by more than a 
fraction of one percent, because a change of s 


1? The characteristic frequencies of mutual vibration of 
the molecules are in the far infrared. C. H. Cartwright, 
Phys. Rev. 49, 470 (1936). 

'*® Raman and Krishnan’s relation is obtained if we put 
s=L,, hence y=0 and 6,=);, B2=be. The latter relations 
imply that the molecules have entirely lost their freedom 
of rotation (g.=1) Raman and Krishnan’s equation gives 
results of the correct order of magnitude because it con- 
tains two errors which partly compensate each other. 
The anisotropy of their Lorentz field is too small (s<L,, 
see Fig. 1), whereas the effect of hindered rotation is taken 
too large. 


(11) 


1—(4 3)aN(s—})[Bi (1—4ry8,N)—Be (1+2ry62,N)] 


alters the value of y=Z,—s in such a way, that 
l,, 1, and R;, remain almost unchanged. 
If we introduce 


p=b, bo, b=4(b,+2d2), 
B=[(1—g2(1—p) ]/(p+2), (12) 
x=129rNb(L,—4)=9(L.,—4)pRo/M, ” 
Ry =4/3xAb. 
Eq. (11) takes the simple form 
R,_/Ryo=B/(1—xB) 
+(1—B)/(1+4x(1-—B)). (13) 


Ro is the molar refraction of the gas, x depends on 
the anisotropy of the quasi-crystalline structure 
and B on the freedom of rotation. Since };, b: and 
Ro can be obtained from measurements in the 
gaseous state, Eq. (13) contains only two un- 
known parameters L, and ge. They are charac- 
teristic of each liquid. ge varies with the tem- 
perature and the density. We do not expect a 
temperature dependence of L., except near the 
critical temperature of the liquid. 
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Fic. 3. The molar refraction of liquids. The variable 
x =9(L.—1/3)pR )/M depends on the anisotropy of the 
Lorentz force, the parameter B on the potential barrier 
and the optical anisotropy of the molecule. The dotted 
curve gives the minima of the curves 


THE MOLAR REFRACTION OF LIQUIDS. BENZENE 


The classical work of L. Lorenz led to the con- 
clusion that the molar refraction is the same in 
the liquid and gaseous state. This fact has fre- 
quently been mentioned as an argument against 
anisotropic Lorentz forces. But more accurate 
measurements have shown that the refraction of 
a liquid is usually slightly smaller than that of its 
vapor, and that it varies with temperature and 
pressure. In addition the dispersion of Ry, is 
larger than the dispersion of Ro. 

Although these differences are relatively small, 
they are nevertheless too large to be explained by 
changes of the polarizabilities of the molecule.'® 
Our theory accounts for all these changes without 
postulating any alteration in the polarizabilities. 

In Fig. 3 Rz/Ro is plotted as a function of x 
for various values of B. Since x is of the order of 
magnitude of 9(L —3)(n?—1) /(n?+2) its absolute 
value is probably never larger than 3. From Fig. 
2 it follows that x is positive for disk-shaped and 
negative for rod-shaped molecules. The param- 
eter B varies between } for ge=}, and p/(p+2) 
for ge=1. Since for most molecules } <p<3 it is 

1° The work of Fajans and Joos (Zeits. f. Physik 23, 1 
(1924)) and the theory of photoelasticity (H. Mueller, 
Phys. Rev. 47, 947 (1935)) show that such changes must 
occur in solids, but in liquids they cannot play an im- 
portant part as shown by the fact that the molar polariza- 
tion of CO, increases with pressure (F. G. Keyes and J. G. 
Kirkwood, Phys. Rev. 36, 754 (1930)) and that the polariza- 


tion of solid benzene is larger than that of benzene vapor 
(C. P. Smyth, J. Am. Chem. Soc. 55, 453 (1932)). 
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sufficient to vary B between the limits 0.15 and 
0.6. Rod-shaped molecules usually have p> 1 and 
hence B>}4, while for disk-shaped molecules 
B<}.™ 

The curves in Fig. 3 have a minimum 9B(1 — B) 
(3B+1) for x=2(1—3B) 3B(1—B) and R_/ Ry 
is 1 for x=0 and for x=(B+1)(1—3B)/B(1—B). 
If the molecules can rotate freely or almost freely 
the anisotropy of the Lorentz force demands 
R,_>R>. It is probable that this condition exists 
in compressed gases, which would explain the 
increase of the molar polarization of CO: with 
pressure.'® If the anisotropy is very large we 
always get R,;>Ro, even if no rotation exists. 
The Lorentz force reaches its largest anisotropy 
in the crystalline state. Hence it is not surprising 
that the molar polarization of solid benzene is 
larger than that of its vapor.'® 

The case R, < Rp occurs only if we have hind- 
dered rotation together with a moderate ani- 
sotropy of the Lorentz force. These are exactly 
the conditions existing in liquids. We notice that 
in all cases of practical importance R,,/ Ry >0.9. 
Hence the approximate validity of the Lorentz- 
Lorenz equation cannot be advanced as an 
argument against anisotropic Lorentz forces. 

Since we have no direct information concerning 
the values of go and L, we are not able to predict 
the refraction of a liquid. Inversely, however, we 
can find approximate values of these parameters 
from the observed values of R; 

We will discuss the procedure for benzene. 
According to Stuart” the polarizabilities are 
b, = 63.5 10-*, bp = b3 = 123.1 10-* hence p=0.516. 
In Fig. 4 R,/Ro is given as a function of x for 
various values of ge. Wasastjerna” has deter- 
mined Ry=27.20 (D-line) and R;,=26.18 (room 
temperature), hence R,,/ Ry=0.962. Fig. 4 shows 
that this value can be understood only if g2 
is larger than 0.775. If we assume a _ potential 
Vo(1 —cos? @) this indicates, according to Fig. 1 

20 This is not true in some exceptional cases, e.g., formic 


acid and acetic acid. For these substances our theory 
predicts R.i>R». 

1H. A. Stuart, Molekiilstruktur (Berlin 1934), p. 173. 
These polarizabilities are determined from the Kerr effect 
of the vapor. Cabannes’ measurement of the depolarization 
of the scattered light gives p=0.505. 

2 J. A. Wasastjerna, Soc. Sc. Fennica, Comm. Math. 
Phys. 2, No. 13 (1924). The corresponding values of the 
molar polarization are P; = 26.61 P)=27.01, hence P;,/P 
= 0.985. The ratio pp of the static polarizabilities is probably 
larger than p. Psotia = 28.5. 
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that the potential barrier must be higher than 
5.2 kT. 

A value for x and a closer estimate of gs can be 
obtained if we consider the dispersion. A change 
of wave-length does not affect ge and L., since 
they depend only on the liquid structure, but it 
does change the value of Ro and hence of x. The 
slope of the curves R, Ro for g2=const. is related 
to the difference between the dispersions of R, 
and Rp». According to Wasastjerna” Rr—Re 
=().74 for liquid benzene, and Rr — Rc = 0.68 for 
the vapor. From these data one estimates 
d(R,, Ro) dx=0.12/x. 
these slopes are shown in Fig. 4. We see that there 
is only one solution for which R;,/ Ro =0.962 and 
for which the slope suits the observed difference 


Lines corresponding to 


of dispersions, namely x = 1.8 and g2=0.875. This 
corresponds to a potential hill of about 8.8 kT, a 
value of the same order of magnitude as those 
given by Debye* and Andrade.* 

With p=0.881, Ry=27.2, W=78.05, we calcu- 
late from x=9(L.—4)pRo/ M (Eq. 12) L.=0.98. 
If we assume a simple hexagonal structure of the 
cybotactic group we find from Fig. 2 that this 
corresponds to c/a=0.56. According to the con- 
siderations of de Boer** a hexagonal structure 
seems a good approximation. He reports that 
recent x-ray work gives c=3.5 to 3.6A. To ac- 
count for the observed density we must have 
M /Ap=ca’*3?/2 or c/a=0.50 to 0.53, which is a 
surprisingly good verification of our result. 
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Fic. 4. The molar refraction of benzene. The dotted line 
gives the observed value. The short lines indicate the slopes 
d(R,/R.)/dx required by the observed dispersions. 


In view of the fact that the value of p is 
accurate to only about 3 percent and that a small 
part of the change of dispersion may be due to the 
shift of the absorption lines,” the above results 
are somewhat inaccurate. The true solution prob- 
ably lies within the ranges 1.4<x<1.8 and 
0.82 <g2<0.9. 

With this solution we can now estimate the 
order of magnitude of the temperature—and the 
pressure variation of R,. If we assume L,=const. 


and abbreviate R,/Ro=r, ¢= Vo(p)/kT we have 


OR, 1dpf Or dr dgo dV» p Or dgs ¢ 

1 Ry = E ¢ |- : (14) 
dT pdTt ax dgedge dp Vy 0ge dg T 
OR, 1dpf dr Or dgs dVy p 

1/Ru——= v—+4 e | (15) 
Op p dp Ox Og,do dp V;, 


With the values (1/p)(dp/d7T)=—1.21 10-3, 
(1/p)(dp/dp)=9.5 107", 
= —0.44 (the latter value estimated from Fig. 4) 
dg2/dg=0.018 (from Fig. 1), ¢=8.8, T=290 we 
get roughly (1, R))dR,/ dT =2 to 3 10+. The ex- 
perimental value is 1.22 10-*. The theoretical 
value of (1, Ro)dR» dp depends largely on Vo(p). 
If we assume” |’)=cp' close agreement with the 
observed value —0.98 10-" is obtained for ¢=3. 


*3E. N. da C. Andrade, Phil. Mag. 17, 497 (1934). 
4]. H. de Boer, Trans. Faraday Soc. 32, 13 (1936). 


xdr/dx=0.12, Or/dge 


But this result depends very much on the value of 
dg2/dg and hence on the choice of the potential 
barrier. At any rate the order of magnitude of 
9R_/dT and 0R_/dp is correct. 

To sum up, the value of the molar refraction of 
liquid benzene, its dispersion, and its pressure 
and temperature dependence can be explained by 
rotation and cybotactic 


assuming hindered 


grouping. The characteristic quantities of CyH, 


* P. Pringsheim, Handbuch der Physik, Vol. 23/1, p. 246. 
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are go=0.88 and x=1.8 or L.=0.98 approxi- 
mately. These quantities can be determined with- 
out making any special assumptions concerning 
the type of grouping or the form of the potential 
barrier. The assumptions are introduced only to 
correlate the parameters with quantities which 
have a simple physical meaning, and to show that 
their values are reasonable. 


THE INTENSITY AND DEPOLARIZATION OF 
THE SCATTERED LIGHT 


The theory of the unmodified scattered light 
for liquids was developed independently by 
King,”® Gans*’ and Ramanathan.” Free rotation 
and isotropy of the Lorentz force were assumed. 
In their original forms the three theories gave 
different results, but the later papers of Ramana- 
than and Didlaukies** have eliminated these dis- 
crepancies. However, the large amount of experi- 
mental data available do not verify the theory.” 

We will follow Ramanathan’s general formula- 
tion of the theory, but correct it by taking into 
account the influence of hindered rotation and 
the anisotropic polarization field. If the incident 
light of intensity 1 is plane polarized with an 
electric vector E, the intensity of the light 
scattered by a unit volume at a distance 1 cm in 
the direction normal to the plane of E is given by 


T= 3(24/d)4*N[Ly (utp) av? + (up?) av — (Mp) av? 
I (16) 


9 


n= 3(24/d)*N (un?) av. 


I, and J, are the scattéred light components hav- 
ing electric vectors parallel and normal to &£, 
respectively. \ is the wave-length in vacuum, V 
the number of molecules per cc, y= NkT8, k 
Boltzmann's constant and £6 the isothermal com- 
pressibility at the absolute temperature 7. yp» is 
the component of the induced dipole moment in 
the direction of E and uy, is its component normal 
to E. (u) ay? is the square of the average and (u”) a, 
the average of the square. To find these averages 
we introduce three Cartesian systems of co- 
ordinates (Fig. 5): A fixed system p, m1, m2 whose 
directions coincide, respectively, with £, with 
the direction of propagation of the incident beam 
and with the direction of observation. The second 
system x, y, 2 has the g axis parallel to the pre- 
ferred direction of a molecule with the x axis in 
the plane (pz). The position of this system is 
given by the angles #, ¢. Since the preferred 
directions are distributed at random the prob- 
ability of finding a molecule with the preferred 
axis in the range dd, dg is sin ddddg/4r. The 
third system of coordinates 1, 2, 3 coincides with 
the axes of the ellipsoid of polarization of the 
molecule. The orientation of the symmetry axis 1 
with respect to the axis z is given by the polar 
coordinates ©, y, the axis 2 being in the plane (21). 
The probability that the molecule has this 
orientation is exp(—V/kT) sin 6dédy/ f exp 
(—V/kT) sin @dédy. The field acting on a mole- 
cule is given by (9). By calculating successively 
| ai mst F,, Fo, F;; Mi, M2, Ma; Mey My, Me ONE 
finds 


(lL (dy sin? 6+, cos? @) cos? y sin? 8+(b;—bde2) cos 6 sin @ cos 8 sin J cos ¥ 


up=(E+4 3x1)! 


+z sin® p sin? 3 ]+/.[(b; cos? 6+, sin? 6) cos? 3 


+(b,:—b:2) cos @sin cos d sin d cosy], (17) 


1,[(b, — be) sin @ sin 3 cos (cos 6 sin 8 cos gy—sin 6 cos 8 cos ¢ cos ¥ 


ttn, = (E+4/3x) 


+sin @sin ¢ sin Y)—bd» sin J cos J cos ¢ | 
(18) 


+1,[(b; cos? 6+. sin? @) sin 3 cos 8 cos ¢ 


+(b:— 2) sin @ cos @ cos 3(sin ¢ sin Y—cos J cos ¢ cos p) }. 


These expressions must be introduced in 


(up )av —_ 


Su, exp (— V/RT) sin 6 sin ddéddd ody 
> (19) 
S exp (— V/kT) sin @ sin ddéddd gdy 

” L. V. King, Proc. Roy. Soc. A104, 333 (1923). 

27 R. Gans, Zeits. f. Physik 17, 353 (1923). 

*8 K. R. Ramanathan, Proc. ind. Assoc. f. cult. sc. 8, 181 
(1923); Indian J. Phys. 1, 413 (1927). 





29M. Didlaukies, Ann. d. Physik 5, 205 (1930). 

%° The difficulties have been discussed by J. Cabannes: 
La diffusion moléculaire de la lumiére (Paris, 1929), pp. 
204-300. 
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and the analogous expression for (u,”) av. (u»)av is of course identical with the average polarization 
per molecule J/N and (u,),v=0. The integration of (19) gives 


(u,?)av = ((E+4n1/3)?/15 


(un®)av = ((E+41/3)2/15]} 


ee ee ene 
+1.2[b:°(3 —5Sg2+2gs) +b2"(3-+3g2+ 2gs) + 2bibo(1 +22 — 24) ] 
| + el .[ 2b1*(g2—g4) +2°(1 +g2—2¢4) +bib2(1 —3g2+4g4)] (20) 
Saati 1 — gy) +2b:b2(gs—ge) ] 
+1?[ bP (1 — g4) +b2?(1 +g2+84) +biba(— 1 — g2+2g4) J 


{ +/,1,[by?(2gs— 2g2) +2?(—1—ge+gs) +bi1b2(—1+3g2—4¢5)] (20’) 


where go and g, are the functions of the potential 
barrier defined in Eq. (2). If we take again s=}, 
the definitions (10) give 


1,=1/(1+x«(1—B)/4), 1,=1/(1—xB), 


where x and B are defined in (12). 
If there is no anisotropy of the Lorentz force, 
x=0, 1,=1,=1, we get the classical relations 


(uy?) av = ((E+4nJ/3)?/15 [3012+ 8b? +45 yb |, 
(un?) av= [(E+4aI/3)?/15 ](b; —be)?. 


This is true for any value of ge and g;. Hence 
hindered rotation alone does not alter the 
classical theory. 

If the incident light is not polarized, one has 
I,’ =I,+T],, I,’=2I, and hence the total inten- 
sity of the scattered light is Jo=J,+3J,, and its 
depolarization A=2/,/(J,+/,). By introducing 
(20) in (16) and using the definitions (12) we get 
To =3(m/d*)?((m? —1)?/N) +2; 

t=[y—14+(4,4+3M2)3/5]; (21) 
A=6M2/[5(y—1)+3(Mi+M2)], (22) 
where 
M,=(12m,4+12m2 
+ 21 .1,m;)(Ro/Rx1)?/(p+2), (23) 
Ms, = (1.2ms3 ; 
+12ms—1,1,m;)(Ro/ Rx)?/(p+2)* 


and 


m,=3+ ¢2(p?+4p—5)+2gs(p—1)?, 
m2= (3p? +2p+3)+g2(—Sp?+2p+3) 

+22:(p— 1)?, 
ms=1+2g2(p? — p) — gs(p—1)?, (24) 
ms=(p?—p+1)+g2(1—p) —ga(p—1)?, 
m;= (1+ p)+g2(2p?—3p+1) —2g.(p—1)?. 


For s=0 we obtain the classical relations 
A=68?/(5y+78) 
where 6&=(p—1)?/(p+2)?=(b;—b2)?/ (6: +262)? 
and t=7+138/5, 
which lead to 
i=6y(1+A)/(6—74). (25) 


For the case of an anisotropic Lorentz field it is 
not possible to retain these relations by intro- 
ducing for example a quantity analogous to the 
anisotropy factor 6. In particular Eq. (25) is not 
justified.*" 

The significance and use of the results can best 
be illustrated by applying them to benzene. 
Since p=0.516 is given, one constructs first a 
table of B and m; for a series of values of go. We 
assume again a potential 9(1—cos*® @) and ob- 
tain for each g2 the corresponding® value of gy 
from Fig. 1. Next one tabulates /],, 1, and 
(R_/Ro)* as functions of x and ge and using (23) 
gets M,(x, go) and M2(x, ge). Finally (21) and 
(22) serve to find the intensity factor i(x, ge) and 
the depolarization A(x, g2). The results are 
plotted in Figs. 6 and 7. They are constructed 

%! The determinations of Avogadro’s number from the 
light scattering in liquids presupposes the validity of Eq. 
(25). This relation is true only if 3M,—4M,=5. In our 
theory (3M,—4M,)/5 is a function of ge and x. The calcu- 
lation for benzene shows that this function is very similar 
to R;,/Ro shown in Fig. 4. Its value never deviates much 
from 1. The fact that the light scattering method has 
furnished approximate, but not accurate values of Avo- 
gadro’s number is therefore not an argument against but 
rather for the anisotropy of the Lorentz force. 

8? Fig. 1 shows that (g2—g,) remains almost constant for 
V./kT <5 and decreases very slowly for higher potential 
barriers. We have verified the same behavior for other 
forms of the potential. Hence for small values of ge and 
also for g2:=1 our curves do not depend on the choice of V. 


Between g:=0.8 and g:=0.95, however, the results vary 
for different potential barriers. 
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Fic. 5. The 3 coordinate systems used in calculating light 
scattering and Kerr effect. p, m, m2 is fixed in space; m 
direction of propagation of incident light; z direction of 
preferred orientation; 1 momentary direction of symmetry 
axis of molecule. 


with y=0.026, which corresponds to the density 
and the compressibility of benzene at 18°. 

If there were no anisotropy of the Lorentz 
force we should find ip =12.0 10-? and Ay=0.57 
whereas the observations at 15° give* 7=8.2 10° 
and A=0.42. These observed values can be 
completely understood if g2>0.8, as Figs. 6 and 
7 show. Our previous results from the refraction, 
that ge is about 0.88 and that x is between 1.4 
and 1.8 are thereby given considerable added 
support.™ 

Recently Turner® has shown for benzene that 
the variation of the intensity of the scattered 
light with the wave-length is not strictly pro- 
portional to (m?—1)*,/A*. Our theory does not 
require this proportionality, because 7 depends 
on x and hence on Rp and the wave-length. 
From \=5460A to \=4350A Ro increases by 
about 3 percent, hence x changes by 0.054. For 
x=1.8, gs=0.88 the slope of the curve is about 
a1 dx=0.2, hence At=(di/dx)Ax=0.011 or an 
increase of i of about 14 percent. Turner’s data 


83Calculated from Cabannes’ value R=10.7 10~. 
Cabannes’ remark, that C,H, scatters more than it should, 
is based on the implication that (25) is correct. (25) gives 
i=7.3 10°. Actually benzene scatters much less than it 
should according to the classical theory. Hence we do not 
think that Rocard’s theory can be accepted. See J. Caban- 
nes, La diffusion moléculatre de la lumiére, pp. 237 and 296. 

54 Figs. 4 and 6 are not appreciably altered if one takes 
p=0.500 or if a different form of the potential barrier is 
chosen. The curves in Fig. 7 however are very sensitive to 
such changes and a better agreement for A could not be 
expected. 

% A. F. Turner, Diss. University of Berlin (1934), Physik. 
Ber. 16, 1860 (1935). The observed increase of i for wave- 
lengths \>5460A could only be explained if our curves 
had a large curvature at their minima or if p changes 
with X. 


give between these wave-lengths an increase of 7 
of about 20 percent. 

From the slope of the curves in Fig. 7 it 
follows that the depolarization should vary by 
less than 0.01 in the range of the visible spectrum. 
Within experimental errors of this magnitude 
the observations have shown no dispersion of 
the depolarization. 

If the potential barrier is low the anisotropy 
of the Lorentz force gives a larger depolarization 
than the classical theory. Hence the ‘‘apparent”’ 
optical anisotropy*® of the molecules can be 
larger than for the gas. This will occur in liquids 
at high temperatures and in gases at high 
pressures. In fact Rao*’ has shown that near the 
critical temperature the apparent anisotropy is 
from 3 to 4 times larger than for the gas, and in 
gases Volkmann**® has found an increase of & 
with pressure. 

Summing up, the theory accounts for all 
anomalies of the intensity and the depolarization 
of the scattered light. It is perhaps surprising 
that the theory of Gans,*?’ which is based on 
statistical considerations, does not lead to our 
results. This is due to the fact that Gans assumes 
that the isotropic micro-distribution corresponds 
to a state of minimum energy. In liquids it 
seems more natural to assume that in small 
volume-elements the optical anisotropy fluctu- 
ates around a finite value. 


THE KERR EFFECT IN LIQUIDS 


The measurements of the electro-optical Kerr 
effect demonstrate the failure of the classical 
theory of Langevin and Born in the most 
impressive manner.*® The Kerr constant K, of 
liquids is always much smaller than its theo- 
retical value Ko. The work of Briegleb,” and 


% The “apparent” optical anisotropy is the value of & 
which is obtained from the relation A = 68?/ (57 +-78*). 

87 R. S. Rao, Ind. J. Phys. 2, 7 (1924). Rao, Stuart and 
Volkmann, reference 2, believe that the large increase is 
not real but caused by experimental difficulties. In view 
of Volkmann's latest results, reference 38, 5 must have a 
maximum where & > 6,?. 

88H. Volkmann, Ann. d. Physik 24, 457 (1935). 

39C,. V. Raman and K. S. Krishnan, Phil. Mag. 3, 713 
(1927) claim the validity of the classical theory for liquids. 
But they compare the Kerr constants with the ‘‘apparent”’ 
anisotropy determined from the depolarization in liquids. 
Hence their considerations are not a verification of Lange- 
vin’s theory. They only demonsirate the approximate 
validity of the relation of Gans between K, and A. 

40G. Briegleb, Zeits. f. physik. Chemie 14, 97 (1931); 
16, 249 (1932). 
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Fic. 6. The intensity factor 7 for the scattered light for ben- 
zene. The dotted line corresponds to the observed value. 


Stuart and Volkmann? makes it seem almost 
impossible to correlate the Kerr constants of 
liquids with the optical properties of the mole- 
cule. 

There exist two possibilities for the formulation 
of a new theory. One might assume that the 
cybotactic groups behave like small anisotropic 
crystals. Even if they are spherical they will be 
oriented by an electric field, and this orientation 
produces uniaxial birefringence of the liquid. 
However it can easily be seen that this assump- 
tion gives much too large values of the Kerr 
constant." Although the groups have a smaller 
than the molecules, they 


themselves much 


optical anisotropy 
would orient more easily, 
because the energy of orientation is proportional 
to their volume, while the temperature energy 
kT is the same as for a single molecule. Supposing 
the groups contain only 50 molecules and are 
spherical a calculation carried out on this basis 
gives for benzene K,>10Ko, while the experi- 
ment gives K,<}Ko. 
Hence we conclude that the molecules in a 
quasi-crystalline group are not held together 
rigidly. Even if the potential barrier is so high 
that rotation is practically impossible the mole- 


w= fm exp (—(V+U)/kT) sin 6 sin ddéddd gdy 


* If the groups are large and not spherical the Kerr effect 
turns out to be exceedingly large and has a large relaxation 
time. This is the case in liquid crystals and in colloids. 

*” This statement does not invalidate our calculation of 
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Fic. 7. The depolarization of the scattered light for benzene. 


‘ 
The dotted line gives the observed depolarization. 


cules still have a considerable freedom of trans- 
lation. They can travel from one _ preferred 
position to another. The cybotactic lattice is a 
rather diffuse structure and does not act like a 
solid.” It is of course possible that the electric 
field changes somewhat the distribution of the 
preferred orientations, but we have no experi- 
mental data to estimate this effect.'® It must be 
small and we shall neglect it altogether. 

We calculate therefore the Kerr effect in the 
usual the 
action of the field on the single molecules. We 


manner® by considering orienting 
introduce again the three systems of coordinates 
(Fig. 5). The axis p is in the direction of the 
static field Ey between the plates of the Kerr 
condenser, ”; is in the direction of the light beam. 
The two other systems are defined as in the 
preceding section. Due the field Ey the 
molecules have an additional potential energy 


U= [ai Fo? + do Fos? + a3 Fo3" 1, 


to 


(26) 


where 4, de, ds are the static polarizabilities of 
the molecule, Foi, Fo2, Fos the components of the 
Lorentz force of Eo. In the system x, y, z the 
components of this field are Fo, =f.(Eo+421o/3)., 
where we place f,=1/(1+x(1—By)/4), f.=1 
(1—xBy), Bo= [1 — go(1 — po) | (pot 2), Po=41/ Qe, 
in analogy to the definitions of /,, /., B and p. 
For a light wave whose electric vector E, is 
parallel to the field Ey the average polarizability is 


(E,+(4 3)r1,) { exp (—(U+V) kT) 


Xsin Osin ddédddgdy. (27) 


the Lorentz force, because it can be shown that Ewald’s 
theory is applicable to ‘diffuse’ lattices, in which every 
atom oscillates about a preferred position. 


43 P. Debye, Handbuch der Radiologie, Vol. 6. 
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Hy is given by Eq. (17). For a light wave whose electric vector E,, is normal to the field Ey the average 
polarizability of a molecule is also given by (27) if we replace u, by wu, where 
L.[ (61 — 62) sin @ cos 6 sin 8 cos ¢(sin ¢ sin ¥Y—cos #8 cos ¢ cos W) 
+ (b; sin? 6+ 5:2 cos? @)(cos? 3 cos® ¢ cos? ¥+sin? ¢ sin? y 
Mn 

= —cos J sin g cos ¢ sin ¥(1+cos Y))+2(cos? 3 cos? ¢ sin? Y+sin? ¢ cos? p 
(FE, + (4/3)rTn) 
+2 cos 3 sin ¢ cos ¢ sin ¥ cos ¥) ]+/,[ (1 cos? 6+, sin® @) sin? 8 cos® ¢ 

+(b,— 2) sin @ cos @ sin 3 cos ¢g(sin ¢ sin Y—cos #8 cos ¢ cos p) ]. 
The rest of the calculation follows the usual method.“ One uses exp (— U/kT)=1—U/kT, finds 
the term proportional to E,*, and from y,=y+dy, one determines the indices of refraction n,=n 
+dn,, n,=n+dn,, where dn,=(4/3)xrN(n?+2)*dy,/6n.“ The resulting Kerr constant can be 





written in the usual form 


K,=(n,—n,)/nE? = (4/3)eN[(n?+2)/2n PL (e+2)/3 20, (28) 


but the value of ©, is now given by 
l.f2[a 1b; 2944+ 2arbo( 1- 22 


4) +4 1b22(g2—g4) +2a2bi(go—g,) | 


—I,f *Laibi(g2—ga) +aebo(1 — gs) +01b2(g2—g4) ta2b1(1 — 2g2+¢:) ] 


—1.f?([aibi(g2— gs) +aebe(1 — gs) +aib2(1 — 2g2+¢4) +a2bi(g2—g:) | 


O,=1/45kT 


+1,f°[a1b12(1 — 2g2+g4) +a2b2(1—g2+gs) +a1b2(—1+3g2—2¢4) +a2b1(-—14+3g2—2g; ] 


+31.f.fy[a@1b1(g2— gs) +d2b2(g2— 25) — dyb2(g2—g4) —d2b;(g2—g:,) | 


+3l,f-fy[aibi(ge— gs) +a2be(g2—2;) —d,bo(g2—g4) —a2b;(g2—g;) |. 


n is the index of refraction of the liquid, « its 
dielectric constant and N the number of mole- 
cules per cc. go and gy are again the functions of 
the potential V defined in Eq. (2). If there is 
no anisotropy of the Lorentz force x=0, 1,=1, 
=f,=f,=1 and Eq. (29) reduces to the classical 
result 6;=(a;—d2)(b; —b2)2kT/45. This is true 
independent of the values of ge and gy. Hence 
the hindered rotation alone does not alter the 
classical equation. 

For most nonpolar molecules the difference 
between molar polarization and refraction is 
small, hence a,;+2a,.=),+2b2. If we assume 
furthermore® p= po, then f,=/, and f,=/, and 
the Kerr constant can be written 

4  (m?+2)? (e€+2)? (b,—b2)! 
K,=-xN————_- ——- ———:x,_ (30) 
3 6n? v) 15kT 


where 


“Our expressions for dn, and dn, do not satisfy the 
relation dn,/dn, = —2, but the deviations are probably too 
small to be measurable. 

* This assumption may introduce a quite appreciable 
error in our results. 


(29) 


x= Ky, Ky= (1Fxitlixe 

+11-(l- +1.) x3 2 1/(p—1)?, (31) 
x1=1+2(p—1)g2+(p—1)*gs, 
x2= (p?— p+1) —(2p?—3p+1)g2+(p—1)%gs, (32) 
xs= —(p+1)+(2p?—5p+3)g2—2(p—1)* gs. 


Ky is the value of the Kerr constant calculated 
from the classical equation by using the observed 
values of N, m and « of the liquid and the polar- 
izabilities b;, b2 of the molecule as found from 
measurements in the gas. The theoretical value 
of x is a function of x and ge. This function is 
plotted in Fig. 8 for benzene, p=0.516. 

Fig. 8 shows that the Kerr constant depends 
on the structure of the liquid to a much higher 
degree than any of the other optical properties. 
If the potential barrier is small the Kerr constant 
rapidly increases*® with the anisotropy of the 
Lorentz force. 

The curves indicate that in liquids the Kerr 
~ 4 No data on compressed gases or liquids near the critical 
temperature seem to be available. Such measurements 


would be very valuable. They would furnish the proper 
interpretation of Rao’s results. See reference 37. 
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effect can be greatly reduced, eventually vanish 
and even reverse its sign. This marked structure 
sensitivity of the Kerr constant is confirmed by 
the experiments of Briegleb.” 

For }\=5460A McComb“ finds for benzene 
the Kerr-constant B= Kn/\=0.63 10-7; Briegleb 
measured B=0.408 10-7. Since all experimental 
errors tend to reduce B, we shall use the larger 
value. We find then K,=2.29 10-" and Ky 
=5.38 10-" from the experimental values of n, 
p, €, 01, be. Hence x = 0.426. 

Reference to Fig. 8 shows that this is only 
possible if again g2 >0.83. The Kerr effect checks 
with our previous result that Vo=8.8k7. Al- 
though in this case x appears somewhat smaller 
than previously this is not serious in view of our 
assumptions. 

A rise of temperature decreases V)/k7 and 
and hence increases the value of yx. Since 


V« 
Zo 


(1/K,)(dK,/dT) =(1/Ko)(dKo/dT) 
+(1/x)(dx/dTF) 


we conclude that the negative temperature 
coefficient of the Kerr constant should be smaller 
than that given by the classical theory. This is 
verified by the results of Hansen.** His data 
furnish (1/K,)(dK,/dT)=—1.6 10-°, while we 
calculate (1/Ko)(dKo/dT) = —4.4 10-*. A rough 
estimate shows that (1/x)(dx/dT) = +2.8 10°* is 
compatible with the theory.* 

If one assumes that (b;— 2) is proportional to 
(n?—1)/(m?+2) the classical theory leads to the 
Havelock relation Bnd/(n?—1)?=h, where fh is 
Havelock’s constant. According to our theory h 
must be replaced by hx. Hence the Havelock 
“constant’”’ should depend on the wave-length. 
However this variation is small. For benzene we 
estimate a change of x of about 5 percent 
between A4=4000A and 6000A. Within experi- 
mental errors of this order of magnitude the 
Havelock relation has been verified. 

The classical theory gives a relation between 
the depolarization of the scattered light and the 
Kerr constant. This relation of Gans cannot be 
verified in our theory, because there is no simple 
relation between M,, Mz and x. Actually for 


7H. E. McComb, Phys. Rev. 29, 525 (1909). 

8 From Landolt-Bérnstein, II. Erg. bd., Vol. 2. 

*°It is probable that in some liquids, especially near 
the critical temperature, the temperature coefficient of 
the Kerr effect should be positive. 
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Fic. 8. The Kerr constant for liquid benzene. The dotted 
line corresponds to the largest observed value. 


benzene the relation is only approximately 
correct. 
CONCLUSIONS 

The value of the molar refraction of liquid 
benzene, its dispersion, temperature—and _pres- 
sure dependence, the intensity and depolarization 
of the scattered light and its frequency-depend- 
ence, the Kerr constant and its temperature and 
frequency variation, all lead to essentially the 
same anisotropy of the Lorentz force and the 
same height of the potential barrier, when 
interpreted in the light of the present theory. 
It is probable that still better agreement could 
be reached by modifying some of our assump- 
tions. Since the potential barrier is very high 
we must expect that wave-mechanics would give 
somewhat different curves. A refinement of the 
theory does not seem worth while before more 
extensive and more accurate data are available. 

It appears then that hindered rotation and 
cybotactic grouping are sufficient to explain all 
deviations from the classical theory. The present 
theory holds only for nonpolar liquids with 
axial symmetry, but it can be generalized for 
more complex molecules. For polar molecules it 
must be modified. 

The writer acknowledges the help of Mr. Jakob 
Leeder in calculating the Lorentz factors for the 
hexagonal lattice, and is indebted to Dr. A. F. 
Turner for illuminating discussions concerning 
the problems of light scattering. 
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The experimental work of Dunning has shown that for 
fast neutrons the capture cross sections of heavy nuclei 
are small compared to their cross sections for elastic 
scattering. While it cannot be concluded from this that 
polarization effects (i.e., contributions to the elastic scatter- 
ing from quasi-stationary states of the nucleus-neutron 
system) are negligible, it seemed worth while to investigate 
the results of a treatment neglecting polarization. We 
use the effective potential for a neutron in the field of a 
heavy nucleus developed by Van Vleck. The values of the 
constants of the inter-particle interactions are adjusted 
slightly from those found by Feenberg from binding 
energy calculations. The problem reduces to that of a 
neutron in the field of a potential well whose depth is a 
slowly varying function of the neutron velocity. The scat- 
tering problem is solved by the method of partial cross 


INTRODUCTION 


XPERIMENTAL work on the scattering 

and absorption of fast neutrons has been 
carried out by several investigators.': 2: * Perhaps 
the most comprehensive investigation of the 
scattering of fast neutrons has been made by 
J. R. Dunning and his co-workers. Using as a 
source the neutrons emitted by the nuclear reac- 
tion between beryllium and alpha-particles from 
radon and its equilibrium products, Dunning has 
obtained values for the collision cross section of 
several heavy nuclei. He has shown further that 
even for nuclei known experimentally to capture 
fast neutrons, as evidenced by the production of 
gamma-rays or of artificial radioactivity, the 
capture cross section is of the order of less than 
ten percent of the cross section for elastic scatter- 
ing. The experimental collision cross sections fit 
a Z} law within the experimental error, Z being 
the atomic number. 

It is to be emphasized that we are here inter- 
ested in fast neutrons, of kinetic energy greater 
than one-half million electron volts. The experi- 
mental results for slow neutrons,*: * or neutrons 

1 Chadwick, Proc. Roy. Soc. A136, 702 (1932); A142, 1 
(1933). 

2 Curie-Joliot, J. de phys. et rad. 4, 21 (1933); 4, 278 
(1933). 

3 Dunning, Phys. Rev. 45, 586 (1934); Dunning, Pegram, 
Fink, and Mitchell, Phys. Rev. 48, 265 (1935). 

4 Amaldi, d’Agostino, Fermi, Pontecorvo, Rasetti and 
Segré, Proc. Roy. Soc. Al49, 522 (1935). 


sections. Anomalous resonances due to the attractive po- 
tential are found in the higher partial cross sections. Al- 
though these resonances have anomalously high peaks it is 
shown that on account of their sharpness they contribute 
little to the cross section for a neutron beam having a 
broad velocity distribution. The cross section is found to 
vary markedly with both velocity and atomic number. 
Even for a broad velocity distribution corresponding to 
Dunning’s experimental source oscillations remain in the 
cross section as a function of atomic number. However, the 
existing experimental points fit the curve within the experi- 
mental error. That the experimental results are in no case 
much larger than those given by the theory would perhaps 
indicate that contributions due to polarization effects 
are small. Possible experimental work to test the theory is 


mentioned. 


having approximately thermal velocities, are 
much more complicated. For these slow neutrons 
abnormally large collision cross sections are found 
for certain nuclei. In these cases the cross section 
for elastic scattering is negligible in comparison 
with that for neutron capture. 

Bohr® has shown since the inception of this 
work that these results can be explained only by 
consideration of the numerous and closely spaced 
energy levels of the nucleus. The chance of the 
neutron giving up some of its kinetic energy to 
other degrees of freedom of the many-particle 
system is very good; robbed of most of its veloc- 
ity, the neutron stays long enough in the nucleus 
to make the probability of capture large. Breit 
and Wigner® have independently dealt with the 
considerations discussed by Bohr in more explicit 
mathematical detail. They consider the effect of a 
single quasi-stationary level of the combined 
system of neutron and nucleus which has an 
energy in the thermal region, and show that it is 
reasonable that the breadth of the level due to 
radiation damping is large in comparison to that 
due to re-emission of a neutron. The ratio of the 
absorption to the extra scattering due to the level 
is correspondingly large. Of course the effect is 
only important near resonance, i.e., when the 


5 N. Bohr, Nature 137, 344 (1936). 
° Breit and Wigner, Phys. Rev. 49, 519 (1936). 
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energy of the system is near that of the quasi- 
stationary level. 

According to Bohr the nucleus plus neutron 
system has such quasi-stationary levels, more 
closely spaced than those at thermal energies, in 
the energy range corresponding to fast neutrons 
Here the absorption due to such levels may be 
small in comparison to the scattering, and iti s 
impossible to conclude from the experimental 
absence of large absorption that the extra scatter- 
ing due to the existence of these levels, or, in 
other language, to polarization effects, may be 
neglected 

In the light of the above remarks the simple 
treatment of the problem dealt with in this paper 
is not to be taken too seriously. The work was 
well begun, however, before Bohr’s paper made 
clear the possibility of large polarization effects, 
and it was thought of some interest to complete it. 
At the present stage of nuclear theory a more 
refined treatment appears impractical. It should 
be said that the vices of our treatment are in- 
herent in the use of the statistical model. Its re- 
sults may be worse than those for binding energy 
calculations with the statistical model, but only 
because scattering calculations are generally 
more sensitive to approximations than energy 
calculations. 


EXTENSION OF VAN VLECK’S TREATMENT OF THE 
INTERACTION BETWEEN A NEUTRON 
AND A HEAVY NUCLEUS 

J. H. Van Vleck’ has considered the interaction 
of a neutron and a heavy nucleus on the basis of 
a statistical treatment of the particles of the 
nucleus. He employs the types of inter-particle 
interaction found necessary to the theory of the 
binding energy of nuclei. Binding energy con- 
siderations lead to neutron-proton interaction of 
the Majorana* (exchange) type. The Wigner? 
(ordinary) potential makes the binding energy of 
heavy nuclei increase too rapidly with atomic 
weight. Similarly neutron-neutron interaction of 
the spin coupling type’ is required; in other 
language than that of the vector model this cor- 
responds to one part Wigner to two parts Majo- 
rana interaction.'° 


* Van Vleck, Phys. Rev. 48, 367 (1935). 

5’ E. Majorana, Zeits. f. Physik 82, 137 (1933). 

* Wigner, Phys. Rev. 43, 252 (1933). 

10 Feenberg and Knipp, Phys. Rev. 48, 906 (1935). 
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The complete Hamiltonian is thus of the form 


IT = Hy) — (h?/8x?M)y,? 
Z \ 
tT (rap) P™(n, p)NtFEN (ran) Sn°Sn, (1) 
j=l 

where //) is the nuclear Hamiltonian, J/ is the 
mass of the neutron (or more properly the re- 
duced mass of neutron and nucleus if the origin 
of coordinates is chosen as the center of gravity 
of the nucleus), J(7,»;)P™(n, p;) is the Majorana 
interaction between the neutron and the jth 
nuclear proton, and N(r,»,)Sa°S,,; is the inter- 
action between the neutron and the ith nuclear 
neutron. P™(n, p;) is the Majorana exchange 
operator permuting the coordinates of the neu- 
tron and the jth proton, s,, is the vector spin 
matrix for the ith nuclear neutron, Z is the 
atomic number, and N the number of neutrons 
in the nucleus. //) of course contains terms similar 
to the last two in Eq. (1), so that the Hamil- 
tonian is symmetric in the neutrons. 

We shall neglect nuclear neutrons not in closed 
shells, so that the last term of Eq. (1) will not 
contribute directly to a_ self-consistent field 
effective potential. However, when exchange is 
taken into account it does enter. From the Dirac 
vector model,'' which shows that the exchange 
effect for identical particles a and bd is formally 
equivalent to insertion of a spin coupling factor 
—(3)(1+4s,-s,), the exchange potential due to 
neutron exchange is here 


— 3 (1/2)P™(n, ni)(1+48n° Sn.) N(ran,) Sn* Sn.» 
or 
— > (3/8—(1/2)8n- Sn) N (ran) P™(n, ni), 


in virtue of the matrix identity 16(s,-s,)* 
+8s,-s,—3=0. Thus the Hartree-Fock equation 
satisied by an approximate one-particle wave 
function for the neutron is 


Z 
—(h? sruveva.)+E | f f eir,) 
j=1 


XJ (rap) ¢)(Xn)¥(X,)dr>p 


N f° 
—(3/8)>- [ | Vi* (Xn) N (Tan, 
IY v7 wo 


Xvi(xX,)¥ (Xn, dtr, = WY(x,), 
1 Cf. Van Vleck, Phys. Rev. 45, 405 (1934). 


(2) 
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where g(x), ¥;(X) are one-particle wave functions 
for the jth proton and ith nuclear neutron re- 
spectively, and W is the kinetic energy of the 
free neutron. 
We need an approximate expression for 

Z 

> $j" (Xp) ej(X , 

y= 1 
the Dirac density matrix” for the protons in 
the nucleus. This may be obtained by suppos- 
ing that phase space is as full as permitted by 
the Pauli principle within the volume r< R, p< P, 
and empty elsewhere. Necessarily (47/3)*P*R® 
=Zh*/2. We assume further that the proton 
wave functions are adequately represented locally 
by plane waves, that is, that the local momentum 
does not vary too rapidly with position within 
the nucleus, and that all momenta p<P are 
equally probable. Then 


Z a) 
) ¢i* (Xp) ¢j(Xn) = (2 mf f 
2=1 P 


exp [(27i/h)(x»—Xn)- p dp, 
os |xp+Xn|/2<R; 
— X,+x,|/2>R, 

the integration being over a sphere of radius P. 


This gives 


Zz 
> ¢)*(Xp) ¢j(Xn) = (sin Rpt ap 
?7=1 
—kofnp COS Rypfnp)/ Pap’, 
| (3) 
Xpt+Xn|/2<R; 
=0, 
Xp>+x,|/2>R, 
where ky =2rP/h=(9nZ/4R*)!. 


A similar expression holds for 
N 
De Vi* (Xna)Wi(Xn), 
i=1 


the Dirac density matrix for the nuclear neutrons, 
with k, replaced by k, =(9rN/4R*)!. It is con- 
venient to assume R the same for both expres- 
sions; since N is roughly equal to Z and k,, in- 
volves N in a low power a further simplification 
may be attained by replacing k, by kp. With 





# Dirac, Proc. Camb. Phil. Soc. 26, 376 (1930). 
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these simplifications our result is what we would 
have obtained had we considered only neutron- 
proton interaction with a J(r) equal to J(r) 
—3N(r)/8; Eq. (2) becomes 


Vn'W (Xn) +R (xn) = (82? M/h?) SSS [J (r’) 
— (3/8)N(r’) JL (sin kyr’ —kpr’ cos kyr’) / x?” 


Xv(x,+x’)dr’, ro<R; =0, r,>R, (4) 


where we have written k?=82°WW,h*. In ob- 
taining Eq. (4) we have replaced the conditions 
on Eq. (3) by the condition that the Dirac 
density is zero when the neutron is outside the 
nucleus and has the given value inside, a proce- 
dure justified by the small range of the neutron- 
proton and neutron-neutron interactions in com- 
parison with R. Consideration of edge effects 
would prove a useless refinement in view of the 
fact that we have neglected them in obtaining 
Eq. (3). 

We proceed to the problem of solving Eq. (4) 
for r,<R. It is of the type 


V2u(x) +k?u(x) = SSS fir’ )u(x+x')dr’. 


Now this apparently troublesome integro-differ- 
ential equation may be reduced to an ordinary 
homogeneous differential equation, as Feenberg!* 
has very kindly shown. Namely, if we try a plane 
wave solution, “=exp 7k’ -x, we find it a solution 


provided 
K+ R= SSS I(r’) exp ike -x’dr’. 


Now on account of the spherical symmetry of 
f(r’) the integral is independent of the direction of 
the vector k’. Since any solution of V?u+k”"u=0 
can be constructed from plane wave solutions 
(having of course the same magnitude of k’), this 
solution is also a solution of the integro-differen- 
tial equation. Thus we may solve Eq. (4) by 
solving 


VV+R"Y=0, r<R; ; 


(9) 
V¥ytky=0, r>R, 
with 
hk” =k? —(82?M/h?) SSS (J(r) —3N(r)/8) 
X[(sin kyr —kpr cos kyr) /x?r* ] (6) 


Xexp tk’ -xdr. 
‘8 Private communication to J. H. Van Vleck. Van Vleck 
had previously obtained the result for an s wave function by 
explicit calculation, which we extended to » and d func- 
tions. 
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The problem is thus reduced to that of a neutron 
in the field of a potential well, the depth of the 
well being a function of the velocity of the neu- 
tron on account of the appearance of k’ in the 
integral in Eq. (6). 

Feenberg'® has succeeded in fitting the ob- 
served binding energies of light nuclei by taking 
Jir) and N(r) of the functional form, 
J(r)=A exp (—ar*), N(r)=B exp (—ar’). Inte- 


same 
gration of Eq. (6) with these gives 
k” =k? — (82°M/h*?)(A —3B/8) E(w.) —E(w_) 
+ (4a/rk”) (exp (—w,*)—exp (—w_*))], | 
(/) 
ws = (k’+k,) 2a}, : 


at 
E(w) = (2/1) | exp (—x*)dx. 
e 


0 


Calculation with (7) shows that the potential well 
becomes gradually shallower as the neutron 
velocity is increased, as would be expected from 
the nature of exchange forces. 

Before considering the scattering problem it 
will be well to examine the numerical magnitudes 
involved. Feenberg'® gives as values of A, B, and 
a best fitting the binding energies of light nuclei 


A=-—74mce, 3B /4=26me’, 8) 
1 7 ‘ ? 
a~?=2.17-10-'* cm. 


We assume the Gamow value R=7.8-107-'* cm 
for Pb and that R is proportional to Z!. The latter 
assumption has an extremely important effect; 
namely, since Z and R enter the expression for k’ 
only ink, = (9%Z /4R*)' thedepth of the well is inde- 
pendent of Z, in line with the fact that the bind- 
ing energy per particle is approximately constant 
for heavy nuclei. (We neglect the slight variation 
of reduced mass with Z.) Now from M=1.670 
-10-*4 g we find k=0.220-10" (W(MEV))! cm™. 
Also ky =1.069-10'* cm~'. With these values k’ 
has the value 1.060-10'* cm~ for k=0. We see 
that even for W as high as 7 MEV (million elec- 
tron volts) & is considerably less than k’, that is, 
the depth of the potential well is very large (of 
the order of 24 MEV) compared to the kinetic 
energy of the neutron. Since 7 MEV is about the 
upper limit of kinetic energy for the fast neu- 
trons from a beryllium-radon source, these neu- 
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trons are not fast in the sense in which one 
speaks of fast electrons in discussing their scat- 
tering by atoms. The Born approximation," so 
useful in many scattering problems, is inappli- 
cable here. The problem may be solved in terms 
of known functions by the method of partial 
cross sections, which though tedious in its com- 
putational application leads in our case to reson- 
ances no simpler treatment could well be expected 


to fit. 


THE METHOD OF PARTIAL CROSS SECTIONS 


The method of partial cross sections'® depends 
on the fact that the wave equation for the motion 
of a particle in a central field is separable in polar 
coordinates, possessing particular solutions which 
are products of functions of the radius and spher- 
ical harmonics, in which more complicated solu- 
tions may conveniently be expanded. We wish to 
express a function representing an incoming 
monochromatic plane wave plus an outgoing 
radial wave, a reasonable idealization of the 
usual experimental arrangement for the study of 
scattering, in terms of these particular solutions. 
On account of the axial symmetry we shall need 
only the particular solutions of the form f;(r)P; 


(cos 6). Our wave equation is of the form 
Vyt+(k?— U(r))y=0. (9) 


If U(r) vanishes not less rapidly than 1/7? for 
large r and behaves no worse than 1/r at the 
origin f;(r) satisfying the boundary condition at 
the origin exists and has the asymptotic form 
(1/kr) sin (kr—In/2+¢,) for large r. The plane 
wave exp (tks) =exp (kr cos @) is an asymptotic 
solution of Eq. (9) and has the asymptotic 


expansion 


exp (tkr cos 6) 
=(1/kr)>>(i)'(2/+1) sin (kr—Ir, 2)Pi(cos 8) 
l=0 


for large r. We write therefore 


4 See, e.g., Mott and Massey, The Theory of Atomic 
Collisions, Chapter VII. 
1 Cf. Mott and Massey, reference 14, Chapter II. 
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exp (7kr cos @)+/(8) exp (tkr) r 


=(1 kr)>-(2)'(21+1) sin (kr—Ir/ 2) Pi (cos 6) 


+ 3° A ,(2/4+1)Pi(cos 0) exp (ikr), kr (10) 


=(1, kr)>— B,(2/+1) sin (kr—lr/ 2+ ¢,)Pi(cos @). 


Equating the coefficients of exp (ikr) and of 


exp (—ikr) gives 
A,=exp (2i¢;) —1, 


so that 


f(0)=(1 k) (2141 (exp (21¢,)—1)Pi(cos 6). (11) 


From (11) we have 


1(6)= | f(@)|?=(1 k)> > (2/+1)(2l’+1) 
l=0 l’=0 (12) 


X sin Li sin gy Cos (¢,;— gv) P (cos 6)P (cos 8), 


the scattered current per unit solid angle when the 
incoming wave is normalized to unit current 
density. It is seen to be determined by the ¢;'s, 
the so-called phase shifts of the radial functions. 

The total cross section Q is obtained by inte- 
grating /(@) over all directions, and is 


QO=>> Q.:=(44 k?)>° (21+1) sin? g. (13) 
i=0 l=0 


k' S14 i(R’R)J, 
h! J" 14 4(R'R)J—1-\(RR) —RJ’_1_\(RR) Jia (RR) 


tan ¢;= —(-1)! 


H. FAY 


Physically Q; may be regarded as the partial 
cross section due to particles of angular momen- 
tum about (/+43)h/27. On this basis one would 
expect Q, to be important only when the classical 
distance of closest approach for particles of the 
corresponding angular momentum is such that 
the particle classically traverses regions where 
U(r) is appreciable. This is not quite true, how- 
ever; in the case of an attractive potential, such 
as we must consider, there may be three classical 
turning points (radii), with virtual states similar 
to those met with in the theory of radioactive 
decay and corresponding sharp high resonances 
which we may refer to as anomalous for want of a 
better term. These anomalous resonances compli- 
cate our problem, but on account of their sharp- 
ness contribute little to the cross section for a 
neutron beam having a broad distribution of 


velocities. 
SOLUTION OF PROBLEM OF SCATTERING BY A 
POTENTIAL WELL 


We now proceed to obtain an expression for 
the phase shifts of the radial solutions of (5). 
Inside the well the /th radial function is 


filr) =e Ji, (Rr); 
outside we have 
filr) =cer Leos gi Ji. ;(kr) 
+(—1)! sin g:J_1-;(kr) }. 
From the continuity of f;(7) and of its first deriva- 


tive at r=R we find 


; (RR) —RI' 14 (RR) Jin (R'R) 
(14) 


This may be reduced to a more convenient form by eliminating the derivatives by means of the 
recurrence formulas satisfied by the Bessel functions, with the final result 


tan gi= 


a formula convenient for computation. From 
Eqs. (13) and (15) it is evident that here Q/R? is 
a function of but two parameters kR and k’R. 
We have computed tables for the evaluation of 
phase shifts and Q/R? for a range of values of 


Ji-\(RR)Jia3 o(R’R) —Jiss/2(RR) Ji i(k’ R) 
(—1)'[J_145(RR)Jiasyo(k’R) — J—1-3/2(RR)Jiv_y(k’R) ] 


(15) 


these parameters sufficient to cover the values 
met with in our problem. These tables will be 
loaned to anyone who desires to use them. 

We wish to make an estimate of the area under 


an anomalous resonance peak of a Q, R® versus 
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kR curve (we do not need to consider the manner 
of variation of k’R with kR, assuming it constant 
over the small range involved); this we may do 
by using the first term of the power series expan- 
sions for Ji,;(RR) and J_)-,(RR) in Eq. (15), 
since RR is small in comparison to / for anomalous 
resonance. It can be shown that 


tan og, = (f9)2!+"/[ (274+-1)(1-3-5- + -(21—3))? 


x (Cf)? —(RR)*) ], (16) 


where ft is the value of RR for resonance. Thus 
for tan g;= +1, the points at which the curve has 
half its peak value, 


(RR)? = (FR) (FM)? /C(21+-1) 
X (1-3-5-+-(2—3))*], 


and the peak breadth at half value is 


A(RR) = (FM)? [ (2 +1) (1-3-5- + -(2l—3))*], 
j=1. 


Since the curve is like a sin? curve the area under 
it is approximately this half value breadth times 
the peak height, or 


(Q1/R?) maxA(RR) 


=4r(fR)?"*/[1-3-5---(27-—3) #, (17) 


a result which should be reasonably accurate for 
ft small in comparison with /. We see that in 
spite of the increase of peak height with increas- 
ing anomaly (decreasing fi) the breadth de- 
creases more rapidly, so that the area under the 
peak, that is, the effective contribution to the 
cross section for a broad distribution of velocities, 
decreases. It is evident that experimental detec- 
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tion of these anomalous resonances will require a 
source of neutrons of quite homogeneous velocity. 


RESULTS 


It was found that the best fit with the ex- 

perimental data was obtained by increasing 
A—3B/8| of Eq. (8) by five percent; a suitable 

adjustment in R would have served about as well. 
k’ is raised from 1.060-10'* cm to 1.074-10' 
cm! for k=0, or less than two percent, by this 
change. 

Behavior of cross section with velocity. Fig. 1 
illustrates the variation of cross section with 
neutron velocity for several elements. Although 
the constants are such for these elements that 
no highly anomalous peaks appear, the curves do 
show important resonances. The cross section 
varies so much with velocity that the effect ought 
to be experimentally detectable by the use of 
some other source than beryllium-radon, having 
a sufficiently different velocity distribution to 
lead to a measurably different cross section. It is 
hoped that such work will soon be undertaken by 
those having the necessary facilities. Of course a 
homogeneous source would be ideal, but its 
experimental realization involves difficulties not 
readily surmounted. It is possible that the neu- 
tron beam from beryllium bombarded with 
deuterons or from deuterium bombarded with 
deuterons may prove sufficiently homogeneous to 
allow the detection of narrow resonances. 

Variation of cross section with atomic number. 
Fig. 2 indicates the behavior of the cross section 
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as a function of atomic number for neutron ener- 
gies of 1, 3, and 6 MEV. Here of course anomalous 
resonances appear. The partial cross sections 
giving rise to the more prominent resonances are 
indicated, as are the heights of those peaks lying 
off the scale. 

The most promising sources of homogeneous 
neutron beams give neutrons of kinetic energy 
around 2 MEV. Fig. 3 shows cross section as a 
function of atomic number for 1.86 MEV and 
for 2.39 MEV neutrons. It indicates what is to be 
expected on the basis of this theory for experi- 
mentally obtainable sources of homogeneous 
neutron beams, such as deuterium bombarded by 
deuterons. 

It is evident from Figs. 2 and 3 that the cross 
section does not even begin to approximate a Z! 
law for any single value of neutron velocity. It 
would be desirable to attempt to locate some of 
the more anomalous peaks with a highly homo- 
geneous neutron beam, as from their relative 
locations one could infer much concerning the 
importance of edge effects. Consideration of these 
in the theory would not change the nature of the 
resonances, but would alter their relative posi- 
tions. 

Results for Dunning’s experimental velocity 
distribution. Dunning, by measuring the ranges 
of protons projected by neutrons from his source, 
has given us a fair idea of the velocity distribution 
of berylliium-radon neutrons. He suggests'® a 


16 | am indebted to Professor Dunning for an interesting 
conversation on this point. 
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distribution of the general nature shown in 
Fig. 4. With this distribution we find a behavior 
of cross section with atomic number shown by 
Fig. 5. It is seen that the broader resonances indi- 
cated in Figs. 2 and 3 persist in the averaging 
process. The peak around Z = 34 can be identified 
with resonance in the first partial cross section. 
The low at Z=45 is near the vanishing of the 
first partial cross section. Resonances in the first 
and second partial cross sections give the peak 
at Z=60. On the high velocity end the fourth 
partial cross section contributes somewhat to the 
peak at Z = 34 and the fifth to the peak at Z= 60 
The crosses on Fig. 5 indicate Dunning’s ex- 
perimental! results for Cu, Zn, Sn, I, W, Hg, and 
Pb, the only elements in the range given for 
which he has made measurements with fast 
neutrons. Except for Cu, the agreement is within 
the possible experimental error, which he esti- 
mates to be about ten percent. However, his 
relative values are perhaps better than this, so 
that we are disappointed not to find the decline 
from Hg to Pb indicated by his results. As may 
be seen from the 6 MEV curve of Fig. 2, however, 
the rise we get comes from the high energy end of 
the distribution, and changes in the assumed 
distribution curve might reduce or eliminate the 
discrepancy. We cannot expect good agreement 
at low Z, so that the bad agreement for Cu should 
not be taken too seriously. It is unfortunate that 
Dunning has measured cross sections for slow 
neutrons for many elements lying at the interest- 
ing portions of our curve without investigating 
their cross sections for fast neutrons. 
Comparison with constant depth well. It is of 
some interest to inquire how the scattering by 
the well of depth a function of the velocity com- 
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pares with a constant depth potential well. The 
essential difference between the two lies in the 
fact that k’ varies more slowly with & for our 
well than for the constant depth well. This has 
the effect of making the cross section vary more 
slowly with velocity; resonances are broader for 
the well of varying depth. For this reason the 
results given by the constant depth well for our 
velocity distribution will show less oscillation 
with atomic number than our curve. At present 
there is no experimental reason for choosing one 
as preferable to the other, as the experimental 
points can be fitted equally well with the constant 
depth hole. 

Importance of small angle scattering. The experi- 
mental work of Dunning does not measure scat- 
tering through angles of less than about 7°. The 
question suggested itself that possibly the theory 
might give such a large small angle scattering, 
particularly near resonances in the higher partial 
cross sections, that it would be necessary to take 


this into account in making comparisons with 
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the experimental results. Calculation in a special 
case revealed that while the small angle scatter- 
ing was of the order of four times the average, 
the contribution to the total cross section could 
be neglected on account of the small solid angle 
of the 7° cone, about 0.3 percent of the total solid 
angle. 

Nature of polarization corrections. The experi- 
mental cross sections for fast neutron scattering 
are in no case much larger than those given by 
our theory. We may be justified in concluding 
from this that polarization effects are small. 
Contributions to the cross section from polariza- 
tion may be expected simply to add to the total 
cross section, although it is possible that inter- 
ference (nonorthogonality) is important. 

It is possible to see in a general way the effect 
of quasi-stationary levels of the nucleus-neutron 
system. We should expect sharp resonances in 
the cross section corresponding to these levels, 
so that the behavior of the cross section as a 
function of velocity or as a function of atomic 
number for a homogeneous beam would be con- 
siderably more complicated. On the other hand, 
it is possible that the cross section as a function 
of atomic number for a broad velocity distribu- 
tion might vary more smoothly than our result 
indicates, on account of the great number of 
resonances to be averaged over. In any event we 
believe that further experimental work on fast 
neutron scattering, a field somewhat neglected 
since the discovery of the exciting properties of 
slow neutrons, should yield results both interest- 
ing and important. 

In conclusion I wish to thank Professor Van 
Vleck for his helpful interest in the problem. 
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Note on the Analysis of Photoelectric Data 
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(Received July 6, 1936) 


the development of Fowler’s' method for 
the determination of photoelectric work 
functions, the assumption is made that ‘“‘the 
photoelectric current J per quantum of light 
absorbed is proportional to N,’’ where Nz is the 
number of available electrons. It has been 
customary, however, to plot current per unit 
intensity of incident radiation assuming a con- 
stant optical reflection coefficient; also the term 
(W.,,—hv)~ has been assumed constant. 

We have made comparisons of different 
methods of plotting, using data on Ba taken in 
this laboratory for the purpose of determining 
the temperature coefficient of the work function. 
Assuming the constancy of the terms mentioned 
above, a plot of current per unit intensity is 
shown in Fig. 1, curve A; curve B shows elec- 
trons/incident quantum plotted as ordinate. 
In curve C the quantity (W,—hv)'!Xelectrons 
per incident quantum is plotted as ordinate. The 
optical reflection coefficient varies so little over 
the range of frequencies used that one cannot 
distinguish between the fits to the theoretical 
curve in plots of electrons per incident quantum 
and electrons per absorbed quantum. It would 
appear from the figure that the experimental 
data fit the theoretical curve about equally well 
in the three cases. However, careful examination 
of large scale plots for different sets of data indi- 
cates that the least accurate fit is obtained for 
the plots of electrons per quantum. 

According to Fowler’s theory, curve C would 
yield the true work function. The restriction that 
W. be known is not a serious one, since a change 
as large as 25 percent affects appreciably neither 
the fit nor the value of ¢ obtained. The value of 
W, used for curve C was computed on the basis 
of two free electrons per atom. Waterman and 
Henshaw’s’ extension to Fowler’s analysis was 
applied to the problem. Plotting electrons per 
quantum, the fit to their theoretical curve was as 
good as for curve C and the work function so ob- 
tained was the same as that found from curve C. 


1R. H. Fowler, Phys. Rev. 38, 45 (1931). 
2 Waterman and Henshaw, Phys. Rev. 44, 59 (1933). 


The values of the work functions obtained 
from the three curves given would necessarily be 
slightly different. In spite of the excellent fit to 
curve A, the accurate determination of work 
functions from such plots is not justified by ex- 
isting theories. Curve B has, as a theoretical 
basis, DuBridge’s* modification of Fowler’s 
theory since his development does not yield the 
factor (W,—hyv)~!. Should further work show 
conclusively that electron-per-quantum plots do 
not fit Fowler’s curve accurately, the deviations 
might well be absorbed by the proportionality 
constant a, which DuBridge assumed to be 
essentially constant. 

DuBridge has interpreted the proportionality 
factor @ as ‘the fraction of the number of elec- 
trons which come up to unit area of surface in one 
second which absorb a quantum of energy when 
the incident light intensity is unity,’ or ‘‘the 
probability that an electron at the surface of the 
metal will absorb a quantum and appear outside 
as a photoelectron.’’* The experimental value of 
«a for Ba based on this interpretation is 4.3 K 10-* 
cm?’ sec. quantum ~'. Simple theoretical considera- 
tions’ show that the optimum value of a, as- 
suming two electrons per atom is 1.4 10~*°. The 
ratio of the a’s 1/325, equals the ratio of experi- 
mental to maximum current possible according 
to this theory. 

Fowler attached no significance to the constant 
of proportionality appearing in the expression 
for the current; however, it may be interpreted 
in a manner analogous to the interpretation of a. 
Writing Fowler’s equation as 

constant 

I= T?o(x), 

(W.—hv)! 
where J is in electrons per quantum, 8 may be 
defined as the fraction of the number of available 
electrons per unit volume which will be emitted 
due to one incident quantum, or the probability 
that an available electron in unit volume of the 


‘L.A. DuBridge, Actualités Scientifiques et Industrielles, 
No. 268 (Paris, 1935). 

*Marvin M. Mann and Lee A. DuBridge, Bull. Am. 
Phys. Soc. 11, No. 4, Paper 24. 
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Fic. 1. Curve A: ordinates, current/unit intensity; horizontal shift, 98.58; ¢=2.508 volts. 
Curve B: ordinates, electrons/quantum; horizontal shift 98.85; ¢=2.515 v. Curve C: ordinates 
(W.—hv)? Xelectrons/quantum; horizontal shift, 98.65; ¢ =2.519 v. Full lines, Fowler's theoretical 


curve. 


metal will absorb a quantum and appear outside 
as a photoelectron. The experimental value of 
the constant 8 obtained for Ba was 5.010-* 
cm*® quantum". The value of 8 computed on the 
basis of Waterman and Henshaw’s analysis 
differs inappreciably. 

A theoretical limiting value of 8 may be 
found by assuming that one electron is emitted 
per quantum and that the number of available 
electrons is equal to 1/2 the total number m of 
free electrons. The latter assumption is made 
because, according to Fowler’s theory, all the 
energy of the absorbed quantum is added to the 
kinetic energy normal to the surface and but 


half the electrons have outward normal com- 
ponents. Hence 


1 electron, quantum= Ben_ 2. 


Assuming two electrons per atom 8 becomes 
6.0 10-* cm* quantum™'. The ratio of the 6's 
1/120, is the ratio of experimental to theoretical 
maximum current. Since the probability 8 (or a) 
depends on the probability of absorption, the 
probability of electron transmission and on the 
number of available electrons, at least a part of 
this discrepancy may be attributed to a lack of 
total absorption of light in that region in which 
photoelectrons originate. 
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A Method for the Determination of Absorption Cross 
Sections for Thermal Neutrons 


The determination of the coefficient of back scattering of 
slow neutrons (absorbed by cadmium) from a thick layer 
enables one to calculate the ratio of the scattering and 
absorption cross sections, o, and og, of the back scatterer. 
As the total cross section (e =0;+0,) can be measured by 
direct methods, the back scattering cross section can be 
used to determine the absorption cross section. If the 
fraction of neutrons scattered backwards is denoted by S, 


the following relation can be derived:' 


2a,.+c, 2a,.+0e;\” 4 
s -(( )-1). 
oC. o 


A convenient way of measuring this back scattering 
consists in immersing a detector foil in a liquid and com- 
paring the activity caused by the neutrons for the bare 
foil with that when the neutrons have been screened on one 
side, in each case subtracting the activity when the foil 
is screened on both sides by cadmium. The ratio R of these 
two activities is given by R=2/(1—S), from which follows 
the simple equation: (R—1)?=¢/o,. For large values of R 
this result becomes identical with an expression given by 
Amaldi and Fermi,’ but in most cases it is not permissible 
to neglect the term —1 with respect to R. 

In actual experiment a correction must be made for the 
absorption in the foil. If a fraction « is absorbed, the 
measured ratio is not R, but R’, which is given by 


2 R’(1—«) 
, Whence R= — 


" f—f/aek” 


, 


mes er” 


The cross sections occurring in these formulae are the 
weighted mean cross sections of the various atomic con- 
stituents in the scatterer. By varying the chemical com- 
position of the scatterer, and especially by studying solu- 
tions and liquid or powder mixtures of various relative 
concentrations, it should be possible to obtain information 
concerning the cross sections of the different atoms. In 
fact, this method of mixtures, in principle, offers the 
possibility of determining the ratios of any two of the 
cross sections for scattering, capture, or collision of any of 
the individual atoms in the mixture. Hence, data obtained 
in this way, taken together with the accurate directly 
determined value of just one particular cross section, 


‘ 


should permit the calculation of the absolute values of all 
the individual cross sections involved. 

For example, if R is measured for a very dilute solution 
of a cadmium salt, the total cross section will hardly be 
different from that of water, as it is primarily due to the 
large hydrogen scattering. A change in R will then be due 
to the change in the mean capture cross section caused by 
the presence of the cadmium. We have performed some 
preliminary experiments with several liquids and solutions. 
In the case of a cadmium solution an effect of the correct 
order of magnitude was observed. Since these measure- 
ments have so far been made with a very weak source, we 
have not yet obtained accurate quantitative results, but the 
method promises to give results of considerably higher 
accuracy than that thus far obtained. 

The experimental arrangement was as follows: A radon- 
beryllium source was placed at the center of a large cylin- 
drical tank of water. A glass cylinder of 8 cm diameter was 
placed with its axis about 8 cm from the source and filled 
with the liquid or solution to be investigated. The silver 
detector foil was placed in a plane passing through the 
source, in order to have a symmetrical arrangement as 
regards irradiation of the foil from both sides. 

In these preliminary experiments we used as detectors 
silver foils of 0.148 g/cm2, for which « should be equal to or 
less than 0.10 as based on the calculations of Amaldi and 
Fermi. The foil was activated for two minutes. The value of 
R obtained in water was 7.3+9.4, which is considerably 
less than the value R=11, given by Amaldi and Fermi. 
They used a thinner foil, placed in a plane perpendicular to 
a radial line from the source and activated it only one 
minute. We have not yet found the cause of this dis- 
crepancy. 

This work is being carried out with support from the 
Rackham Fund of the Graduate School. We are also 
indebted to the Saskatchewan Cancer Commission for the 
radon used in these experiments, and especially to Mr. G. 
Wrenshall who prepared the sources. 

C. T. ZAHN 
E. L. HARRINGTON 
S. GouDsMIT 

Department of Physics, 

University of Michigan, 

August 21, 1936. 


1D. L. Bayley, B. R. Curtis, E. R. Gaerttner and S. Goudsmit, 


Phys. Rev. 50, 461 (1936). 


2? E. Amaldi and E. Fermi. Ricerca Scient. 7, 3 (1936). 
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LETTERS 


The Scattering of Gamma-Rays 


The experiments described by the writer in these pages 
some time ago' have been extended. The apparatus has been 
modified so that a single gamma-ray counter and a single 
electron counter are used. Figs. 1 and 2 show the relation 
of the scatterer and the counting tubes used in the experi- 
ments. The gamma-ray counter G has a copper wall of 
sufficient thickness to absorb the recoil electrons coming 
from the paraffin scatterer. The electron counter R has 
the recoil electrons to enter 


very thin walls which allow 


easily. 








¥ ri pee Ke 
\ 
.@ ~ 
OR 
Fic. 1. Position of gamma- Fic. 2. Position of gamma- 
ray and recoil electron ray and recoil electron 


counters for cases A and B. counters for cases C and D 

The thickness of the scatterer (paraffin 410.05 cm) 
was such that the recoil electrons in emerging from any- 
where save near its edge would suffer multiple scattering. 
Thus the angular positions of the counters were not 
especially critical. The inhomogeneity of the gamma-ray 
source used (84 mc of Ra C) would also tend to make the 
exact positions of the counters uncritical. 

Observations were made with the counters in the two 
positions shown. In each case readings were taken both 
with the paraffin in place and removed. The experiment was 
carried on in a vacuum drum to eliminate scattering by the 
air. The results were as follows: 


A Position 1 with scatterer in place Rate =69.1+5.4 hr. 
B Position 1 with scatterer removed Rate = 10.7+4.4 hr.~ 
C Position 2 with scatterer in place Rate =10.1+5.0 hr. 
D Position 2 with scatterer removed Rate = 12.5+5.6 hr.~! 


In cases C and D where both counters are on the same 
side of the scatterer no coincidences would be expected 
except between gamma-rays and associated recoil electrons 
which had suffered a deflection of approximately 180°; but 
the thickness of the scatterer is such that the great majority 


of these must be absorbed. The observations show that 
within the experimental uncertainty no increase in the 
coincidence rate results from the introduction of the 


scatterer. 

The results A and B obtained with the arrangement 
shown in Fig. 1, however, show a decided increase in the 
recorded coincidences when the paraffin is introduced, 
indicating that there is a coincidence in time between the 
appearance of a recoil electron and the scattered gamma- 
ray which liberated it. 

The residual coincidence rate of about 10 hr.~' is due to 
the finite resolving time of the circuit and to cosmic-ray 
showers. 


RoBERT S. SHANKLAND 
Case School of Applied Science, 
Cleveland, Ohio, 
August 24, 1936. 


'R. S. Shankland Phys. Rev 


49, 8 (1936) 
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Behavior of Slow Neutrons at Different Temperatures 


In the study of slow neutrons the reflection coefficient 
of hydrogenated substances for thermal neutrons plays an 
important part because it gives information on the average 
number of collisions N occurring before the neutron is 
“‘albedo”’ 


the 


captured. The reflection coefficient, called the 
by Amaldi 


assumption of isotropic scattering of the neutrons at every 


and Fermi,' can be calculated? under 


collision and is given by 


B=(VN-1)/(/N+1) =1-2/VN+::: 1) 


N is of the order of magnitude of 109). The quantity NV 


is connected with the mean free path \ of the neutrons in 


paraffin, their mean life + in paraffin and their mean 


velocity v, by the equation 
N=vr/n. 


Assuming r and \ to be independent of the velocity, as 
one would expect for neutron energies very small compared 
with the vibrational energy Av of the hydrogen atom in 
paraffin, the albedo should decrease with temperature 
because v varies as \/T. 

We performed some experiments to measure @ in paraffin 
at liquid-air temperature (90°K). The experimental pro- 
cedure was similar to that used by Amaldi and Fermi.! 
We used two cylindrical paraffin blocks 12 cm in diameter 
and 10 cm in height. The source was at a distance of 2.5 cm 
from the detector and was kept at constant temperature 
by a small Dewar flask. The detectors were silver foils 
5 X 5 cm with a mass of 0.06 g/cm? and the cadmium 
screens had a mass of 0.35 g/cm?, so that their absorption 
in the thermal region was a total one. Activities were 
measured with an ionization chamber and electrometer. 

At room temperature we found 8=0.82 and the same 
result was obtained at liquid-air temperature. 

This result shows that the assumptions of isotropic 
scattering and of constancy of \ and 7 are not valid. There 
may be many reasons for this: first, as was shown by 
Fermi? for thermal neutrons the mean free path \ depends 
on the ratio W/hv of the energy of the neutrons to the 
quantum of the elastic bond of hydrogen atom in paraffin, 
in such a way that a decrease of energy gives rise to a 
decrease of A. A second fact, which must be considered, is 
the coherence of the directions of the motion of a thermal 
neutron before and after an impact, which was neglected 
in the deduction of (1). This coherence depends also on 
the ratio W/hy and is zero only for W<hyv. Both effects 
are in the sense of increasing the albedo for low tem- 
perature. 

Though it is now 
explanation of the observed fact on account of multiplicity 
of factors involved, still the experimental result is of 
interest in connection with the interpretation of experi- 


impossible to give a quantitative 


ments on thermal neutrons. 
E, AMALDI 
E. SEGRE 


Columbia University, 

August 29, 1936. 
!' E. Amaldi and E. Fermi, Ricerca Scient. [7] 1, 454 (1936) 
? E. Fermi, Ricerca Scient. [7] 2, (1936) 





~ 


STTERS TO 


nN 


Crystalline Symmetry and Shear Constants of Rochelle 
Salt 

The equations' for the converse piezoelectric effect in the 

hemihedric and hemimorphic classes of the orthorhombic 


system are: 


Hemihedric Hemimorphic 
Czz = Cyy =€,,=0; C22 =3;F:; 
Cyz =A Ez: Cyy = d32F;; 
€:z =d25E,; Cz: = G33; ; 1) 
Cry =AzcE:; €yz =du Ey; 
€zz = 414E;; 
é,, =U 


In the hemihedric class alternating electric fields applied 
along X, Y, or Z can excite, respectively, the yz, sx and xy 
shear modes of vibration.* The situation is quite different 
in the hemimorphic class. In it longitudinal modes whose 
major displacements are executed along X, Y or Z may be 
excited by E,, while the yz and zx shear modes are excited 
by EZ, and E,, respectively. Each of the above modes of 
vibration gives rise to a characteristic interference pattern® 
in the simple interferometer. Hence the above crystal 
classes are easily differentiated by observing how the 
amplitude of vibration of the shear or longitudinal modes 
varies with the direction of the applied field. 

Since there has existed some doubt‘ as to whether 
Rochelle salt is invariably hemihedric, it is desirable to 
check by independent means the conclusions of Staub‘ 
from x-ray data. The crystallographic d, b and ¢ axes of 
Rochelle salt have been identified with X, Y and Z, 
respectively. 

Twelve rectangular Rochelle specimens were studied in 
the interferometer at frequencies ranging from 20 to 200 
kc/sec. Their crystalline symmetry was found to be hemi- 
hedric with no transitions to the hemimorphic class in the 
temperature range from 0 to 40°C. These results confirm 
those of Staub. 

Representative frequency-temperature curves in Fig. 1 
show negative temperature coefficients for all three shear 
modes. All curves show anomalies near 23.5°C (compare 
with the curves of Davies. These anomalies may consist of 
discontinuities, points of inflection or the presence of 
maxima or minima. The presence of the maximum in the 
curves for the yz and xy shear modes of vibration depended 
upon the previous thermal and electrical treatment of the 
crystal. A more detailed study of this effect is being made. 

The shear constants C44, C55 and Css were evaluated at 
22°C from the measured frequencies of the yz, sx and xy 
shear modes, respectively, by means of the equation’ 

F vin = 43(Ci;/p) ¥(m2/x 32 +n? /x;*)4, (2) 
where C23 =¢44; Cis =¢s5; Ci2=Cee; and p was taken as 
1.7605 g/cm’. 

It was found, however, that (2) is far from accurate in 
Rochelle salt. The apparent elastic constant thus deter- 
mined varies parabolically when plotted against the ratio 
of x; to x;. The maximum value of the constant lies at the 
vertex of the parabola, where x; =x;. The maximum value 
of the shear constants so determined are compared below 
with Mandel’s® values obtained from static experiments. 
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Fic. 1. Variation of frequency of the shear modes with temperature 
Change in frequency in kc/sec. as ordinates; temperature in °C as 
abscissas. Frequency yz shear mode is 168.9 kc/sec. at 1°C, sx shear 
mode 28.3, xy shear mode 52.8. 


Writers Mandel 
448.72 X10" dynes/cm? 16.36 X 10" dynes/cm? 
55 2.46 3.24 
ces 7.95 12.43 


The discrepancy between these dynamically and 
statically determined values is so large that the writers 
plan to make a redetermination of the static values. 

A large coefficient of viscosity,’ b4;, was apparent in the 
yz shear mode. This was manifested by broad tuning in 
large plates of the resonant fundamental frequencies 
Thus at frequencies of 50 kc/sec. the amplitude of oscilla- 
tion would remain strong over a frequency range as large 
as 10 kc/sec. With the application of static electric fields 
the shift of the interference fringes showed a time constant 
of more than one minute. The viscosity coefficients 5 
and be. must be considerably smaller than b4,. Evidence 
indicated that Eq. (2) fails in Rochelle salt partly since it 
does not take into account viscosity. 

The writers wish to express their thanks to the Brus) 
Laboratories for gifts of both cut and uncut Rochelle salt 
crystals. 

R. TASCHEK 
H. OSTERBERG 
University of Wisconsin, 


August 27, 1936. 


1 Voigt, Lehrbuch der Kristallphysik, p. 830 

? Osterberg and Cookson, Rev. Sci. Inst. 6, 347 (1935 

3 Osterberg and Cookson, Rev. Sci. Inst. 6, 247 (1935). 

* Valasek, Phys. Rev. 20, 639 (1922). 

5 Staub, Helv. Phys. Acta 7, 1 (1934). 

* Davies, Phil. Mag. 16, 97 (1933). 

? Osterberg and Cookson, Phys. Rev. 47, 781 (1935). 

* Mandel, Proc. Roy. Soc. 116, 623 (1927) 

* Thompson, Phil. Trans. Roy. Soc. London 231, 365 (1933) 
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Some Theoretical Considerations on Anomalous Diffrac- 
tion Gratings 


I wish to publish some theoretical considerations about 
the phenomena shown by certain gratings in the experi- 
ments of Wood! and Strong.? Both consider as a “typical 
case” the following phenomenon exhibited by their grating 
C. Let us consider the intensity of a continuous spectrum 
as a function of the increasing wave-length. Then: (a) at 
certain characteristic wave-lengths (Rayleigh’s wave- 
lengths) the intensity falls nearly to zero and shows a sharp 
edge; (b) after having been nearly zero for an interval of the 
order of 100A (dark band), the intensity rises to abnormally 
high values and shows a bright band without sharp edges. 

Part (a) may be explained by Rayleigh-Voigt’s theory.® 
I wish to point out that the derivative of that second-order 
term, which carries the irregularity in the distribution of 
intensity, becomes infinite for Rayleigh’s wave-lengths. 
Hence a sharp edge in the intensity must be present in 
every grating, but only the geometrical characteristics 
determine whether the irregular term is the important one 
for the intensity, or not. 

In connection with part (b) of the phenomenon the 
following should be noted. (I) A suggestion of Langer 
related by Wood: Since the position of the bright bands 
depends upon the optical constants of the metal of the 
grating, these bands would be connected with that part 
of the diffracted waves which travels with different ve- 
locity along the surface of the grating. (II) Rayleigh’s 
condition for the phenomenon (a) means that the wave- 
lets with velocity c diffracted by a groove must reach the 
neighboring grooves in resonance with the incident wave; 
similarly when slower waves were diffracted by the grooves, 
one would expect effects of resonance for larger wave- 
lengths. (III) Applying the principle of Huygens to the 
diffraction, we must remember that an oscillator on a 
metallic surface sends also Sommerfeld’s superficial waves,* 
which, to be noted, are polarized with their magnetic 
vector parallel to the surface, as the anomalously diffracted 
light. 

Starting from these considerations I made an attempt to 
extend Rayleigh-Voigt’s: theory, including superficial 
waves on the grating. This extension is not quite formally 
satisfactory, but it seems to me to be sufficiently justified. 
It explains the essential features of the phenomenon (b); 
the width of the bright bands is due to the absorption of 
the superficial waves, and their distance from Rayleigh’s 
wave-length depends upon the strong variation of the 
optical constants with wave-length. Wood's result that a 
negative image of the bands is found in the spectrum of 
the regularly refracted light, is easily understood from the 
point of view of the approximation method. 

The theory was checked by a rough numerical evaluation 
starting from suitable data. In carrying out the calculations 
it appears clearly that the whole phenomenon of the bands 
is extremely sensitive to changes of the form of the grating. 
Hence one has to think that every anomalous grating will 
show its own type of bands. It follows also that the study of 


grating C was very useful in distinguishing phenomena (a) 
and (b) which seem to be due to quite different physical 
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mechanisms. I do not think, however, that it will be easy 
to define the characteristics of a grating showing the same 
bands as grating C. The presence of strong bands seems to 
be due to accidental relations between the coefficients of 
Fourier’s development of the profile of the grating; one 
can only say that the presence of narrow grooves or ridges 
in the profile, that is, a prevalence of high harmonics, is a 
necessary condition. A detailed account will be published 
in the Nuovo Cimento. 
Uco Fano 
Istituto Fisico della R. Universita 


Rome, Italy, 
July 21, 1936 


Wood, Phys. Rev. 48, 928 (1935 
2 Strong, Phys. Rev. 49, 291 (1936 
' Lord Rayleigh, Proc. Roy. Soc. 79, 399 (1907); Voigt, Gétt. Nachr 
41 (1910 
‘Sommerfeld, Ann. d. Physik 28, 665 (1909 


The Paschen-Back Effect—*S *P Multiplets 
in Strong Fields 


In an article published in the Physical Review on the 
Paschen-Back effect of 2S ?P multiplets Green and Loring! 
have stated that the fine structure of magnetic triplets in 
“strong fields” had never been experimentally established, 
the theory, however, having been expounded long before. 
The authors then explain how they have studied the 
phenomena for the 6708 line of Li “in reversal” and 
managed to resolve the s components. 

Presumably the authors had not had information about 
my experimental work on the same subject carried out at 
the Laboratoire du Grand Electro-aimant de |’ Académie 
des Sciences at Bellevue, under the direction of Professor 
A. Cotton. Nevertheless the results had been published as 
early as January 1935 in the Comptes Rendus Acad. Sci. 
and republished in the Jubilee Book of Peter Zeeman ('S 
Gravenhage, Martinus Nijhoff, 1935) in an article by F. 
Croze and P. Jacquinot. 

In this work the 6708 line has been studied up to 44,000 
gauss and I stated explicitly that the s components are 
resolved into two sub-components and that, moreover, the 
components on the short wave side give a larger separation 
than those on the other. The values of the separation | 
found are in satisfactory agreement with the values 
calculated after K. Darwin. The results had been obtained 
for the Li emission line given by a high frequency tube as 
described in the two notes previously quoted, and are 
quite the same as those obtained by Green and Loring in a 
38,000-gauss field with the line ‘in reversal.” 

It was with great pleasure that I had the opportunity of 
reading the confirmation of my own experiments in the 
article by Green and Loring whom I am always personally 
happy to quote. 

P. JACQUINOT 


August 31, 1936. 


! Green and Loring, Phys. Rev. 49, 630 (1936 
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The Nature of Energy States in Solids 


Using the methods of quantum mechanics and group 
theory Van Vleck,! Penney and Schlapp? and Amelia 
Frank* have developed a theory for the behavior of energy 
levels of an atom when subjected to the electric fields of 
other atoms in a crystalline solid. They have treated in 
detail the hydrated sulphates of neodymium, praseo- 
dymium and samarium. 

They applied their theory to the experimental results of 
Gorter and de Haas‘ on the magnetic susceptibility of Pr 
and Nd salts and were able to reproduce the experimental 
curves from 14°K to room temperatures. While the agree- 
ment here was very good, unfortunately, a number of other 
discrepancies made their appearance. Thus the potential 
field constant D in the equation V=D(x*+y'+<2*), which 
one would expect to be very similar for the various rare 
earths due to their similar chemical and crystallographic 
properties, is about four times as large for neodymium as 
for praseodymium. Further, while one would expect Er 
and Nd to behave somewhat similarly since they have the 
same spectroscopic term only inverted, the actual sus- 
ceptibility data are not alike, for the 1/x, 7 diagram for 
Er gives almost a straight line while for Nd it is distinctly 
curved. It should be mentioned, however, that the experi- 
mental work of Gorter and de Haas is not in agreement 
with two room temperature points determined by Zernicke 
and James,® and Cabrara.*® Finally, the agreement with 
samarium is extremely poor. Spedding and Bear’ have also 
found experimentally at least five low lying levels while 
the maximum number permitted by theory is two in a 
cubic field and three in a combined cubic and rhombic field. 

In view of these discrepancies I have thought it advisable 
to study the energy levels of the rare earths spectro- 
scopically. In collaboration with Dr. Nutting and Mr. 
Hamlin, I have found that the lowest state of neodymium 
is split into at least three levels with a separation of 76 and 
260 cm™ from the lowest state. Theory as applied to the 
work of Gorter and de Haas demands about 0, 243 and 
834. Our work is in definite disagreement with these results. 
However, if the experimental work of the above authors is 
disregarded the relative separations are the same as pre- 
dicted by theory. With our separations most of the dis- 
crepancies disappear. The magnetic susceptibilities as 
calculated from these levels give a 1/x, T diagram which 
is almost a straight line in good agreement with the experi- 
mental work on erbium. Further the calculated suscepti- 
bility agrees within experimental error with the room 
temperature values of Zernicke and James, and Cabrara. 
Finally the field constant D, being only one-third as big as 
before is in good agreement with that of praseodymium. 
The discrepancy in the case of samarium seems to be real 
and bears out the author’s contention that there must be 
an additional electronic level lying near the basic one in 
that case. 

In exactly the same way that Penney and Schlapp 
treated neodymium and praseodymium, I have treated 
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erbium and dysprosium using the field constants as deter- 
mined by Nd after correcting them for the additional 
screening according to Slater’s method and found the basic 
state of Er*** would be split into five levels separated by 
19 cm—, 38 cm™, 83 cm™ and 87 cm™ from the basic state: 

Nutting and Meehan have advised me by letter that 
they have observed levels at 0, 19, 41 and 85 cm™. This 
agreement is remarkable when one realizes that there are 
no arbitrary constants involved so that the relative 
splitting cannot be altered. As there is a slight leeway in 
the value of the screening constants used, the over-all 
spreading could be varied slightly and has been chosen to 
give the best fit. The value used, however, is well within 
the experimental error of determining screening constants 
by any other method. 

The levels calculated for Dy while of the right order of 
magnitude were not in good agreement with the experi 
ment, but this is not surprising for the basic state of Dy 
is the same as that for Sm, except that it is inverted, and 
if an extra level exists there one would expect it to be 
present in Dy also and while it would probably be farther 
away it would still perturb the splitting of the basic state. 

The excellent agreement in the case of Nd and Er seems 
to indicate that Penney and Schlapp’s assumption of a 
cubic field is justified. While the magnetic susceptibility 
would not be very sensitive to the existence of a rhombi 
term in the potential energy equation, the splitting of the 
levels would be extremely sensitive. From our data the 
rhombic term could not be more than 2 percent of the cubic 
term as there is absolutely no evidence of the rhombic 
splitting in our photographs. 

From the above results it can be seen that a study of the 
absorption spectra offers a powerful tool for determining 
the position, symmetry and distance of the atoms in the 
first row about the ion being studied. It, therefore, com- 
plements x-ray studies in that this method is most powerful 
in the region where the x-ray measurements are least 
reliable. It also should give considerable information about 
the potential and electric fields which exist about an ion 
in a crystal or in solution. 

I wish to thank Dr. Nutting and Mr. Meehan for very 
generously permitting me to use their values for erbium in 
advance of publication. I wish also to thank Professor 
Kirkwood for valuable assistance on the theoretical side 
of the problem. 

This work was carried out under the George Fisher 
Baker Research Fellowship at Cornell University. 


FRANK H. SPEDDING 


Cornell University, 
Ithaca, New York 
September 3, 1936. 


1Van Vleck, The Theory of Electric and Magnetic Susceptibility 
(Oxford Press). 

2? Penney and Schlapp, Phys. Rev. 41, 194 (1932). 

* Amelia Frank, Phys. Rev. 48, 765 (1935). 

* Gorter and de Haas, Leiden Comm. 218% 

* Zernicke and James, J. Am. Chem. Soc. 48, 2827 (1926). 

* Cabrara, Comptes rendus 180, 669 (1925). 

7 Spedding and Bear, Phys. Rev. 46, 975 (1934). 
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LETTERS TO 


On the Chemical Differences Between Nitrogen Isotopes 


We wish to report on some preliminary experiments 
concerning the possibility of a chemical separation of the 
nitrogen isotopes. Although our work in this field has not 
yet been completed, a publication by E. Ogawa' causes us 
to publish a short summary of our conclusions at this time. 

The analyses were performed with a mass spectrometer 
in the laboratory of Professor Bleakney in Princeton 
University as described by Wahl, Huffman and Hipple.* 
As these authors found indications of a small difference in 
vapor pressure between N'*H; and N!°H; and because 
distillation, if effective, would be by far the easiest way to 
obtain a change in the isotope rates, we first experimented 
on this possibility. We used the method by which Wahl and 
Urey investigated the vapor pressure ratios in water.’ 
Distillations were carried out without boiling at different 
temperatures, the final residue varying between 1/15 and 
1/70 of the amount of starting material. No definite trend 
of the vapor pressure ratio with temperature could be 
distinguished. As an average value for the pressure ratio 
P(N'4H;) : P(N!5H;) we obtained 1.0025 in six experiments, 
with an average error of 0.003. It may be mentioned that 
Wahl, Huffman and Hipple give 1.0052 as their best figure 
for this proportion. It is probable that N'4H; is slightly 
more volatile than is N'°H;, but we feel convinced that the 
vapor pressure ratio differs from unity only a few tenths of 
a percent at most. Therefore the nitrogen separation by the 
distillation of ammonia does not seem promising. Next we 
investigated the reaction of N®H;+N"“H,OH@N"'H 
+N“H,OH for which Wahl, Huffman and Hipple had 
reported a possible but uncertain concentration of the 
heavy isotope in the NH,OH. Purified air was passed 
through a bottle containing concentrated ammonia till 
only a small fraction of the original amount of ammonia 
was left in the solution. From the isotopic composition of 
the original material and the final sample we calculated the 
equilibrium constant from the Rayleigh formula. We found 
that the N'® possibly concentrates slightly in the liquid but 
the effect is so small at any rate, that it could not be 
recognized with certainty. 

Then we turned to the exchange reaction N'°H;+ N'“H,* 
=N'"H;+NH,", already considered by Urey and Greiff.* 
We started from a concentrated solution of ammonium 
chloride, to which at intervals we added small amounts of 
sodium hydroxide solution. A current of purified air 
removed the NH; formed after each addition. At the 
moment we have the analysis of only one run in which we 
started with 0.61 mole NH,Cl and finished with 0.052 mole. 
During this process the ratio of N'*N™ decreased from 
250.2 to 233.8. (These numbers are not corrected for diffuse 
background and hence the analysis of the original material 
does not agree with the generally accepted value). This 
gives for the fractionation factor the value 1.027. We hope 
to check this value soon with other samples we have 
prepared, but we have been forced to abandon this matter 
temporarily because it will not be possible to obtain 


analyses by the mass spectrometer for some time. 
The last experiment described seems to be the same as 
one performed by Ogawa, an extract of which publication is. 
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printed in Nature.’ From the figures given fractionation 
factors ranging from 1.7 to 1.8 can be calculated which are 
very much higher than ours. We feel quite certain that as 
large an increase in concentration of N'® could not remain 
undetected by the mass spectrometer and that Ogawa's 
values, which are obtained from atomic weight determi- 
nations, are too high. 

The importance of finding the most favorable reaction is 
clearly seen from the calculations by Urey and Greiff on 
the efficiency of a column as depending on the separating 
factor of the process used. 

Finally, it may be mentioned that a preliminary experi- 
ment on the exchange HCN'®+CN'*-@HCN'"+CN™ 
indicated that N'® was concentrated in the HCN gas to 
quite an appreciable extent. 

We are indebted to Dr. W. W. Lozier for the analyses 
made for this work, and also to the E. K. Adam's Fellow- 
ship Fund for,a stipend for one of us (A. H. W. A.). 

H. C. UREyY 
A. H. W. ATEN, JR. 
Department of Chemistry, 
Columbia University, 
New York, N. Y 


1 E. Ogawa, Bull. Chem. Soc. Japan 2, 428 (1936) 
? Wahl, Huffman and Hipple, J. Chem. Phys. 3, 434 (1935 
3 Wahl and Urey, J. Chem. Phys. 3, 411 (1935) 


4 Greiff, J. Am. Chem. Soc. 57, 321 (1935 
5 E. Ogawa, Nature 138, 294 (1936 


Effect of Hydrocarbons at Low Temperatures on Slow 
Neutrons 


Investigations of the properties of slow neutrons passed 
through paraffin at various temperatures have established 
the fact that temperature equilibrium must be attained at 
least partially.'-> The somewhat remarkable character of 
this fact, considering the size of the energy losses and the 
probable mechanisms involved led the authors to investi- 
gate the phenomenon in more detail, first with paraffin and 
later with pure normal heptane. 

The apparatus used is shown in Fig. 1. It consists of a 
metal Dewar (brass and monel metal walls 1 mm thick) 
surrounded by a large cylindrical block of paraffin in 
which the radon-beryllium neutron source (100-220 mil- 
licuries) was imbedded. The inside diameter of the Dewar 
was 14 cm and the height one meter. A soft glass Dewar 
(3 cm outside diameter) was mounted along the axis of the 
metal Dewar, to allow irradiation of the cylindrical silver 
targets without alteration of the temperature of the 
hydrocarbon. The cylindrical copper container shown was 
used for the heptane only, and was the same length as the 
silver targets (11 cm), and 4.5 cm thick, with a wall thick- 
ness of 0.3 mm. Copper-constantan thermocouples were 


Taste I. 
Thickness of Temperature Activity 
paraffin cylinder (°K) (% 300°K activity 
1.5 300° 200 100 +2 
100 11° 11342 
2.5 70 9442 
21° 8622 
3.5 00 103 +2 
55° 7142 
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soldered to the top and bottom of the container to measure 
the temperature of the hydrocarbon. The targets were 
silver cylinders 13 mm in diameter, 11 cm long, and 0.25 
mm thick, which could be slipped tightly over the thin 
wall tube counter used for measuring the activities. In the 
work with paraffin, the copper container for the heptane 
was removed and replaced by cylinders of paraffin of 
different thicknesses. 

The procedure was to cool the contents of the Dewar 
to 20°K by keeping liquid hydrogen around the container 
for approximately 20 minutes. The hydrogen was then 
allowed to boil away and readings were taken when the 
temperature of the hydrocarbon had risen to 25°K, thus 
insuring the absence of any liquid hydrogen in the Dewar. 
Readings of the silver activities were taken continuously, 
as the Dewar warmed, over a period of eighteen hours. 
The results for normal heptane are shown in Fig. 2 in which 
the circles show the measured activity (140 seconds period). 
The solid curve was calculated from the equation 


, 1 
nT?) 


ww? *(1—e 


A = be 
in which A is the silver activity, 5, y and 7 are constants, 
and T is the absolute temperature of the hydrocarbon. The 
assumptions on which the equation is based are (a) that 
the neutrons are at the temperature of the hydrocarbon, 
(b) that the proton capture cross section is proportional to 
the reciprocal of the neutron velocity, and (c) that the 
silver capture cross section obeys the same law. The con- 
stant 7 is evaluated from the thickness of the silver targets 
and the cross section, 36X10°** cm?, from the data of 
Amaldi and Fermi‘ for the capture by silver of the neutrons 
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exciting the 140 seconds silver activity. The data seem to 
be in essential agreement with these assumptions. 

The data obtained in some preliminary work with 
paraffin cylinders, are presented in Table I. 

The 1.5 cm cylinder caused an increase in the activity, 
which appeared between 200° and 100°, and was relatively 
constant from 100° down to liquid helium temperatures 
(4°K); the 2.5 cm cylinder showed a small decrease, slightly 
dependent on temperature at 70° and below; and the 3.5 
cm cylinder showed a rapid decrease between 90° and 60 
Though none of these results were obtained over as wide a 
range of temperatures as those for the heptane, they are 
all in agreement with the heptane results, and the assump 
tions listed above. 

All of the above results apply to the thermal neutrons 
(those absorbed by 0.3 mm Cd) which excite the 140 
seconds silver activity. A series of experiments with a 0.3 
mm thick cadmium shield around the silver targets showed 
that at least 92 percent of the activity disappeared on 
interposition of the cadmium shield, roughly independent 
of the temperature (20°-300°K) of the surrounding hydro- 
carbon cylinder. The activity with the cadmium shield 
interposed was too small to measure accurately. This large 
cadmium absorption effect may be connected with the 
fact that the majority of the neutrons pass through approx- 
imately 2 mm of soft glass and 2 mm of brass and monel 
metal before striking the target. 

The authors wish to express their gratitude to Professor 
G. N. Lewis for his invaluable advice and criticism, and to 
Mr. J. W. Stout for the preparation of several liters of 
liquid helium. 

W. F. Lipsy 
Ear. A. LoncG 


Department of Chemistry, 

University of California, 
Berkeley, California, 
September 3, 1936. 


1 Dunning, Pegram, Mitchell and Fink, Phys. Rev. 47, 888 (1935 
2? Moon and Tillman, Nature 135, 904 (1935). 
Dunning, Pegram, Fink, Mitchell and Segré, Phys. Rev. 48, 704 
1935). 
‘Westcott and Niewodniczanski, Proc. Camb. Phil. Soc. 31, 
1935). 
Moon and Tillman, Proc. Roy. Soc. A153, 476 (1936) 
* Amaldi and Fermi, Ricerca Scient. [VI], 2, 344 (1935) 
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Bond Effect in the Action of Protons on Neutrons 


The marked ability of protons bound in hydrocarbons 
to slow neutrons down to thermal equilibrium with the 
hydrocarbon, as indicated in the preceding letter and 
earlier work'~* induced the authors to undertake a syste- 
matic investigation of different types of hydrogenous 
molecules to discover any dependence of this phenomenon 
upon the character of the molecule and the bond. The 
experiments with liquid hydrogen here described reveal 
the existence of a remarkable difference between the H,. 
molecule and the hydrocarbons. 

The apparatus used is shown in Fig. 1 of the preceding 
communication with the exception that the copper con- 
tainer used for the heptane was removed, leaving the space 
between the inner soft glass Dewar and the outer metal 
Dewar empty. The procedure used was to determine the 
activity of the silver targets (1) with the Dewar empty, 
(2) with the Dewar filled with liquid or solid hydrogen to 
a level approximately 5 cm above the top of the silver 
target, (3) same as the preceding but with the liquid 
hydrogen nearest the inner Dewar replaced by a 1.5 cm 
cylinder of paraffin, and (4) at room temperature, the 
Dewar being filled with a solution of methyl alcohol in 
acetic acid (21 percent alcohol by volume) which contained 
the same number of protons per cubic centimeter as does 
liquid hydrogen. Table I presents the results corrected 
for radon decay. The solid hydrogen was obtained by 
exhausting the Dewar until an equilibrium vapor pressure 
of 3.8 cm of mercury was obtained. 

The comparison of Experiments 2a—) either with 
Experiment 3 or with the experiments described in the 
preceding communication shows the great difference 
between liquid or solid hydrogen and a hydrocarbon. 
A comparison of Experiments 2 and 4 shows a nearly 
identical activity whether the Dewar is filled with hydrogen 
at very low temperatures or with a solution of the same 
proton density at room temperature. Moreover, the same 
result had been obtained in a preliminary experiment in a 
smaller Dewar in which the thickness of the liquid hydrogen 
layer around the soft glass Dewar was 3.5 cm instead of 
5.5 cm. It would not be correct to conclude, however, 
that the liquid and solid hydrogen do not cool at all. 
The experiments indicate, rather, that the low speed 
neutrons have been cooled by liquid hydrogen at 20°K, 
and by solid hydrogen at 12°K, to a mean temperature of 
perhaps 150°K and that the increase in silver activity 
which should be thus produced has been nearly offset by 
the increase in capture by protons. This view is in agree- 
ment with the results of Fomin, Houtermans, Leipunsky 
and Schubnikow,® on the backward scattering of slow 
neutrons from liquid hydrogen and from water at room 
temperature. 

In interpreting Experiment 3 we assume, therefore, that 
the 300°K neutrons which enter the Dewar have been 
partly cooled by the liquid hydrogen and then brought to 
very low velocities and consequently absorbed by the 
thin layer of paraffin. It must in any case be concluded 
that at very low temperatures paraffin has a much greater 


effect either in decelerating or absorbing these neutrons 
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TABLE I 

: Activity of Ag 
Experiment Description 140 seconds 

1 Empty Dewar (300°K 593+15 

2a Liquid hydrogen (20°K 769 +15 

b Solid hydrogen (12°K 7582415 

3 1.5 cm paraffin cylinder +liquid 579+15 

20°K.) hydrogen 
4 Acetic acid solution (300°K) 776415 


than cold hydrogen. Unless this difference can be ascribed 
to the carbon nucleus itself, and experiments with which 
we are now engaged should soon settle this question, the 
difference in behavior between very cold hydrogen and 
paraffin must be ascribed to the chemical bond. 

Our experiments are analogous to those of Westcott and 
Niewodniczanski* on paraffin and hydrogen at 20°K. 
They performed two sets of experiments and found it 
impossible to reconcile their first set with their second. 
Our results agree closely with their first set. 

The authors again wish to thank Professor G. N. Lewis 
for his interest and assistance in this work. They are also 
indebted to Mr. Samuel Ruben for assistance with the 
experimental work, and to Dr. L. R. Taussig and the 
University of California Hospital for donation of the 
radon used. 

W. FF. Lippy 


EARL A. LONG 
Department of Chemistry, 
University of California, 
Berkeley, California, 
September 3, 1936. 
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(1935). 
Moon and Tillman, Proc. Roy. Soc. A153, 476 (1936 
* Fomin, Houtermans, Leipunsky and Schubnikow, Physik. Zeits 
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Resonance Energy of Cadmium for Neutron Capture 


The experiments of Rasetti, Segré, Fink, Dunning and 
Pegram! on the absorption of slow (thermal) neutrons by a 
layer of cadmium moving with a velocity relative to the 
neutron source have indicated that the collision cross 
section remains practically constant in an energy interval of 
about 20 percent of 3k7/2. On the other hand Powers, 
Fink and Pegram? have reported a 7 percent increase in the 
absorption coefficient on lowering the neutron temperature 
from 300°K to 95°K. 

It is of interest to consider these two sets of experimental 
data in terms of the Breit-Wigner formula for the cross 
section here involved. For this purpose the latter formula 
can be written in the form: 


o =1/E'|(E+1)?+B*}, 


where E and B are measured in energy units Eo; and @, in 
arbitrary units. Eo represents the resonance energy; and B, 
the half-width of the virtual level responsible for the cap- 
ture. The plus sign refers to the case of positive resonance 
energy; and the minus sign, to the case of negative reso- 
nance energy. 
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For positive resonance energy it is easily seen that when 
B<2/¥v5 the cross section has both a minimum and a 
maximum value; and when B=2/¥y75, these two extremal 
values coalesce at a horizontal flex, for which E=3/5. For 
still broader resonance o decreases steadily with increasing 
energy, without showing a resonance maximum. These 
variations in shape of the cross section curve are very 
clearly illustrated by a set of curves in o and E for given 
values of B. 

The failure in the first of the above-mentioned two sets 
of experiments to detect an appreciable variation of « with 
E, throughout a 20 percent range around 3k7/2, could be 
explained qualitatively, either if one of the extremal values 
of o occurred at E=3kT/2Eo, approximately, or if the point 
E=3kT/2E, lay in the long flat portion of a curve with a 
horizontal flex. 

It is, of course, not accurate to regard the thermal neu- 
trons as a homogeneous group in energy, but one should 
rather consider them as having approximately a Maxwellian 
distribution, and interpret the experimental results in terms 
of appropriate integrals. At the present stage of experi- 
mental uncertainty it does not seem worth while to carry 
out such integrations; but by considering the neutrons as 
homogeneous one should be able to determine the order of 
magnitude of the resonance energy Eo and the half-width 
BE). Each of the two sets of experiments permits an esti- 
mate of B asa double valued function of 3k7/2E 9; and both 
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taken together, an estimate of Ey and B. Thus, taking 
3kT/2=0.036 volt, the two relations lead to a resonance 
energy of 0.05; volt and a half-width of about the same 
value. In view of the rough approximation in the method 
of interpretation assuming homogeneous thermal neutrons 
and of the wide limits of relative experimental error (due 
to the great difficulties connected with experiments of the 
type here discussed, especially the thermal experiments 
the above values should be considered as giving little more 
than the order of magnitude of the quantities involved. 

The possibility of a negative resonance energy would, of 
course, be excluded by the first set of experiments, since ¢ 
would have no extremal values for positive neutron ener- 
gies, but would decrease more rapidly than 1/y¥ E. 

Before making a more careful analysis of such data it 
would be desirable to have data similar to the above men- 
tioned, in the first case with the neutron source at low 
temperatures, and in the second case for different tempera- 
ture intervals. 

E. SEGRE 
C. T. ZAHN 
Department of Physics, 
University of Michigan, 


Ann Arbor, Michigan, 
September 3, 1936 


! Rasetti, Segré, Fink, Dunning and Pegram, Phys. Rev. 49, 104 
(1936 
2? Powers, Fink and Pegram, Phys. Rev. 49, 650 (1936); see also Till- 


man, Proc. Phys. Soc. 48, 642 (1936 











